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Abstract 

Using the (nearly) Ricci-flat Swiss-Cheese metric of pQ, in the context of a mobile space-time 
filling Z?3-brane restricted to a nearly special Lagrangian sub- manifold (in the large volume limit, 
the pull-back of the Kahler form close to zero and the real part of the pull-back of e~ l8 , 8 = -| times 
the nowhere- vanishing holomorphic three-form providing the volume form on the three-cycle) of the 
"big" divisor with (fluxed stacks of) space-time filling D7-branes also wrapping the "big" divisor 
(corresponding to a local minimum), we provide an explicit identification of the electron and the 
u-quark, as well as their S'[/(2)L-singlet cousins, with fermionic super-partners of four Wilson line 
moduli; their superpartners turn out to be very heavy, the Higgsino-mass parameter turns out to 
be large, one obtains one light (with a mass of 125 GeV) and one heavy Higgs and the gluino is long 
lived (from a collider point of view) providing a possible realization of "/x-Split Supersymmetry" . By 
explicitly calculating the lifetimes of decays of the co-NLSPs - the first generation squark/slcpton 
and a neutralino - to the LSP - the gravitino - as well as gravitino decays, we verify that BBN 
constraints relevant to the former as well as the requirement of the latter to be (more than) the age 
of the universe, are satisfied. For the purpose of calculation of the gravitino relic density in terms 
of the neutralino/slepton relic density, we evaluate the latter by evaluating the netralino/slepton 
(co-)annihilation cross sections and hence show that the former satisfies the requirement for a Dark 
Matter candidate. 
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1 Introduction 



It is challenging to provide a suitable choice of vacuum (local minimum) of string theory to infer 
various cosmological and phenomenological issues, and other fundamental physics. One of the most 
exciting unresolved issues in particle physics and cosmology is the nature of dark matter (DM) in the 
Universe. Over the years, astronomical and cosmological observations have put significant constraints 
on its expected properties. The recent Wilkinson Microwave Anisotropy Probe ( WMAP) observations 
provide the relic abundance of a cold dark matter (CDM) [2] to be H-cDMh 2 = 0.1109 ± 0.0056. 
Theoretical status of DM generally emerges in the context of theories beyond Standard Model. It is 
well known that supersymmetric models with conserved R-parity contain one stable neutralino which 
is a candidate for cold dark matter. However, in models coupled to gravity and even in other scenarios 
of supersymmetry breaking mediation including gauge mediation, the gravitino ( the supersymmetric 
partner of Graviton) stands out as probably the most natural and attractive candidate for the LSP 
and DM although it generically suffers from a well-known "cosmological gravitino problem" [3] the 
resolution of which depends on whether one considers gravitino to be LSP or unstable particle. The 
long-lived gravitino's are generated in the early Universe. In the standard big-bang cosmology, they 
were in thermal equilibrium and then because of their weak gravitational interactions, frozen out 
while they were relativistic. In this case their abundance might overdose the Universe. Even if 
during inflationary epoch of the universe, primordial abundance of gravitino is completely diluted, the 
problem can not be solved because gravitinos are regenerated in the thermal bath after the reheating 
if the reheating temperature is high enough though overclosure of universe by dark matter constrain 
the exact value of reheating temperature [H [6]. However, in addition to it, the production of 
gravitino also depends on non- thermal production mechanism, the abundance of which is independent 
of the reheating temperature, i.e fioj^ = ^3/2 + ^3/2 SP ' Generically in case of heavy gravitino, re- 
heating temperature is low enough to produce appropriate abundance of Gravitino. Therefore sufficient 
number of gravitinos can be produced after NLSP decays to gravitino after its decoupling from the 
thermal plasma [7j. In other words, annihilation density does not depend on freeze out of gravitino 
(which causes a gravitino problem), instead depends on freeze out of NLSP. Hence relic density of 
gravitino is given by relic density of NLSP [71 [8] according to the relation = ^3/2- Therefore 
it is interesting to study the thermal cross-section and hence relic density of the (co-)NLSP(s) to get 
the the right order of relic density of gravitino. 

In spirit of above discussion, the purpose of present work is to give the signatures of gravitino as 
a potential dark matter in the context of the large volume limit of type IIB ("big divisor") D3/D7 
Swiss cheese phenomenology. The consistent compactification scheme and phenomenological features 
of the model were initiated in [9j [10] and "/i split-SUSY scenario" -like realization of the same was 
shown to be possible to be realized in [11] where in addition to getting one light Higgs (one of 
the primary motivation of (i split SUSY Scenario), we evaluated the life time of gluino which came 
out to be long and hence satisfied one of the important phenomenological features of \x split-SUSY 
scenario. In general, stability of LSP is governed by conserved R-parity. However, in this paper, we 
have considered cosmology of the gravitino LSP scenario based on R-parity violation. We argue that 
in addition to non-zero R-parity violating couplings, high squark masses helps to reduce the decay 
width and hence the lifetime of the gravitino decays become very long, typically of or larger than the 
age of the Universe. The explicit calculation of matrix amplitudes and hence life times of N(LSP) 
candidates requires the complete identification of the SM particles. Therefore we build up a set up 
which is able to provide an explicit identification of the first generation SU{2)l leptonic and quark 
doublets as well as their S'?7(2)i-singlets, with fermionic super-partners of four Wilson line moduli. 
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Further, soft SUSY parameters obtained in the context of gravity mediation provides sleptons and 
gaugino/lightest neutralino (being almost degenerate in mass) as co-NLSPs while gravitino naturally 
appears as Lightest Supersymmetric particle (LSP) which traditionally can be considered as viable 
dark matter candidate. The more explicit realization of the same requires life time of the LSP to be 
around and preferably more than the age of the universe. 

The organization of the paper is as follows: In section 2, we start off with details of an improved 
version of the large volume scenario set up discussed in |9j. Here, we construct four harmonic distri- 
bution one-forms supported on a sub-locus of big divisor localized along the mobile space-time filling 
Z)3-brane, and by utilizing the geometric Kahler potential of pQ we show that the aforementioned 
sub-locus in the large volume limit is nearly a special Lagrangian sub-manifold of the "big" divisor. 
By calculating the intersection matrix valued in the Wilson line moduli sub-space y 

appearing in Kahler coordinate Tg, we write (the Kahler sector of the) Kahler potential that includes 
four Wilson line moduli 01,02,03,04 and two position moduli of a mobile space-time filling D3-brane 
restricted to the abovementioned (nearly) SLag (corresponding to a local minimum); for the purposes 
of evaluation of "physical" /normalized Yukawa couplings, soft supersymemtry breaking parameters, 
etc., this Kahler potential is diagonalized to produce the Ai, -Zj-basis. Further, the estimate of the 
Dirac mass terms appearing in the N = 1 supergravity action of [12] calculated from superpotential 
and Kahler potential suggest that the fermionic superpartners of A\ and ^3 correspond respectively 
to the first generation leptons: and en, and the fermionic superpartners of A2 and A4 correspond 
respectively to the first generation quarks: ul and ur. We also provide explicit bi-fundamental rep- 
resentations for the four Wilson line moduli and two L>3-brane position moduli (after having turned 
on appropriate two-form fluxes on the D7-brane world volume decomposing adjoint- valued matter 
fields into bi-fundamental matter fields) and speculate about the possible modification in the relevant 
M = 1 chiral coordinate in the presence of the same. We also show that the (effective) physical 
Yukawas change only by 0(1) under an RG-flow from the string to the EW scale. Building up on 
this identification , we have evaluated the contribution of various three-point vertices in the context 
of gauged N = 1 supergravity action, in section 3. In 3.1, we study decay width and life time of 
two-body and hadronic three-body decays of gauginos (NLSP), life times of which, except for gluino- 
to-Goldstino decay, ensure that energy release from both gauge boson and hadronic decays do not 
spoil the the predictions arising from Big Bang Nucleosynthesis (BBN) and cosmic microwave back- 
ground, etc.; gluino is hence ruled out as an NLSP. In 3.2-3.3, we consider R-parity violating decays 
of the co-NLSPs - (Wino/Bino-dominated)neutralino and sleptons/squarks to ordinary particles (even 
R-parity) which are harmless because they rapidly decay to ordinary particles and do not disturb 
the beautiful predictions of BBN. In section 4, in order to meet the requirement of an appropriate 
DM candidate, we calculate life time of three-body R-parity violating decays of gravitino LSP which 
comes out to be more than the age of the Universe. However, more exact determination requires that 
relic density of dark matter must be compatible with present observations and should not overdose 
the Universe. Therefore, in order to be able to perform a reliable comparison between theoretical 
predictions and improving measurements of the relic abundance from underground DM searches, in 
section 5, we calculate the relic density of co-NLSP's and hence relic abundance of gravitino. Strictly 
speaking, we consider very conservative approach where one assumes that decays of NLSP's account 
for almost all of the dark matter of the Universe . Utilizing extensively the results of [13j . |14| for 
providing the exact analytical expressions of various annihilation channels of neutralino and slepton, 
we obtain estimates of thermal cross-section and hence relic density of N(LSP) candidates. Section 
6 has the conclusions. There are two Appendices; appendix A has some details of the goemetric 
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Kahler potential calculation and appendix B has details of calculations of various soft supersymmetry 
breaking parameters and the fact that its possible to obtain around 10 2 GeV masses for W/Z vector 
gauge bosons and 125 GeV for the light Higgs. 



2 Setup 

From Sen's orientifold-limit-of-F-theory point of view corresponding to type IIB compactified on a 
Calabi-Yau three fold Z-orientifold with 03/07 planes, one requires an elliptically fibered Calabi- 
Yau four-fold X4 (with projection tt) over a 3-fold B 3 (= CY3— orientifold) where B 3 is taken be an 
n-twisted CP 1 -fibration over CP 2 such that pull-back of the divisors in CI3 automatically satisfy 
Witten's unit-arithmetic genus condition. For n = 6 [15], the Calabi-Yau three-fold Z then turns out 
to be a unique Swiss-Cheese Calabi Yau in WCPj^ 1 1 6 g i [x\ : X2 : £3 : X4 : X5] given by a smooth 
degree-18 hypersurface in WCP 4 [1, 1, 1, 6, 9]; the exceptional divisor corresponding to resolution of a 
Z3-singularity x\ = X2 = X3 = 0, via the monomial-divisor map, is encoded as the (pXiX^x^cf) being 
an involutively odd complex structure modulus) polynomial deformation in the defining hypersurface. 
Also, throughout we will be working in the non-singular coordinate patch xi = 1 in the large volume 
limit of the Swiss-Cheese Calabi-Yau. 

In the presence of a space-time filling D3-brane and a space-time filling _D7-brane wrapping E#, 
the Kahler potential relevant to all the calculations in this paper (without being careful about 0(1) 
constant factors): 



K ' 2ln \[a B + Cr B - 7-^geom) 2 - (ffif + &S ~ 7-^geom) 2 ] 



3 , . 3 



+| { l~ T) ~ 2 4 E n ° E ( f~ r) " cos(mk.B + nk.c) 

m,nezV(o,o) (2ip\m + riT\3 peH - {C Y a ,x) m,neZV(o,o) ( 2i ) 3 \ m + nr l 3 

(1) 

will be rewritten in terms of the N = 1 chiral coordinates: T a= s ,5, 
definitions of all) - in particular: 

<7« ~ T a — (iK abc c b B c + K a + - - , K abc g b (g c - g c ) 

- 3? - i - - \ 

+i5^K\ l i 1 l 2 C I a J a I aj + jSarQf + ifi 3 l 2 (oj a ) i]Z l {z> - -z a {f>a)\z l )\ , (2) 

where 

• K4 is related to four-dimensional Newton's constant, \jl 3 and /X7 are D3 and L>7-brane tensions; 

• K aab s are triple intersection integers of the CY orientifold; 

• c a and b a are coefficients of RR and NS-NS two forms expanded in odd basis of H^_(CY); 

• C J a J = / sB i*uj a A A 1 A A J , u a e H^f(CY 3 ) and A 1 forming a basis for H [ ^\z B ) - immersion 

map is defined as: i : S s CI3, aj is defined via a Kaluza-Klein reduction of the U(l) gauge 
field (one-form) A(x,y) = A l ,(x)dx fl P-(y)+a I (x)A 1 (y)+aj(x)A J (y), where P-{y) = 1 if y € T, B 
and -1 if y G a(S B ); 



3 



K = In 



• z a ,a = 1, h 2 ! (CYs), are D = 4 complex structure deformations of the CY orientifold, (Vs,)] = 
pjp^*' (Xa) M j, i.e., V : rcy 3 (1,0) — ► rCy 3 (0,1) via the transformation: z< c - s -^ orm g i + 
\z a ('Pa)] z 3 1 z l are scalar fields corresponding to geometric fluctuations of D3-brane inside the 
Calabi-Yau and denned via: z(x) = z l (x)di + z l (x)di, and 

• Qf = 12 J S fl / A /, where / G ff£ = coker (h 2 (CY 3 ) ^ H 2 _{^ B 

The closed string moduli-dependent Kahler potential, includes perturbative (using [IT]) and non- 
perturbative (using [18J) a'-corrections. Written out in (discrete subgroup of) SL(2, Z) (expected to 
survive orientifolding)-covariant form, the perturbative corrections are proportional to x(CY^) and 
non-perturbative a' corrections are weighted by {n®}, the genus-zero Gopakumar-Vafa invariants that 
count the number of genus-zero rational curves (3 6 H 2 (CYs,Zi). In fact, the closed string moduli- 
dependent contributions are dominated by the genus-zero Gopakumar-Vafa invariants which using 
Castelnuovo's theory of moduli spaces can be shown to be extremely large for compact projective 
varieties [?] such as the one used. 

Based on the study initiated in |10|, 111), one us (AM), inspired by the Donaldson's algorithm - see 
[19J - had obtained in pQ, an estimate of a nearly Ricci-flat Swiss Cheese metric for points finitely 
separated from based on the following ansatz: 

h Z '~ Z h 4 2 Z 2 4 + h Z ^ Wz A Z A + ^My23/36 (zi +2l + Z2 + + ftfZyll/lBfa^ + ^ + ^ + ^ 



+ h z 4 z 4 yll/18 ^2 + ^ + + Z2 2 + z 2 + z x z 2 ) + h Z ^V 7/12 (z 2 Zi 2 + z\z X + z\z X + Z X Z 2 2 + Z 2 z\ + z\z 2 ) 
+h z*zl v 5/9 ^2-2 + z 2-2 + ^2-2 + ^2-2) + h 44jy^ + ^) + h^V 17 ^ [z X Z A + Z 2 Z A + Z X Z A + Z 2 Z A ) 

+h 44y 4 / 9 (z A z x 2 + z\z A + z\z A + z 2 2 z A ) + h z ^V 5/12 {z x z A z x 2 + z 2 z A z x 2 + z\z A z x + z 2 z\z A + z x z 2 2 z A + z 2 2 z 2 z A ) 

+ h Z ^V 5/1S (z 1 Z 1 Z A Z A + Z X Z 2 Z A Z A + Z X Z 2 Z A Z A + 22^2-24^4) + h Z ^V 11/36 ((z X + z x ) ZA Z A + (z 2 + Z 2 )z A Z A ) 

+h^l W (Z A 2 + Z 2 ) + ( 24 2 fl + Zi - Z 2 + Z2 -2 + ^^2) + ^^5/18 (^2 + -2^2 + ^2-2 + ^2-2) 

W (z A Z A 2 + Z 2 Z A ) + h Z ' S W 5 ^ 6 ((z X + Z 2 )Z A Z A 2 + fa + Z 2 fe) + /l^V^^ + z lZ ^ 



+ \zi\ 2 (z A + z 4 ) + z x z 2 (z A + z A ) + \z 2 \ 2 (z A + Z A )) + 



(3) 



From (lA3p . one sees that /r ~ h z * Zi V-s. Using (j25h . one can show that: 

***** ~ 2 , — „,2 - ( 4 ) 

' 3 



Solving numerically: J2 n=0 a n yi z i z -ij V z = 0, as was assumed, one (of the eight values of) h z ^, up 

— — r 

to a trivial Kahler transformation, turns out to be V » , V ~ 10 . Curiously, from GLSM-based analy- 
sis, we had seen in [S] that on T,b(z a = 0), the argument of the logarithm received the most dominant 

1 2-2 

contribution from the Fl-parameter r 2 ~ Vs which using this value of /i z 4 2 4 ; is in fact, precisely what 
one obtains even now. Using this value of h 2 * 2 *, one obtains: -R 2i z 4 > Rz 4 z 4 ~ 10 . Further, as has 
been assumed that the metric components g Zl 2 z 4 are negligible as compared to g Zi z, - this was used in 
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[lOj in showing the completeness of the basis spanning Hj for a large volume holomorphic isomet- 
ric involution restricted to z\ — > —Zi-,z 2i3 — >■ 22,3, which is now born out explicitly, wherein the 
latter turn out to about 10% of the former. One does not need the aforementioned restriction on the 
geometric Swiss-Cheese metric for large volume involutions restricted to ©, of the following type. In 
the u 2 7^ 0-coordinate patch, the defining degree-18 hypersurface in WCP| 1 1 6 9 [u\ : u 2 ■ u 3 : u 4 : it 5 ] 
can be written as: 1 + z{ 8 + zf + *f - i>z\zwi - $z\z\ = 0,zx = ^,z 2 = ^,z 3 = ^,z 4 = j^. 

The same can be rewritten as: {iz 4 ) 2 — tyz\z 2 z 3 z 4 = z| + 4>z\z\ — z^ 8 — z^ 8 - This is can therefore be 

T 2 (z 3 ,z 4 {z 3 )) -> WCPi >1X6i9 ( Zl ,z 2 ,z 3 ) 
thought of as the following Weierstrass variety: 1 tt . Now, 

CP 2 (zi,z 2 ) 

assuming that the complex structure modulus ip has been stabilized to an infinitesimal value so that 
one can disregard the polynomial deformation proportional to ip, and defining xi 2 = ^i2> this elliptic 

T 2 (z 3 ,z 4 (z 3 )) -> WCPf )ljl)6 ' 9 (xi,X2^3) 

4 7T 

fibration structure can thought also of as: (Complex curve) C(xi) — * CP 2 (xi,X2) , 

Itt' 
CP 1 ( X2 ) 

corresponding to (iz 4 ) 2 ~ z| + </>XiX2 — Xi — X 2 - Alternatively, one can also think of the Weier- 

T 2 (- Xl ,z 4 (-xi)) -»• WCPf il)li6)9 (xi,X2^3) 

4 7T 

str ass variety as: (Complex curve) C{— z 3 ) — > CP 2 (— 2:3, X2) 1 corresponding to (iz 4 ) 2 « 

Itt' 
CP 1 ^) 

(— Xi) 3 + <^XiX2 + z i — xi- Therefore, near ([5]), one can define a large volume holomorphic involution: 
Xi O — 23, — >■ — 4>(the fact that complex structure modulus <p is involutively odd, is also used in (pj). 

After Zj — )• + V^e and aj — > aj + (a/), the EW symmetry is spontaneously broken to U(l) em . 
We therefore require to construct four harmonic distribution one-forms Aj localized on a sub-locus of 
S B : 

C 3 : \z\\ ~ Vse, |^ 2 | ~ Vze , 1 23 1 ~ Vs. (5) 

This is a toroidal three-cycle and the Calabi-Yau can be thought as a T 3 (swept out by 
(argzi, argz 2 , ar gz 3 ) )-fibration over a large base (|zi|, |z 2 |, | ^3 1 ) ; precisely apt for application of mirror 
symmetry as three T-dualities a la S(trominger) Y(au) Z(aslow) [20]. Interestingly, ([5]) is almost 
a s(pecial) Lag(rangian) sub-manifold. The requirements for the same are that f*J = 0, f*il = 
e l9 vol(C 3 ), where / : C3 — > CY 3 [21] . Let us see if these requirements hold up. Using the geometric 
Kahler potential of (|25p . the geometric metric near C 3 is estimated to be: 



0(10- 



/ -J— 1 1 \ 

/ If "W v 7 / 36 \ 

1 1 1 

-T^f -T^y yTpm 



1 1 



(6) 



which justifies the assumption made in [10] about the block-diagonal form of g^. The same was 
necessary in explicitly demonstrating h 1 ^ 1 7^ in [10] for the large volume involution: z\ — > —Z\, z 2j3 — > 
212,3- The most non-trivial example of involutions which are meaningful only at large volumes is mirror 
symmetry implemented as three T-dualities in |20] (Also see the first paragraph of section 4.3 of [22J 
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for implementation of SYZ proposal and importance of large base for doing so); mirror symmetry 
could be thought of as an involution since the mirror of a mirror is the original manifold. 



1. Using one sees that: 



/*(ds§)~0.05 



V is (dzidzj) +V a\dz 3 



0.05 



J, d(argz i )d(argz- j ) + d(argz 3 )d(argz 3 



(7) 



This implies /* J ~ 0.05, which is small. 



2. In the u 2 / 0- coordinate patch with z\ = z 2 = ^ 5 z 3 = £4 = by the Griffiths residue 



formula, one obtains: 



«5 

•A' 



n 



dz\ A dz2 A (iz4 ctei A dz2 A d&i 



(8) 



The coordinate £4 can be solved for using: 



l\)Z X Z 2 Z 3 ± ^^ 2 {z lZ2 Zz) 2 - 4(1 + Z \* + Zf + if - 0(Z!Z 2 ) 6 ) 2 1 



(9) 



which implies that the T 3 will never degenerate as the roots will never coincide: 

{2^z\z\z 3 - \2z\) 



dzi ~ dz 3 ( — ^— - ± 



2 4^( Zl z 2 z 3 ) 2 - 4(1 + 4 s + zf + 4- <f>( Zl z 2 ) 6 ) 



+ (...)dzi + (...)dz 2 



± 1.5i 



^1+^2 + 



dz 3 + (...)dzi + (...)d2 2 , near ©, |^| << 1, \(f>\ « 1. 



(10) 



Therefore, restricted to (|5|) and substituting Q and dTO]) in (jHJ), one obtains: 



18 ± 1.5«Vl2 



3V3 - 



2 1 



dzi A d£ 2 A dz 3 



tbV 18 _i\ 

± 0.5iV 4 d^i A d_z 2 A dz 3 . 



(11) 



Hence, if one could estimate the pull-back of the nowhere vanishing holomorphic three-form as: 



f*Q ~ ^ty^Ts + TS+T, + 0.5iV~k + T8 + h J d(argzi) A d(argz 2 ) A d{argz 3 ) 
~ (0.2^ ± 0.3z) d{argz\) A d{argz 2 ) A d(argz 3 ), 



V~10 6 



(12) 
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and: 

vol(C 3 ) ~ <ff*gf*{dzx A dz 2 A dz 3 ) 



~ V 7 (0.05) 3 V is + e d(argzi) A d(argz 2 ) A d(argz 3 ) 

v~io 6 

~ 0.2d(arg.Zi) A d(argz 2 ) A d(argz 3 ), (13) 

then relative to a phase e~ %e ,6 = |, 9m (/*e _ie fi) ~ Q.2ij)d{argz\) A d(argz 2 ) A d(argz 3 ), which for 
1^1 << 1, is close to zero; (equivalently)Re ( y f*e~ l9 Q) ~ vol(C3). This implies that C3 is nearly /almost 
a special Lagrangian sub-manifold [23]. 

Following our previous constructions in [9J [10] , the harmonic distribution one- forms can be con- 
structed by integrating: 

dM = (Ps fl (21,2,3))' cbi A dz2, (14) 

with 7 = 1, 2 (done), 3, 4 where 

A/ ~ 6 (\z 3 \ - 5 (\zi\ - 5 (\z2\ - V^) [coi(z 1 ,z 2 )dz 1 + w/(zi, z 2 )dz 2 ] . (15) 



From (I14p . one sees that ^4/ is harmonic only on Eg and not at any other generic locus in the Calabi- 
Yau manifold; ([15|) shows that Aj are distribution one- forms on Eg localized along the L>3-brane 
which is localized on the three-cycle C 3 of ([5]). 

Writing Ai{z\,z 2 ) = ui{z\,z 2 )dz\ + u>i(zi, z 2 )dz 2 where ui(-zi,z 2 ) = u(zi, z 2 ), u>(—zi, z 2 ) = 
—Cj(zi,z 2 ) and d\ui = —d 2 ui, one obtains: 

Ai(zi, Z 2 , Z 3 ~ V*) ~ -zf^^l + 4 9 4 8 ^2, 

1 / z 19 \ ( z 19 \ 

A 2 (z u z 2 ,z 3 ~ Ve) ~ - ( +^ 8 2 2 J dzi + ( +z 2 18 zi j dz 2 , 

*S, * ~ V*) ~ - 3^ (VV + zf) +z 2 (VV + zf) \ Z -l(3V + 6VVz{* + 3zf) ) dz x + (1 2) 
zfzf dzi + (1 O 2), 

^ 4 (zi,2 2 ,z 3 ~ V5) ~ -zfzf dz x + z| 6 zf dz 2 . (16) 
Hence, 

Cu ~ y J z\*zfzfzf8 (\zi\ - V^) <5 (jz 2 | - V&) ctei A dz x A cfe 2 A dz 2 , (17) 
which using cfei A dz\ A A dz 2 ~ |zi|d|zi| A d{argz\) A |z 2 |d|£2| A d(argz 2 ) yields: 

Cn~y y |zi| 36 |zi|<5(|zi| -V^)d|«i| /" N 38 N<5(N - V^) d|z 2 | 

37,39_ 1 10 

~V 36 36 = V 9 ; 

Q22 ~ y y zi 8 ^! 8 ^ (kil ~ ^) 5 ^|2 2 | - ctei A dz\ A d«2 A dz 2 

~^y ki| 36 ki|«y(ki|-v^)d|^i| y |^ 2 | 2 N 2 |5 (k 2 | - vs?) d|z 2 | 



C12 ~ y J z\ % z\ 9 z\ s z 2 5 (\zi\ - Vse^ (5 (jz 2 | - Vse j dzi A ctei A dz 2 A cte 2 
I^I^I^I^I^I-Vm)^! y N 2 °N<5 (|z 2 |-VM)d|z 2 | 



37 i 21_1 11 
~ V 36 " r 36 = V 18 i 



Ci3 ~ y y 4 8z 2 9 ^i 8 ^2 7 ^ (kil - ^ (M - V 3 ^) dzi A dzi A cZ2 2 A cfe 2 

~^y kii 36 kii5(izii-v^)dizii y iz 2 i 56 k 2 i5(iz 2 i-vi6)diz 2 i 

3I_l5I_1 29 

~ V 36 36 = V i g ; 

C33 ~ i y z\ s zfz\ 8 zf5 (|zi| - V 3 ^^) 5 (\z 2 \ - dzi A dzi A cZ2 2 A dz 2 

3^ d\zi\ y |^ 2 | 74 |2; 2 |(5 [\z 2 \ - d\z 2 \ 



y I kir|:i|^(|--j.| - V- 



37 1 75_1 19 



C23 ~ y y z\ 8 z 2 zl 8 z 2 7 6 (\zi\ - V 3 ^) 5 (\z 2 \ - dzi A cbi A <iz 2 A dz 2 

~^y kii 36 ki^(kii-v^)d| 21 i y M 38 N5(M-v^)diz 2 i 

37 1 39_i 10 

~ V 36 36 = V 9 ; 

C44 ~ y y zfzfzfzfd [\zi\ - V 3 ^) 5 (\z2\ - V^j dzi A dz x A cZ2 2 A dz 2 
~iy |zi| ra |zi|j(|zi|-V^)d|zi| y k 2 r 4 k 2 |<5(|z 2 |-V^)d|z 2 | 



73, 75_i 28 



C14 ~ y y ^i 8 ^ 9 ^? 6 ^ 7 ^ (kil - - V 3 ^) dzi A d^i A d2 2 A dz 2 

~y J kil^kil^^il-V^dl^il y \z 2 f% 2 \5(\z 2 \-V^)d\z 2 \ 



£2_i_ o£ _ 1 



~ V 36 36 



19 

.18 „ ^36-37< 



C24 ~ y / z™z 2 zfz 2 5 [\zi\- V^) 5 [\z 2 \ - V 3 ^ ) dz\ A dz\ A cfe 2 A dz 2 



hi 

y y |zi| 54 |zi|5(|zi|-V^)d|zi| y |z 2 | 38 |z 2 |,5(|z 2 |-V 3 i3)d|z 2 | 



-1 



29 

V is : 



~ V 36 36 

C34 ~ — I z\*zfzfzf5 (\zi\ - 5 (\z 2 \ - V® ) (lit '•• >7-.i - d: 2 A r/:_> 



y y I^I^I^I^I^I-Vm)^! y |z 2 r 4 |z 2 |<5(|z 2 |-V^)d|z 2 | 



~ V 36 36 



47 
VT8. 



(18) 
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Lets now look at the term quadratic in the L>3-brane position moduli in T a= B t s which is given by: 

W^^-j^^)- (19) 

where 

(J>~y = WVfcmj ^ van —jiKimi \/va/hml ^Q) 



p||2 9 i2h 9 j3 ~ j3 9 k2k ^ i2j2k2 n hh - k2 



As aj a forms a basis of -fT^C^) and Xa forms a basis of H^^CY^), this does not therefore depend 
on the choice of the divisor a. By Griffith's residue formula, £1 = Q,i 2 ±dzi A dz 2 A dz± where 

UlM ~ dP(zuZ2,z 3 ,Z4) - 3Z% - lf) Zl Z2Z A ' ^ 

OZ3 

From P(z±, Z2, Zs, Z4) = 0, one obtains: 

^Z\Z 2 Z/± 

Z 3 ~ j- 

3 



( z 18 + z 18 + z 2 _ ^6) + ^/^ ZlZ2Z4 y + ( z 18 + z 18 + ^2 _ ^6)< 
+ (" (4 8 + 4 8 + *2 - 4>zfzl) + yJ^Z^f + (z« + Z 18 + z | _ ^6)^ 3 ^ (22) 

which for zi ; 2 ~ Vse , z 4 ~ V 6 yields Z3 ~ Vs . Hence, 

f} 124 ~V~i (23) 

The polynomial deformation coefficient ^ is a complex structure deformation modulus and given that it 
must be of the type z a , i.e., odd under the holomorphic isometric involution a : z\ — > —Z\, 2:2,4 — > ^2,4, 
a : ip — > —if;. So, given that <7*Q = — hence a : O124 — > ^124- 
Now, 

1 ^ sj (Xa)m v ^ 



- {Va) l ' — " - — 
2 ^124 ^ 



~ O(10). (24) 
v~io 6 



Assuming the complex structure moduli z a are stabilized at values: 0(10) z a ~ O(l) for V ~ 10 6 , 
one sees that contribution quadratic in z l z^ goes like {oj^f^z 1 ^). ojb,s are Poincare-duals of T,b,s 
respectively. Hence, ojb,s = <5(-P£ SiS )d-P£ SjS A ^(PT, B ^ s )dPT, BS . 

. Noting that »* Cfe 3 ~ ^4^dz l+ z ld z 2 )-^ dz^dz,) ( . Eb(1+ 18 + 18 + z | = ^6) ^ C y 3)) 

(^«4- 2l 18 -zl 8 -l)3 



which 



1 17 

near 2:^2 ~ Vss implies cfe3 ~ V^{dz\ + dz 2 ). Hence, near zip ~ ^7=-, ~ ^rr(^i + 



dz2) A (efei + ti^) ~ (^1 + dz 2 ) A (dzi + dz2). 

lv~io 6 

Noting that i* s dz 4 ~ Mf|(£^^H^+4!^) fa . £ 5 ( 1 + ^18 + ^2 = ^ ^ ^ 

12 1 17 

which near z\p ~ V» implies c?Z4 ~ Vs(cfei + 0^2). Again, hence 2^2 ~ ^J"' 6 ^ ~ "^-(^1 + 



dz 2 ) A (dzi + ^2) 



~ (dzi + dz 2 ) A (dzi + dz 2 ). 

V~10 6 
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K 
Mf 



-2ln 



i 

We hence see that {uJB)fj ~ i^s)ij ~ C(l) near zi j2 ~ -^=-. 

Therefore, near c a , b a ~ ^ < l(in units of M p = 1), n ~ 0(1), /c a ~ 0(10), 

|2 



Tb + Tb 



r2 {\zi\ 2 + \z 2 \ 2 + z ± z 2 + z 2 zi} ,,io|ai| 2 .,11 (aia 2 + /i.e.) ,,I Kl 2 ,,22 (oia 3 + /i.c.) 



+V1T 



(a 2 a 3 + /i.e. 



+ V^ 



12 a 3 



+ V s (0104 + 0401) + Vis (a 2 a4 + a4a 2 ) +V« (0304 + 0403) + V 9 |a4| 
3/2 



3/2 



+ ;2 {>i| 2 + N 2 + ^ 2 + ^i} > , V- 0/ 



Mr, 



M P 

One can argue that near 

J_ 1 _2 _1 _13 _11 4 

Zl,2\ ~ VssMp, |^3 1 ~ VeMp, |ai| ~ V 9M p ,|a 2 |~V sM p ,|a3|~V isM p ,|a4|~V 9 M p , C = 0; 

M p I , one obtains a local meta-stable dS-like minimum corresponding to the posi- 



(25) 



G a ~ o(i)k a {~o{ia)) ±VJ -p ) 
tive minimum of the potential 

e K G Tsfs \D Ts W\ 2 stabilizing vol(S B ) ~ Re(a B ) ~ vi,vol(S 5 ) ~ Re<r s ~ Vts such that Re(T) 5 ~ 
Vis and in the dilute flux approximation: 



Re(T B ) + lnlP(J: s )\ D3ls \+lnlP(X s 



lD3|j 



9j=su(3) or 5(7(2) 

+C (u(l) - Flux 2 ) ~ Vn. 
Indeed, near the aforementioned stabilized values of the open string moduli, 

I Ci 1 1 a 1 1 2 1 ~ V "8" , \C X 2 (aio 2 + h.c.)\ ~ V^, |C 1 g (aia 3 + /i.c.)| ~ V§, 

ICxj (aia 4 + /i.c.)| ~ V» , |C 22 |a 2 | 2 | ~ V~i, |C 23 - (a 2 «3 + h.c.)\ ~ Vis 

— 1 9 1 — — 

|C 2 2 (a 2 a 4 + /i.c.)| ~ Vis, I C*33 1 a 3 1 I ~ V 3 , |C 34 - (a 3 a 4 + h.c.)\ ~ V» , 

1 9 1 2 

IQjM I ~ Vs, 

we see that there is the possibility that: 

Vol(E B ) + C l3 aiaj + h.c. ~ V^, 



(26) 



(27) 



(28) 



and hence (1261) could be implemented. 

To obtain the "physical" /normalized Yukawa couplings, soft SUSY breaking parameters, vertices 
later in the paper, etc., one needs to diagonalize the matrix generated from the mixed double derivative 
of the Kahler potential, the same given as under: 



1 

y2/3 
1 

y2/3 



^373 a 13y57l2 a 14 v il/l2 ffl 15 WV ai 6 V 7 / 12 ^ 

a 23y57T2 a 24 v n/i2 a 25 



• V 



1 



034 IS 



^35 



V 17 / 18 a3 fi V 13 / 9 



/ an 

an y 

a 13y57T2 a 23y57T2 °33 

°14 v ii 1 /i2 a 24 v ii/i2 034^ 

ojbW «3 5 V 17 / 18 a 45 V 4 / 9 a 55 V 13 / 9 a 56 V 35 / 18 

V«ieV 7 / 12 a 26 V 7 / 12 a 36 V 13 / 9 a 46 V 17 / 18 a 56 V 35 / 18 a 66 V 22 / 9 / 



,s- a 44 ^ a 45 V 4 / 9 a 46 V 17 / 18 



(29) 
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With V ~ 10 6 and some fine tuning in an, 12, 22, 14,13,23,24,36,44,46,56,66 > the numerical eigenvalues are 
estimated to be: 

10 12 , 10 7 , 10 4 , 1(T 2 , 1(T 3 , 1(T 5 , (30) 
with the corresponding eigenvectors given by: 

_ 3 _ 6 _ 9 _ 2 

^4~o 4 + V sa 3 + V sai + V ^>a 2 + V [z\ + z 2 ) ; 
A3 ~ —03 + V~^a 4 — V^tai — V - sa 2 + V~t (zi + Z2) ; 

_3 _i _6 _6 

Ai~ai-V sa 3 + V a 2 - V sa 4 + V ?{zi+z 2 ); 
A 2 ~ — 02 - V _1 ai + V _ 5a 3 - V - ^ (z\ + z 2 ) ; 

(Zl + Zo) 6 3 8 o 

^2 ~ ~ r ~ V~sai + V-5 fl2 + V^sa 3 + V" 2 a 4 ; 
v2 

(Zl — Zo) 6 3 8 , 

Zi ~ v 21 - V*ai + V"a 2 + V"5a 3 + V~ 2 a 4 . (31) 
v2 



The system of equations ([31]) can be solved to yield: 



Zl „ (^3) _ v -l Al _ V -I A2 + v -f ^ 3 + v -2^ 4 . 
v2 

ai ~ V^t^i + V~t,Z 2 - Ai - ^ + V _ §-4 3 + V~f A±] 

3 Z> j\. 7 9 

a 2 ~ V~i Zi + y - -j- + -4 2 - V~^ 3 + V~f ^4; 

a 3 ~ V~^i + V _ tz 2 - V - f .Ai - V~^A 2 - A 3 + V"f ^ 4 ; 

^ 11 »/4. 9 3 

a 4 ~ -i + V~^£ 2 + -f - V"^U 2 + V~Ia 3 + A 4 ; 



(32) 



The non-perturbative superpotential we will be using is given by: 

W ~ (l + 4 8 + 4 8 + (3(^0^4 " zl 8 - zf Y - e -n°vol(Z s )-v 3 (a sZ l +l 3 sZ l + ~ /s z lZ2 )_ (33) 

For n s = 2, vol(Yls) and were stabilized at around Vis and Vse. Now, considering fluctuations in 
Zi : Zi — > Vse + 5zi, with the fluctuations expressible in terms of ^1,2,^1,2,3,4 using ([32]) . and with 
the understanding that Vis (1 + as + fts + Is) ~ ^ n ^, one obtains the following non-zero physical 

effective Yukuwa couplings :Yq c c = , 1 3 k , Q being an open string modulus which for 

us is 2, 5^.1,2,3,4, where, e.g., y^A 1 2^3 4 * s ^(-^i ^coefficient in the following mass term in the 
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M = 1 SUGRA action of [12]: efv A D A Wx Al/3 X A ' 2/4 ■ Using ([32]), and: 

efv zi D Zl W = V~%, e%V Zl D ai W = V"W, ef V Zl D a , W = V"^, 

25 17 A" 71 

eTrD jei D aa W = V-w, eT^^W^ = V" w , e~D ai D ai W = , 

efp a2J D a2 T^ = V"£, efp a2J D a3 ^ = V"f, ef P a2J D a4 lF = V"§ , 

ef P a3 L> a3 lK = V^, efp a3 D a4 Ty = yl, ef V a4 D a4 W = (34) 

one can verify that e%V Al D Aj W ~ ef P a/ D aj Ty. Now, (ef X^gPF) x£xg = (eTV t DjW} X X l Xr, 
where A/E index ^4/,Zj and X/J index a/,Zj. What is interesting is that using (f3Tj) . ([32]) . e.g., 

e * V Ah2 D A34 W) Xl X R ~[e 2 V- ai2 D- a , iA W } x L X R = J J [e * V- ai2 D- a ^W j % L X X R 

~ (ef V- ah2 D- a3A w) xT 2 xT ■ (35) 



This implies that the mass terms, in the large volume limit, are invariant under diagonalization of the 
open string moduli. 



G(Zi) term in e~V M D M W _ ~ eff 



\J K Z l Z i K AxAx K MM 



1X4^3- 10 ~ d xV-5, (36) 



which implies that for V ~ 10 5 , (Zi) ~ 246GeV, {Zi)Yz lAlAa ~ MeV - about the mass of the 
electron! 

which implies that for V ~ 10 5 , (.Zj) ~ 246GeV, (2j)iz i> A l ,A2 ~ WMeV - close to the mass of 
the up quark! 

There is an implicit assumption that the vev of the L>3-brane position moduli, identified with the 
neutral components of two Higgs doublets, can RG-flow down to 2A6GeV - that the same is possible 
was shown in [10]. The RG-flow of the effective physical Yukawas are expected to change by 0(1) 
under an RG flow from the string scale down to the EW scale. This can be motivated by looking at 
RG-flows of the physical Yukawas Iasa in MSSM-like models: 

= Ta^ea + 1%Y AkA + 7a^as k , (38) 
where the anomalous dimension matrix 7^, at one loop, is defined as: 



7 ^ = ^2 ( W^ T -2^^ a)($A) ^j ; (39) 

(a) 
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where (a) indexes the three gauge groups and A, etc. index the diagonalized basis fields Ai, Zi. Using 
(|33p . one can show that at M s (string scale): 



_1 a _79 a _31_ - 143 

YziZjZ,, ~ V s , Y ZiZiAl ~ V « , Y ZiZiAa ~ V « , y M/4 ~ V 40 , 

107 a _163 a 23 

Xkj^a ~ V 40 , ~ V 40 , Y ZiAlM ~ V s , 

191 *> 227 ■*> 67 

143 a _179 - 219 

_51_ - _55 ~ 247 

YziA 3 A4 ~ V 8 > Vzi^UAi ~ ^ 8 . ^i^i^i ~ V 40 > 

199 a _55 „ 311 

W ~ V 40 » X4iA^3 ~ V 8 ) ~ V 40 , 

151 a _ 227 a 263 

Y AlA2A2 V 40 , Y AlMM ~ V 40 , ~ V 40 , 

339 a _75 „ 103 

^42^4^5 ~ V 40 , Y^AiJ^ ~ V 8 , y^ 2 ^ 2 ^ 2 ~ V 40 , 

179 a _43 „ 299 

^-4*43-43 ~ V 40 . X4 2 ^2^4 ~ V 8 , *A*4a*» ~ V <» , 

263 a _ 367 a 283 

Y A2AaAi ~ V 40 , y4 4 ^ 3 ^ 3 ~ V 40 , y^,^^ ~ V 40 , 

327 a _ 403 a _331 a _439 

Ya 2 MM ~ V 40 1 ^.AA ~ V 40 , ^43^3^3 ~ V 40 , ~ V 40 . (40) 

From (|40p . one sees that Yz^z^Mg) ~ V~s is the most dominant physical Yukawa at the string 
scale. Hence, 

F a 1a . 2E( )4,) cr («)(*r)^rEA , . 

7a^tsa ~ ^2 YkZiZiYz.z.zYz&A ^ • ( 41 ) 

Now, the first term in (|4ip is volume suppressed as compared to the second term at M s ; let us assume 

a 2 (t) F 21 

that the same will be true up to the EW scale. Using the one-loop solution for j"^ = 1+ ^ a) t in ([4"Ijl 
and therefore (1551). one obtains: 



whose solution yields: 

3 -2(c (a) (A)+C (a) (S)+C (a) (A)) 

Y ASA (t) ~ Y A2A (M S ) [J (l + /3 (o) t) 6 W . (43) 

(a)=l 

The solution (j43j) justifies the assumption that all Yaeas change only by 0(1) as one RG-flows down 
from the string to the EW scale. 

A similar argument for the RG-evolution of Y^ AlAj 8 would proceed as follows. We will for 
definiteness and due to relevance to the preceding discussion on lepton and quark masses, consider 

dYS. B . . dY cfi 



and _^MM_ NoW; 



ul Z lAl A 3 A -Oeff I _ A vcfi i A i>eff ( AA \ 
J t = I/Z1YAA1A3 + 7Ai Y ZiAA3 + 4 Y 2iAA' ( 44 ) 
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where: 



yZ x Z^Yiz x Z V etc - 



(45) 



implying: 



aY z 1 AiA 3 
dt 



9 Yz 1 z 1 z 1 Yz iZi z 1 Yz 1 AiA 3 +Yz 1 z 1 z 1 Y£ 1 z 1 z 1 Y3£ AlA3 + Yz 1 z 1 z 1 YX 2 z 1 z 1 y X^AiA 3 



cff 



^eff 



Yz 1 z 1 z 1 YX i z 1 z 1 Y XAiA 3 + Yz 1 z 1 z 1 YX i z 1 z 1 YA A AiA 3 ^ 

dY cS 

and a similar equation for — Zl ^ 2Ai ■ Using (|40p and: 

ef V ai/2 D as/4 W 3 V~l (zi - Vs) + V§ fai - V—tt ) + V~~ { a 2 - V~~ 

i / 13 \ r / ii 

+Vs ( a 3 - V" is J + Vs ( a 4 - V T 

and (1151) , one sees that at M s : 



(46) 



(47) 



..A v-eff | „ A \>eff , .A \>eff = 
TZi r A^!^ 3 + 7^! r ^iA^ 3 + 7,4 3 r 2i^iA 3 



_,A v-eff i _,A v" e ff i „A i> c ff =s 
7-Zj r A^ 2 ^4 + 7^2 r -ZiA*4 4 + ^A i Y Z 1 A 2 A = 



_ 1_1_4 
V 8 8 9 



V K Z 1 Z 1 K AiAi K A 3 A 3 



233 

~ V 180 , 



V~10 5 



1/889 



\/ K Z 1 Zx K AiA2 K AaAa 



43 
V~36. 



(48) 



v~io 5 



Hence, at the string scale, the anomalous dimension matrix contribution is sub-dominant as compared 
to the gauge coupling-dependent contribution, and by a similar assumption (justified by the solution 
below) / reasoning: 



3 -2(c (a) (A)+C (a) (S) + C (a) (A)) 

*AEA(*) ~ YH^{M S ) J] (1 + {a) t) "(a) 
(a)=l 



(49) 



This suggests that possibly, the fermionic superpartners of A\ and ^3 correspond respectively to ej, 
and and the fermionic superpartners of A2 and .A4 correspond respectively to the first generation 
ul and ur. 

We will consider four stacks of Z)7-branes: a stack of 3, a stack of 2 and two stacks of 1. The 
matter fields: L-quarks and their superpartners will be valued in the bifundamentals (3, 2) under 
SU(3) C x SU(2)l; the L-leptons and their superpartners will be valued in the bifundamentals (2, — 1) 
of SU(2)l x U(l)y- Before the fluxes (|50l) are turned on, the Wilson line moduli are valued in the 
adjoint of U(7). With the following choice of fluxes: 



/ fx 

























fx 

























fx 

























h 

























h 

























h 







V 

















h 


) 



(50) 
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the U(7) is broken down to J7(3) xU(2)x U(l) xU(l). So, the bifundamental Wilson line super-moduli 
Ai will be represented as: 



A! 



[(a, &) = !,..., 7], 



(51) 



where (e a6 U = <5 oi %. Hence, 



A = 



-A 2 



/ 























°\ 











































































+ t 
























e + t 


)e 











P e 


+ 6v e 


g + 8e 












V o 























o / 




/ 











u + 6u 








o \ 




























































u 


+ M 














































































{ 




















o ) 





and 



^3 



Aa 



( 

















\ 








































































































en + 9e R 


V o 














e R + ee> R 


o / 


/ 
















ur + 0u R \ 


















































































^ + 















\ u R + 6u l R 







I 



(52) 



(53) 



(54) 



(55) 



Unlike the Ai which correspond to matter fields corresponding to open strings stretched between two 
stacks of D7 branes (with different two-form fluxes turned on their world volumes), the Higgses would 
arise as a geometric moduli corresponding to the fluctuations in the position of the mobile space-time 
filling D3-brane. 

Now, we will see if one can construct appropriate a/ (bi- fundamental) and Z{ (mimicking bi- 
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fundamental fields) such that (j3T|) . ([52 



and 



m, 




- 


m, 


as well 


as: 








1 





















\ 







































































+ 9H U 









































+ 0&1 













v° 





















0/ 


( o 





















\ 
































































































+ 9Hd 

















H d + m\ 










v° 
























(56) 



(57) 



are satisfied. Guided by (|5ip . we will make the following ansatze for a/ and Zi~. 



ai 



V 












af 


a} 5 


a} 6 











af 


af 


af 











af 


af 


af 


af 


-af 


af 








af 


af 


af 


af 








af 


af 


af 


af 


af 


af 





af 


af 


af 


af 


af 


af 











4 4 


-15 


-16 











*f 


-25 


-26 











-34 

■S 


-35 

*i 


-36 

■*i 


z^ 


Z 2A 


? 34 








-46 

H 


-15 
H 


? 25 










-56 
<*« 


5I6 

H 




536 


? 46 


? 56 





z 17 


z 27 


? 37 


-47 


557 


567 



A 7 

-.57 

,67 



7 

/ 

27 



2 

37 
i 

47 

i 

57 

j 

67 



(58) 



/ 



Using (|58p in (|3ip . one obtains 17 x 6 = 102 equations in 102 variables. The same can be solved (using 
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Mathematica) to yield the following solution: 



ai = 



a 2 = 



a-3 = 



tl4->,--I " 



-1-7 11 - 



( 



( 



£ 4 V-s»l 







0.3£ 2 14 S L 







S3 v 















£ 













r d6 











0.3£ 2 14 n L 
















4l 6 











e f + v~*h. 





( 







i v 11 

Zl 7 v-*u R 









\ 



si 6 (- 



13 - 
5 UR 



















o o el 6 









~h L + V f>H, 



e!4u L 
S3 V 









'.) 





















Hd 



{-V-le L + V"f 5^V-^H 







btfV-r>H d 


-^v-Ur 





; (-V-te L +V-fiJ u ) 













Si? 



\ 



67 V-fe R 



M6 
S3 





o 
o 
o 

0.3£fe fl 



04 



0.08^l 4 u L 





/ o 





0.08^ 4 u L 



o o o tfv-ii L efv-Ti^ 

V 0.2£l 7 u R 





























J 



(59) 



17 



and 



21 



a 







14 

1 T 










Z2 



4 _ 

5 UL 































af 



u i v 







V 5e L 



V2 











af^ 



a 



16 







9 







V 5 ul 









a 



v- 



a— tt 

K p T ±±M 

5eL 72 











_056#d 

"2 



56 

1 73 




Si. 1 

V2j 



17 — — 

5a{ V 5 u_r 







afV~ 



afV" 









af 



6V- 



a 2 



e-L - 

Hd 
V~2 



J 

17-. <_M 
5«2 V 5 










\ 5al 7 V~u R 





67 efl 
2 



a 



67 e fl 
2 v=7 



(60) 



In (|59p and (|60|) . l<a<6, 1<6<7 are 0(1) numbers. Now, from (|18|). the intersection matrix 
can be written as: 



C 



u 



( cnV 10 / 9 ciaV 11 /^ ci 3 V 29 / 18 C14V 19 / 9 \ 

Ci 2 V n /18 C22 ^7 C23V 10/9 C24V 29/18 

C13V 29 / 18 c 23 V 10 / 9 C33V 10 / 9 C34V 47 / 18 

V C14V 19 / 9 c 24 V 29 / 18 C34V 47 / 18 C44V 28 / 9 / 



(61) 



a=i,...,fti- 1 (cr 3 ) 



where c a t, 1 < a, 6 < 4 are 0(1) numbers. Assuming that the complex structure moduli z 
are stabilized at very small values, which is in fact already assumed in writing (J5|) which has been 
written upon inclusion of terms up to linear in the complex structure moduli, let us define a modified 
intersection matrix in the aj — Z{ moduli space: 



1 = 1, J = J: 



C XJ = C 
C IJ = fi 3 {2-na'f (u a f , X = i,J 
C XJ = 0, l = I,J = j, etc.. 



Now, from (1591) and (1601). one can show that for 



0.8V5S, V ~ 10 5 : 



(62) 



C XJ Tr{MiMt) ~ C ai ~ ai \e L \ 2 + C 



" 2a2 \u L \ 2 



C a ^\e R \ 2 + C a ^\u R \ 2 + M3 (2W) a |i/n| 2 ,(63) 



where A^x = aj,Zi. Now, using (|3ip . one can show that in the large volume limit: C-^-^j ~ 
C aia J ,C ZiZ j ~ C ZiZj . Clubbing together the Wilson line moduli and the _D3-brane position mod- 



(64) 



uli into a single vector: M.\ = Ai,Zi, then one sees from (|52j) - (|57|) . and (j63|) : 

C A2 Tr (.Ma-MQ ~ C XJ Tr(7W x 7W^). 
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In the large volume and rigid limit of T,b{C, A = which corresponds to a local minimum), perhaps 
C A ^Tr ( MhML J is invariant under moduli transformations in the (Aj, Z{)/(aj, Zj)-subspace of the 



open-string moduli space, which would imply that C IJ aiaj + 1^3 (a 1 ) (ojb)^z 1 z^ for multiple D7- 
branes, in a basis that diagonalizes g^xM - a ^ stabilized values of the open string moduli, is replaced 



by ([631). 

For the purpose of evaluation of (N)LSP decays in the subsequent sections, we will be using the 
following terms (written out in four-component notation or their two-component analogs and utiliz- 
ing/generalizing results of [IB]) in the M = 1 gauged supergravity action of Wess and Bagger [12] with 



the understanding that m moduli/ / modulini << m K K I ~ 



V~10 5 / 6 



\\) n G<\~ I . A/, = ^ 



V~10 5 / 6 



10 15 GeU, and that for multiple £>7-branes, the non-abelian gauged isometry groud_|, corresponding 
to the killing vector 62^4^7 (2tto.') Qb9t b , Qb = (2na') J s i*ws A P-f arising due to the elimination 

(2) 

of of the two- form axions D B in favor of the zero- form axions p B under the KK-reduction of the 
ten-dimensional four-form axion 116] (which results in a modification of the covariant derivative of T B 
by an additive shift given by §in\p7 (27ra') Tr(QBA^)) can be identified with the SM group (i.e. 
is the SM-like adjoint- valued gauge field |12j): 



C 



9YM9T B j Tr \ xTB xi^g, r) +i9xj Tr 



xl 



0x1 + n^a M x£ + - {d aM K$a M - c.c.) xl 



K_ 

+ ^ pxDjW) Tr {xlxi) + g TB T B Tr [{8„T B - A^) (d»T B - A»X T ^ 

+g TB jTr (X^A^xiYl^, r) + V'L, ^ X -f^, L F pX + $ L , ^7%, L W+W^ 



+Tr \ St L 4 ( 6nl^(^a')Q B K + 

KqT b T b 



\2n\p 7 {2-Ka')Q bv 
V 



B 



K 9, L 



+ - 



6iKi(2Tra')Tr QsA^d^ ( k^^ttc/ ) 2 C IJ aiaj 



+ h.c, 



(65) 



3 NLSP Decays 

The generation of mass scales of superpartners and realization of p split SUSY in the context of Big 
Divisor D3 — D7 type IIB string compactifications provides us with the gravitino as the L(ightest) 
S(upersymmetric) P(article) and sleptons/squarks as N(ext-to) L(ightest) S(upersymmetric) P(article)s 
with (Bino/Wino-type)gaugino/(Bino/Wino-type)gaugino-dominant neutralino, in the dilute-flux ap- 
proximation, only differing in their masses by an 0(1) factor. This helps in shedding some light on 
identifying a viable dark matter candidate in this framework. Decays of the neutralino, are, in gen- 
eral, driven by any of the trilinear R-violating couplings which necessarily involve squark or slepton 
as propagators. In MSSM, stability of LSP is governed by conserved R-parity however same is not a 
restrictive condition in " p split SUSY" . Some of the R-parity- violating interactions are not necessarily 

3 As explained in [IS], one of the two Pecci-Quinn/shift symmetries along the RR two-form axions c a and the zero-form 
axion p B gets gauged due to the dualization of the Green-Schwarz term / R i. 3 dD^ A A coming from the KK reduction 
of the Chern-Simons term on Es Uct(E_b) - being an RR two-form axion. In the presence of fluxes ([50} for multiple 
D7-brane fluxes, the aforementioned Green-Schwarz is expected to be modified to Tr (q b J r i,3 dDjp A A*j, which after 
dualization in turn modifies the covariant derivative of T B and hence the killing isometry. 
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negligible and the consideration of LSP to be a viable dark matter candidate needs the contribution of 
the same to be evaluated. Here, the large squark masses in /x-split SUSY helps to suppress the decay 
width. In BSM/2HDM models, one considers R-parity violating coupling generated from R-parity 
violating superpotential. However, here we consider the possible R-parity violating couplings in the 
effective ftf = 1 gauged Supergravity action and by evaluating the same, we explicitly calculate the 
decay widths and hence life times for neutralino/gaugino and slepton decays. The neutralino mass 
matrix, for the four- Wilson-line- moduli setup, similar to [TT], can be diagonalized and the smallest 

2 

eigenvalue V37713/2 corresponds to the neutralino: 

xl A + / (til + H^j + charge conjugate, (66) 

where Hf 2 are the Higgsinos. We will first evaluate the lifetimes corresponding to two- and three-body 
decays of Wino/Bino-type gauginos and the gluino. We will then evaluate the lifetime for three-body 
decays of the neutralino followed by the lifetime corresponding to two- and three-body decays of the 
slepton. 

The Big bang nucleosynthesis predicts the universal abundances of D, 3He, 4He, and 7Li, essentially 
fixed by average lifetime t ~ 180sec. In our study, we are considering both radiative as well as 
hadronic decay modes of N(LSP)'s. Therefore, if life time of NLSP's are more than 10 2 sec, the 
high energy photons emitted via radiative decay might destroy the abundance of light elements by 
inducing photo-dissociation of same, hence reducing the standard baryon-to-photon ratio and resulting 
in baryon-poor universe . Similarly, during the hadronic decay, the released hadronic energy can 
produce mesons/charged pions causing interconversions between background proton and neutron (p 
/ n conversion) that alters the neutron-to proton ratio, resulting in the change of 4He abundance 
and hence destroying this success of standard BBN. Therefore, if one assumes that all of the dark 
matter density are produced by non thermal decays of NLSP decaying into LSP, it is first important 
to determine whether life time of these decays do not effect standard abundance of D, 3He, 4He, and 
7Li, and hence predictions of BBN. We will see that one-loop gluino decay into gravitino violates the 
BBN requirements, and hence the gluino can not be an NLSP. 

3.1 Gaugino Decays 

In this sub-section, we discuss two- and three-body decay modes of the Wino/Bino-type gauginos as 
well as the gluino. 



3.1.1 Two- and Three-Body Gaugino(Wino/Bino) Decays 



The relevant gravtino-gaugino-photon/^-boson vertex in the M = 1 SUGRA action of |12| is obtained 

from^V^VA^ 
is given by (See |24J): 



from j]-il) p (T flp 7 P A^ F^y , (a) corresponding to the three gauge groups. The decay width for B — > 



m l/2 M P 



4to ~ V"3M p . (67) 



So, (|6"7|) implies a lifetime of around 10" 



-30, 
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Figure 2: Three-body gaugino-decay diagram 



Similarly, the decay width for B — > ip^ + Z is given by (See |24j ) : 
T ( B —?■ + Z 



cos 2 9 w rn° § 



48ttM| m\ j2 



m 



1 



m" 



3/2 



1 



3/2 



B 



1 2 t" + 










m| ( 


f 3+ ^/2 






m| 


I "4 







+ 



m 



B 



m z 
ml 



12 



x z 
4 



m 3/2 m| 



(68) 



"3/2 P 

So, (|68p implies a lifetime of around 10 _30 s. 



Given that the gauge boson-quark-anti-quark vertex will be accompanied by — v N . 2 - 



K 



MM, 



InV, which for V ~ 10 is approximately /V 3 , using further the results of [24] one obtains 



the following result for the decay width A — > ipu + u + 



u: 



T ^B -)• ip^ + u + u 
which yields a lifetime of around 10~ 13 s 



9ymP v 3m |o 
32 (2ir) 3 m\ /2 M% 



/ 2 V" 4 M p , 



(69) 
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(a) (b) 
Figure 3: Three-body gaugino decays into the gravitino 



We will now discuss three-body gaugino decays into gravitino mediated by squarks. To avoid 
channel overlap, for simplicity, we will assume that the gaugino decays mediated by vector bosons and 

those mediated by squarks, involve different gauginos. Utilizing the conditions = 7 M Wt = 0> 

the amplitudes for the above two diagrams can be written as: 



1 



M {a) ~ 2p^ ( 4 + W«(P?) I ! 



p~ — m~ + ze 



u(p q )Vx- q -q\l 



(+) 



M (b) ~ 2 Pl Jf ( u (Pq)^i } (PV>MP. 



1) 1 I 2 2 1 

p~ - m| + ze 



(70) 



V\-q-q ~ /V 5 being the gaugino-quark-squark vertex. From ([7U|). one obtains the following helicities 
averaged sum: 



E 



E 



Tr 



m, 



M b) M\ b) ~ \V M rj^Tr V${t,+ 



m. 



Tr 



+ m x - g [-p. + m, 



(71) 



and utilizing u{pq)ijjfi ijPip) = — V> \Ptl>) v (Pq)-> v (Pq) = — u (P<j)Ag to rewrite My,) and therefore 
obtain: 



In (f70|) - ([72]) . the positive and negative energy gravitino projectors are given by: 



(72) 



*3/2 ± m 3/2 



_ P3/2,^3/2,u \_ 1 / _ P3/2,/J>3/2, g W _ P3/2,XP3/2,X \ a X 



"3/2 



"3/2 



"3/2 



(73) 
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Figure 4: Contact-vertex three-body decay diagram 



A typical term that one would need to calculate is: 




, 7] = ± corresponding to p v = Pqjpq 

{VP3/2 ■ P V ± ™3/2™»?) 



Pq -P3/2P 3 / 2 \ ft] 



"3/2 / \ m 3/2 j 

Using results of [25], to get an estimate of the decay width, one sees that: 



tPv -Pa/2 ±fn z / 2 m n 



(74) 



r w 



Max 



<? + <? + 

dS23 



-S23 



S13 = 



ds 



13 



UM\V\-q-q\ 2 X ( PXg ■Pq+m X .m q ) 



Max(l0~ 32 ,10~ 30 ,10 



! 23 

-32 



( s 23 



m 2 ~) 2 ,(s 2 



)M p 



10~ 28 M 



implying that the corresponding lifetime would be 10' 



13 



-15, 



13 



ml) 



(75) 



In the approximation, 



i»ir 



"13/2 



0, the decay width for W — >■ tp^ + W + + W is given by 



(See [24J): 



r (w° -> % + w + + w~ 



2 9 
9YM m w 



34560 {Ivy m^M* 



(76) 



which implies a lifetime of around 10 s. Treating the Wino/Bino gauginos as the co-NLSPs, we 
hence conclude that the two- and three-body gaugino decays into the gravitino (LSP), respect the 
BBN bounds. If one were to calculate similar neutralino decays into the gravitino, then for the 
Higgsino contribution to the nuetralino-gauge-gravitino vertex, one will have the following suppression 



factor: / ( fV •'• ) x 0{zi) term in g\M 



9T R 



x w hich is / 2 V~i"f x lot. For V ~ 10 6 , this 



is approximately, 10 . The corresponding gaugino-gauge-gravitino vertex will, without worrying 



about the gamma matrices, be of O ( jj- ) . Now, even at low Z-momenta, this is around 10 for 



16 



Mz 10 2 GeV. Hence, the neutralino and gaugino decay (e.g., into gravitino and Z boson) widths, will 
approximately be the same. 



23 




Figure 5: Three-body gluino decay diagrams 
3.1.2 Gluino Decays Revisited 

In this sub-section, we will readdress the gluino decays in the light of previous discussion done in 
The earlier discussion was made with an ambiguity in identifying the quarks with right set of wilson 
line moduli. In the four Wilson line moduli framework where the first generation left handed quarks 
can be identified with Wilson line modulus a 2 (strictly speaking A2) and the first generation right 
handed quarks get identified with 04 (strictly speaking A4) as mentioned above, by recalculating the 
moduli space metric, we will briefly estimate Gluino decay life time results for tree-level as well as 
one-loop decays into neutralino and Goldstino. 

i qqxn 

We first discuss gluino three-body decays that involve the process like g — > qqxn - 9 being a 
gaugino, q/q being quark/anti-quark and Xn being a neutralino. More specifically, e.g., the gluino 
decays into an anti-quark and an off-shell squark and the off-shell squark decays into a quark and 
neutralino. The Feynman diagrams involve gluino-quark/antivertex- squark vertex and neutralino- 
quark/antiquark-squark vertex. The gluino-quark/antiquark- squark in gauged supergravity action as 
also given in (|148p . is proportional to: 

G q ~^ A ' 2 ~ fV~2qx q ^Xg,Xg corresponds to gluino (77) 
and from (|149p . contribution of physical neutralino-quark/antiquark-squark is proportional to 

X f ~ f~V-iqX q/q xi (78) 

Using RG analysis of coefficients of the effective dimension-six gluino decay operators as given in [26], 
it was shown in [11] that these coefficients at the EW scale are of the same order as that at the squark 
mass scale. 

The analytical formulae to calculate decay width for three-body tree-level gluino decay channel as 
given in |25j is:- 

2 f 1 

r(<? -> XnliQj) = 256^3 Af 3 2 / dsi3ds23 2 Ue { Ai ^ S2 ^ + Bi i( s ^) ~ 2 e n e g (s 2 3, S13)) (79) 

9 i,j 
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where the integrand is the square of the spin-averaged total amplitude and i,j = 1,2,..,6 are the indices 
of the squarks mediating the decay. The limits of integration in (|79p are given in . The Aij terms 
in (|79p represent the contributions from the gluino decay channel involving gluino— )>squark+quark and 
squark— )-neutralino+anti-quark, whereas the Bij terms come from channel gluino — )>squark-|- anti-quark 
and squark— meutralino+quark. The same are defined in [IT]. Utilizing the results as given in (|77|) 



and dZSJ), Aij \ Tr 
C\ Tr f 













( Tr 




Tr 







/4 V " 



ir~ B l3 \Tr 



Tr 



G L G V X I X V 



/ 4 V " 5 . Substituting these values for various vertex elements and 



solving equation (|79|) by using the same limits of integration as given in [TTJ , dominating contribution 
of decay width is: 



2567r 3 m 3 

2 
5: 



?4-> 5 -22 4, ,12 , 74. ,_ 6i 10 4 i4,._M 10 4 



256vr 3 V 2 mf 

2 



(80) 



~ O(10^ 2 )/ 4 Ge1/ 

Considering lower bound on f 2 ~ 10~ 8 as calculated in in dilute flux approximation , decay width 
of gluino i.e equation ([80]) becomes: T(g -> XnQiQj) ~ O(10~ 2 )/ 4 < e>(10~ 18 )GeF. Further, Life time 
of gluino is given as 

10" 34 J sec 10~ 23 

(81) 



n 



r 10- 2 / 4 GeF f 4 



> W~°sec 



g -»■ X3 + S 



Relevant to Figs. 4(a) and 4(b), the gluino-quark-squark vertex and the neutralino-quark-squark 
vertex will be given by ([77]) - ([78]) . As explained in [11] , the quark-quark-gluon vertex relevant to figure 
4, from [12] is given by: 



9l jX J a ■ A Im (X B K + iD B ) X 7 , 



ffyM^a^X" 42 ^ • A I 6Ki^ 7 27ra'Q B K + 



V 



(82) 



X^ 2 and x^ 2 to quarks and antiquarks. X = X b 8b = —Vlma! k^h-jQb9t b corresponding to the 
killing isometry vector and D term generated is given by: 

B _ itra 1 'k\^Qbv b 



D 



which utilizing the fact (see equatior il24p : 



V 



9MM ~ V 9 ' 



as well as w ~ Va , ~ V z (2ira') f, yields for the quark-quark-gluon vertex: 



A' 



.A2.A2 



0(10- 



O(10) 2 V^s^ 2 ^ 2 ) /> • e ~ V~^fa • e; for V 



10° 



(83) 
(84) 

(85) 
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Figure 6: Diagrams contributing to one-loop two-body gluino decay 



The gauge kinetic term for squark-squark-gluon vertex, relevant to Fig. 4(b) will be given by 



G aBCB \7 TbSJ Tg. This implies that the following term generates the required squark-squark- 



gluon vertex: 
Uk\w2-kcIQ b G Cb ° b 



Vae- [2k- {p^o +p- g 



K 



A2A2 



O(10 2 )fV-9 2e-k-e-{p^+p~ g 



~ /V«,for V ~ 10 E 



(86) 



In it has been discusssed that behavior of Wilson coefficients corresponding to two-body 
gluino decay does not change much upon RG evolution to EW scale. Using the vertices calculated 
above relevant to Figs. 4(a) and 4(b), and the Feynman rules of [27], one obtains for the scattering 
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Figure 7: Different vertices relevant to one-loop gluino decay 
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amplitude: 



M ~ f 3 M p 



x V* 



d 4 k 



X V~2 



ia ■ k 



k 2 — m 2 + ie 



V" 



(k - Pg) 2 -m 2 ~ + ie 
( 



s , , i d k , , 3 



\ 



V 



V 



ia ■ k 



k 2 — m 2 + ie 



x V so a ■ e 



IG 



Using the 1-loop integrals of [28J: 

1 f d^k (^/i) kfjjtii) 



(87) 



(2tt) (k 2 — m 2 + ie) (k + pi) 2 — m 2 , + ie {k + p\ + p 2 ) 2 — m 2 + ie 



4vr" 



Pl^Cu + p 2fl C 12 ,p lfl p lly C 21 +P2nP2vC 22 + {pinP2u + P\vPl») C 2 3 + V^ C 24 



(a) m 1 = m q ,m 2 = m 3 = m f ; p 1 = ~P^o,p 2 = -p g /, 

(b) m 1 = m 3 = m q ,m 2 = mf, p\ = p^o , p 2 = -p~ g . (88) 

The one loop three point functions C«'s relevant for cases (a) and (b) have been calculated in [11] 
and are given as under: 

C$ = O(10 6 ),c£ ) = O(10 6 ); 

C-W = ^ x 10- 21 ~ 0(1) x W- 23 GeV- 2 ; 

C§ } = (9(10) x 10- 16 Ge^ 2 ; 

C§ = O(10) x 10- 16 Gey- 2 ; 

c® = of = of ~ 0(1) x io- 10 o e y- 2 ; 
c^ 6) = ^ x 10- 21 ~ 0(1) x io- 23 G e y- 2 . 



Now, equation ()87p can be evaluated to yield: 



3-1,-2.1 



' Pxg^H + °" ' p 9 C V2 } CT • efJ • + CT " Pxg a ' ea ' p xi C 2l + & ■ Pg° ■ ■ PgC$ 



a ■ p^oa ■ ea ■ p g + a ■ p g a ■ ea ■ p^o ) + a^a ■ ea^C^ 



\u(Pg), 



(90) 
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which equivalently could be rewritten as: 



PHP**) 



a • A + a • P^oa ■ ea ■ B\ + a ■ p 9ii a • ea ■ B 2 + D^a^a ■ ea^C. 



00 
21 



U(pg), (91) 



where 



A H = y-2.1 

BJ* = V~ 3 



Pg^ U 12 + °23 ~~ °22 J + ^x| ^ 21 22 ~~ 23 



Y (o) ^(o) 



B£ = V- 3 



X3 - v 11 ^ 21 y^^ c ^ v 12 ^ 23 



*(o) 
y 24 



-3 



(92) 



Replacing tt(p^o)a • p^o by m^ou(p^o) and cr • p g u(pg) by m^, and using e • p^o = Q.(f9~T]) be simplified 



to: 



~ pu(p^o) Aa • e + £?i<r • eu • p^o + B 2 a • p^cr • e + D\a ■ p g a ■ ea ■ p^o + D 3 a^a ■ ea^C^} u(p g ), 

(93) 

where 



B 1 ^V- 3 (c®+2C® + C®)m* , 
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'12 I- ^21 J ""x" 3 

D 2 = V- 3 (cfi + C§) m- g , D 1 = V-" ( - C; - ( '. 



3 / M b ) _ 

12 u 23 



(94) 



Utilizing values of C's calculated in 
A ~ 0(1(T 3 ), Si ~ O(10~ u )GeV- 2 , B 2 ~ O(10~ 14 )Gey- 2 , £>i ~ 0(1(T 25 )GW- 2 ; D 3 ~ V" 3 ; (95) 



£ |.M| 2 ~/ 6 Tr(a 

l v5 SDins \ 



Act • e + f?i<7 • ea ■ p^o + -B2<7 ■ Pg& ■ e + -Di<7 • PgO" ■ • P^o + -^3°> cr • ea^C^ 



g and X3 spins 



x a -p g 



Aa • e + -£>i<7 • ecr • p^,o + B 2 a ■ p g a ■ e + D\a ■ p g a • ea ■ p^o + D^a^a ■ ea^C. 



'21 



(96) 



which at: 



yields: 



P 1 c H" m li 



(97) 



I 6 
256' 



:m~ c 



6mg 2 (B 1 + .Dim^) 2 + {8^! + 16D 3 C 2A + m g + y/tj B x + 8B 2 + (5 + V^J flpj) }' 



(98) 
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Figure 8: Three-body Gluino decay into Gravitino 



in the rest frame of the gluino. 

Incorporating results of ([M|) in equation (j98|) . one gets 



]T |M| 2 ~0.2/ 6 A 2 m? 

g and s P ms 



(99) 



Now, using standard two-body decay results (See [29]), the decay width T is given by the following 
expression: 

E^and x« spins |>[| 2 (m|-m|o) 



167r/im~ 



(100) 



Using result of |)5]); D 1 ~ O(10- 25 )Gey~ 2 and m g ~ V _ Im p GeV, ~ m 5 + GeV, two 



body decay width is given as: 



r 



(0.2)/ 6 L> 2 m? O(10- 2 )/ 6 .10- 54 .m/(2mg./ 2 ^' 



"a- v 



io- 3 / 6 G e y 



167rm~ 

Considering / 2 < 10 -8 as calculated above, T < 10~ 27 GeV. Life time of gluino is given as: 

h W~ u Jsec 10~ 22 



(101) 



T W- 3 f 6 GeV f e 



n 

-sec > 10 sec 



(102) 



Gluino(g) decays into Goldstino(G') 

We now consider the three-body decay of the gluino into Goldstino and a quark and anti-quark: 
g —7- G + q + q. The Feynman diagrams corresponding to this particular decay are shown in Fig. 8. 

The gluino- (anti-)quark-squark vertex will again be given by (|148p . The gravitino-quark-squark 
vertex would come from a term of type ip^qLQL, which in M = 1 gauged supergravity lagrangian of 
[12] is given by: 

(103) 



l K 



9iJ ( d^A J ) xV%^ - -e^ (DjW) X VVv + h.c. 
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r 18 ^/ 6 > ~ O(10 2 )V"§ ~ V-f for V ~ 10 5 



,K A 2 A 2 \j K A 2 A 2 



Figure 9: The Goldstino-quark-squark vertex 



From [30] . the gravitino field can be decomposed into the spin-^ Goldstino field G via: 

i>v = Pv + °uG, G = -^cr M W 



(104) 



p u being a spin-| field. Hence, the Goldstino-content of (|103p . using o v a ii o v = —2a u , is given by: 



29JJX 1 (d^A 1 ) o»G + |ef (DjW) X 'G + h.c. 



Now, utilizing: 



(105) 



9A 2 A 2 ~ V 9 and e 2 D A2 W 
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A 2 ^A 2 +V 3 



one obtains: 



( * • + (2?^ W) x'6 + h.c. 



7 tog 



^ 2 G + V~§x / ) A>0' -I- V-^A->\'G ~ V--.4,\ V; 



for m W2 ~ O(10) AfeV, m G ~ and V ~ 10 5 

The physical Goldstino-quark-squark vertex will be given as: 
For this particular case:- 
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(106) 
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Limits of integration as given in [TTJ are: 
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Figure 10: Diagrams contributing to one-loop Gluino decay into Goldstino and gluon 



Sl3max = V^Uls ~ S 2 3,Sl3min = 



where m g = V 3777,3 ~ 10 GeV, = 0. Utilizing the values of vertex elements calculated above and 



from (|79|) . we can calculate decay width for Gluino in this particular case. 



r(<7 -> XnQiQj 
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The life time of gluino is given as: 



n 



10" 34 Jsec 



10" 



r 10- 21 / 2 GeF / 2 



> 10 4 sec 



We now consider the two-body decay of the gluino into a Goldstino and a gluon: 
The matrix element for the above will be given by: 
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As discussed in [11] , the Wilson coefficients corresponding to Gluino-Goldstino- Gluon coupling 
do not change much upon RG evolution to EW scale. The matrix amplitude in equation (|110|) can be 
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written as: 



p v -3 

+f 2 v- 4 
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(*>) 
21 



111) 



which equivalently could be rewritten as: 



a ■ A + a ■ PqO ■ ea ■ B± + a ■ p g ^a ■ ea ■ B2 + D^a^a ■ ea^C. 
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where 
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This time around replacing u(p G )a ■ p G by and a ■ Pgu(pg) by rrigu(pg), f 1 11 2 j) can be rewritten as: 



P u (Pg) \ A 2<? ■ e + -830- • p g a ■ e + D 2 a ■ p- g a ■ ea -p G + D^a^a • eo^og) ) u(p~ g ), (114) 



where 



a 2 ee v- 3 c^;b 3 ee v- 4 M 5 (e^ + or 2 j ) 
d 2 ^v- 4 (-eg) + eg)) ;J D4 = v- 4 



(115) 



Results of various C's functions required for this particular case similar to the ones given in [TT 



arc: 



Cg) = eg } ~ 0(1); 

eg) ~ o(io) x io- 22 e e y- 
eg } ~ o(io) x io- 22 c e y- 
eg } ~ o(io) x io- 22 c e y- 

cg } = eg } ~ o(io) x io- 22 e e y- 2 ; 
eg) = o(io) x io- 22 c e y- 2 . 



(116) 
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Utilizing dUSD, one gets: A 2 = 0(1) x 10~ 15 , B 3 = O(l) x lO-^GeV -1 , D 2 = 0(1) x W^GeY" 2 , 
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which at: 



m ~g/ 2 ,P G =P 2 G =P% 



m-g 
' " 2v/3' 



yields: 



jT4 2 



m; 



f 4 A 2 2 m s 



So, using results from [29], the decay width comes out to be equal to: 

V - \M\ 2 
r = ^ andG f pins ' ' ~ 0{lQ- i yfA 2 2 m- q ~ 10- 20 / 4 OeF. 
167rnmg 

since / 2 < 10 -8 as calculated above, T < 10 -36 GeV. Life time of gluino is given as: 

h 1(T 34 J sec 1(T 5 _ n 



-sec > 10 sec 
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r W- 20 f 4 GeV f 4 
3.2 R-Parity Violating Neutralino Decays 

In the following we will be discussing the decay width of tree level decay diagrams of neutralino 
mediated by sleptons/squarks, Higgs and gauge Bosons. 
Decays mediated via squarks and sleptons 

Now, we evaluate all possible neutralino decays which are mediated by squarks/sleptons and also 
involve R-parity violating vertices as shown in Fig. 12. The life time calculation requires the evaluation 
of various matrix amplitudes corresponding to various decay channels. In particular, to evaluate the 
contribution of neutralino-squark-quark/neutralino-slepton-lepton, we will explicitly work out squark- 
quark-gaugino vertex replacing the gaugino by —^3 with mass half of that of the gluino and also the 
squark-quark-Higgsino vertex replacing the Higgsino by fx% and then add these contributions. To 
calculate the contribution of various interaction vertices shown in Fig. 12 in the context of N = 1 
gauged supergravity, one needs to consider the following terms of gauged supergravity action of Wess 
and Bagger [12] . 



£ = aYM9 a jX a x J ^ g + wax 
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a ■ dx' + r%a ■ da L X j + (p aL Ko ■ a L - c.c. 
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e 2 



+ ~— (ViDjW) xV + h.c, 



:i22) 



where, X a corresponds to the components of a killing isometry vector. From here, one notes that 
X a = —QiK^pbjQa, where a = B, Q a = 2na' J Tb i*ui a A P- f where P_ is a harmonic zero-form on Eg 
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Figure 11: Feynman diagrams for the R-parity violating decays of Neutralino 



taking value +1 on S# and —1 on a (Tib) - a being a holomorphic isometric involution as part of the 
Calabi-Yau orientifold - and / £ Hi (T B ) = coker (h 2 (CY 3 ) 4 Hl(T B )). Also, 



PiD^ W = fl^fljW + (d^-lT) VF + diKD jW + djKDiW - (diKdjK) W - T k tj D k W. 
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These expressions will be useful for all subsequent calculations. 

Higgsino-fermion/anitifermion-sfermion + Higgsino-fermion/anitifermion-sfermion + 
Fermion-sfermion-fermion/anitifermion vertices 

• As discussed above, the neutralino-lepton-slepton vertex corresponding to Fig. 11(a), is given by 
considering the contribution of gaugino-slepton-lepton vertex with a small admixture of higgsino- 
lepton-slepton vertex as given in equation (|66p . Since neutralino is of Majorana nature, in two- 
component notation, the contribution of Higgsino-lepton-slepton vertex in gauged supergravity 
action of Wess and Bagger [12] , is given by 



K_ 

e 2 



{v Zl D Ai w) x Zi x c 1 +igiA i r 



1 



a-d X c 1 +T^ A a-dA lX c 1 + 7 {d M Ka ■ A 1 - c.c.) X 



X Z is SU(2)l higgsino, x Al corresponds to SU(2)l electron and A\ corresponds to left-handed 
squark and gj A = 0. 

Strictly speaking, SU(2) EW symmetry gets spontaneously broken for Higgsino-lepton-slepton 
vertex, however the effective Lagrangian respects SU(2) symmetry. Therefore to calculate the 
contribution of same, basic idea is to generate a term of the type IlIlH^Hl wherein x 1l an d 
Hl are respectively the SU(2)l quark and Higgs doublets, II is also an SU(2)l doublet and iT£ 
is SU (2) l Higgsino doublet . After spontaneous breaking of the E W symmetry when H° in Hl 
acquires a non-zero vev (H°), this term geneiates:(H°) H^IlIl- Now, in terms of undiagonalized 
basis, consider 

ViD- ai W = (did- ai K) W + diKD- ai W + d- ai KDiW - {diKda.K) W, 
where ai,Zj correspond to undiagonalized moduli fields. 

2 

1. Considering a\ — > a\ + V sM p , using equations (|29j) , (|33|) and further picking up the com- 

2 

ponent linear in z\ as well as linear in fluctuation [a\ — V~sM p ), we see that: 
ef ((did ai K) W + diKD ai W + d ax KD % W - (diKd ai K) W) tx ai ~ V"Szi(ai - V-5M p ).(131) 
As shown in ([35]) 

efViD Ai W ~ 0(l)efViD- ai W, 
Higgsino-lepton-slepton vertex will hence be given as: 



e^ViD Al W X Zi X cAl 



efViD^Wx 2 ^ 1 



A\ corresponds to left-handed slepton. The contribution of physical Higgsino-lepton-slepton 
vertex after giving VEV to Zj will be given as : 



V"18 < Zi > 



3 

V"2. 



(132) 
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2. The interaction vertex corresponding to gaugino-lepton-slepton vertex is given by 

9YM9jf B x * B X J ^g- 

Utilizing (HMD andX B = -6«|u 7 Q Ts , «!^7 ~ ~ V'^,Qt b ~ V^(2Tra') 2 f,g Y A49T B 

_ 2 

—V sai and 

for iJj —7- (2^) ~ V 36 (in M p = 1 units), the dominant contribution to the physical gaugino- 
lepton-slepton vertex is proportional to : 

cwi „ ^M^^-v-j/ ^ 7 ( 10 -»y-3) „ / f v -s} . (133 ) 

3. Keeping in mind the fact that physical neutralino eigenstate X3 is largely a gaugino with a 
small admixture of Higgsino in our set up, by adding the contribution of ([1320 and ([1330 
as according to ([660 . the physical neutralino-lepton-slepton vertex is given as: 

CxSfefe : / (V~§) <Uix V 1 + /V-WiAgX* ~ / (v"§) x%Il for V ~ 10 5 . (134) 

The quark- quark- slepton vertex corresponding to Fig. 11(a) in gauged supergravity action is 

K 

given by the following term: [T) a D Aa W') X^X ■ In terms of undiagonalized basis, consider 

V a - 2 D a4 W = (d a - 2 d a4 K) W + d a - 2 KD a4 W + d a4 KD a - 2 W - (d a - 2 Kd a4 K) w. 
Using equations ([290 and (1330 : 

ef {(d a - 2 d a4 K) W + 8 a - 2 KD a4 W + d a4 KD a - 2 W - (d a - 2 Kd a4 K) W) ~ V"§ + V*5oi. (135) 
From ([1350 . picking up the component of fluctuations linear in a± and from ([350 

e^A,^ ~ 0(l)efv a - 2 D a4 W, 
one gets the contribution of quark-slepton-quark vertex as : 

e^V M D A4 W X A h A4 ~ ef^D^^X^ 4 ~ (vsMi) X^V 4 - (136) 
The physical quark-slepton-quark vertex is given by 

C d i l ' LUL V ~ (10) _ xiH ~ V - f for V ~ 10 5 . (137) 

KA2A c 2 K AlAl k AiA4 

The contribution of Higgsino-quark-squark relevant to the Neutralino- quark- squark vertex of 
Fig. 11(b) is given by: 



K 

e 2 



— {V l D Ai W)x L X M +i9l M X I 



a ■ 3 X M + T% Ai a ■ dAAX M + \ {d M Ko ■ A 4 - c.c.) X M 
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Working in diagonalized basis, ffi^ 4 =0. 
Now consider 

ViD a4 W = 8id ai W + (did ai K) W + diKD a4 W + d a4 KDiW - (diKd a4 K) W - T k iai D k W, 
where a 4 , Z{ correspond to undiagonalized moduli fields. 

1. Considering a 4 — > a 4 + V~~vM p and utilizing (|125p . (|126[) and (|130p . one gets: 

rg* ~ vf ~ v%6a.4,r% ai ~ v-i +vifo 4) r^ a4 ~ v§ + v§fo 4) r^ 4 ~ v~& + v^5a 4 , 

23 7 

r^ a4 ~V-86+Vs5a 4 . (138) 
Next, using equations (f29|) . ([33]) . (|138j) . one obtains: 

(y%a A D Zi W + r^ 4 L> ai W + T% a4 D a2 W + I^A^W + T% aA D aA W) ~ (v~w + vS«5a 4 ) (139) 
and 

f ((c^if) + ^iO^W + d a4 iOW - (diKd a4 K) W) x l X aA ~ V"M + V^Wx" 4 - (140) 

Utilizing equation fll39f) and (I140p . picking up the component linear in fluctuation <5a 4 , 
contribution of Higgsino-quark-squark vertex will be 

efv i D A4 W X Zl X M ~ ef AD^W^X- 44 ~ V^<U 4 xV 4 
and physical Higgsino-quark-squark takes the form as below : 



71 



19 71 



11 



(10)~"3"Vra ~V~T5. (141) 



2. The interaction vertex corresponding to gaugino-quark-squark vertex is given by 

V *B c A i \ 

For a Kahler moduli space, 9_a 4 t b = i- e the gaugino-quark-squark vertex does not con- 
tribute to this particular vertex. 

3. In view of above, the physical neutralino-quark-squark vertex as according to equation (f66l) 
is given as: 

. />-^a 4 x a V] ~ fV-£ X ° 3 d c L d R (142) 

The lepton-squark- quark vertex corresponding to Fig. 11(b) is given by considering the COntri- 
if , K A 

bution (Vfa D_a 2 W) x c 1 X" 42 + {PAi ^Ai W) X Al X c 2 ■ I n terms of undiagonalized basis 

V^D^W = (da.da.K) W + d^KD a2 W + da.KDa.W - (da.Kda.K) W 

and 

V ai D a - 2 W = (d ai d a - 2 K) W + d ai KD a - 2 W + d a - 2 KD ai W - (d ai Kd a - 2 K) w 
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Utilizing — > 04 + V s and equations (f29j) , ([33]) , on solving, 

ef (da-^K) W + d dl KD a2 W + ^iTDoiW - (da-.Kda.K) W X dl X a2 ~ (V"T + V-^a 4 ) x a V 2 

ef (5 ai 5 a - 2 K) W + d ai KD^W + 8 a - 2 KD ai W - (d ai Kd a - 2 K) W X ai X a ~ 2 ~ (v~T + V-*<5a 4 ) X ai X ai - (143) 

From (|143p . one gets the contribution of lepton-squark-quark vertex as : 

ef D^Z^WVV 2 + efv ai D a - 2 Wx ai X a z ~ (v"*^) x a V 2 + (v"*^) X ai X ai (144) 

and the physical lepton-squark-quark vertex is given by 

, s 5 
C l L d R u L + ~y-a (145) 



K AxA-l K AiAi K AaAa \ K AiAi K A 2 A2 K AaAa 



The Neutralino-quark-squark vertex of Fig. 11(c) is given as considering the contribution of 



K 



a-d + r** A & ■ dA 2X cM + \ (d M Ka • A 2 - c.c.) 



X 1 is SU(2)l higgsino, x Ai corresponds to SU(2)l quark and A 2 corresponds to left-handed 

squark. and gjA 2 = 0- 

In terms of undiagonalized basis, consider 

ViDu 2 W = (did- a2 K) W + diKD U2 W + d- a2 KD t W - (diKd- a2 K) W, 

where a 2 , Z{ correspond to undiagonalized moduli fields. 

1. Considering a 2 — > a 2 + V~zM p , using equations (|29 |) . ([33]) and picking up the component 
linear in z\ as well as fluctuations of a 2 i.e {a 2 — V~»M P ) , one gets 

ef ((did- a2 K) W + diKD- a2 W + d^KDiW - (diKd- a2 K) W) ~ (V~f ) Zi 5a 2 . (146) 
Using pS]) : 



the Higgsino-quark-squark vertex will be given as 

e-ViD A2 W X Zl X c -e-ViDa-.Wx^x" 2 ~ V~~ Z^A^x" 5 

^2 correspond to left- handed squark and the contribution of physical Higgsino-quark- 
squark vertex will be given as : 



2(1 



C Hlu L u L _; ^ : _ y- t ( U7 ) 

K ZiZi K A2A 2 K A 2 A2 



4 
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2. As discussed for gaugino lepton-slepton vertex in Figure (a), the gaugino- quark-squark 
vertex is given by 

9ym9jb*X* x 

. Utilizing [El and X B = -UkImQ b ,k\h7 ~ ^,9ym ~ V~^,Qb ~ V^{2na') 2 f, This 
time, 

9YM9B5.2 -> -V~*5a 2 - 

Now, 

9ym9t b A\ ~ 0(l)g Y M9T B a 2 ->■ V" 45^2 • 
The dominant contribution to the physical gaugino-quark-squark vertex is proportional to 

— — J 

C A S ^«1 _ V / ~ / fv~i) . (148) 



^2^2 J ( y K A 2 A 2 

Following equation no (|66p , the neutralino-quark squark vertex will be given by : 
c xWu L ^ c x- g u h u L + j c Hiu L u L ^ j / v -*j sA 2 x I X cA2 + f (V~i) <W 2 X A A 5 ~ /V-t x 0nLUL. (149) 

• The lepton-squark-quark corresponding to Fig. 11(c) is given by considering the contribution of 

k . 

(T>£ Dj\, a W) x c 1 X^ 4, i n gauged supergravity action. In terms of undiagonalized basis, 

V a - x D aA W = (da.da.K) W + da.KDa.W + d a ,KD Sl W~ (da.Kda.K) W. 

Utilizing equations (|2"9"j) and one obtains: 

ef {{da-^K) W + da-.KD^W + datKD^W - {d dl Kd ai K) W) ~ V-Ts + V-Uo2. (150) 
From (|150p . picking up the component of fluctuations linear in a 2 and from (|35p using: 

one gets the contribution of lepton-squark-quark vertex as : 

ef V Ai D M W X ^X M ~ ef P^VU^V 4 ~ (v~5M 2 ) X^V 4 - (151) 

The physical lepton-squark-quark vertex is given by 

_ i 

C d i lLlL V ~ (10)~tV~s ~ V~t for V ~ 10 5 . (152) 

Now, working in two component notation, in order to calculate the decay width, we will be using the 
decay width formula as given by H.Dreiner et al. in |27j . Because of Majorana nature of neutralino, 
we are considering both right handed as well as left handed incoming neutralino though the dominant 
contribution occurs in case of vertices corresponding to right handed neutralino. Henceforth, we will 
be using incoming right handed neutralino in the matrix amplitude calculation. 
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diPd) 
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V 3d c L l L u L 



X°s(Pi) 



d.R 




d.R 



lipl) 



V U L d R u L 



X3(Pi) 



UR 



UR 




U{p u ) 



d(Pd) 



V 3l LX d l x l 
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Guided by their notations, right handed incoming and left handed outgoing are denoted by wave 
function y\ = y ] {Pi, Aj), x\ = x*(p u A;), xt = x^(p u , X u ), and x\ = x^(p d , X d ), respectively . Using the 
numerical estimate of vertices as calculated in set of equations (|134p . (|137p . (|142p . (|145p . (|149p . (|152p 
for all three feynman diagrams, the corresponding contributions to the decay amplitude is: 



3 . . 5 



M > n; ' n; 'M (K _ t ,),_ m J ^M4 (153) 



ii 



4 5 



iM 3 = fV-i V-t I —2 - y\x\x\x\. (154) 

Neglecting all of the final state fermion masses, one can express kinematic variables in the form as 
given below: 

zi = 2pi-ki/m 2 x o = 2Ei/m x o,z d = 2 Pi -k d /m^o = 2E d /m x o,z u = 2p i -A; 11 /m^o = 2E u /m x o(155) 
and the total matrix amplitude can be rewritten as: 

M = c x y\x\x\x\ + c 2 y\x\x\x\ + c 3 y\x ] u x\x\, (156) 

2 1 

where for m^o ~ V'^mz and mi, ~ m,- ~ m,-,, ~ V^mz 
x 3 2 L d R L 2 



ci = fV-'h-l/[m? L -m? x o(i- Zl )] = -fV-lv-l/ml[V-Vl(l-z l )] 
c 2 = /V-sV-i/[m| H - m 2 x o(l - z d )} = -f\r^V~^ /m| [V - V*(l - z rf )] 

= / v -| v -|/ [m 2 _ m 2 = _/v-tv-|/mi[V-Vt(l-z u )] (157) 

A 3 2 

As explained in [27J, by applying Fierz identity, one can reduce the number of terms: 

M = (c 1 - c 3 )y\x}xlx d + (c 2 - c 3 )y\x\x\xi. (158) 

and 

\M\ 2 = |ci - c 3 fy\x\xpyix\x\x d x u + |c 2 - c 3 \ 2 y\x\x d y i x\x 1[ u x u x l 

-2Re[(ci - c 3 )(c2 - c%)y\x\x x x u x\x\x d y^ , (159) 

Summing over the fermion spins by using following identities, 

J> Q (p, S )*+(p, S ) = +p-a a p,^2x^(p,s)x fS (p,s) = +p-^,Y,y j '"(P,s)y tS (p,s) = +p-^ 

s s s 

^2y a (p,s)y^(p,s) = +p-<J a p (160) 

s 

one obtain: 

|7W| 2 = |ci - CspTr^-api-o-jTr^-afcu-a] + |c 2 - c 3 \ 2 Tr[k d -api-a]Tr[k u -akra] 

spins 

-2Re[(ci - c 3 )(c* 2 - cl)Tr[ki-Wk u -ak d -api-a]] . (161) 
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Applying the trace formulae and using 

+ 2ki-k d = (1 - z u )m 2 x o, +2k r k u = (1 - z d )m^o, +2k d -k u = (l-Z/)m£o. (162) 



^ |A4| 2 = 4\a - c 3 \ 2 pi-kk d -k u + 4\c 2 - c 3 \ 2 pi-k d k r k u 



spins 



4Re[(ci - c 3 )(c* 2 - c* 3 )](k r k u pi-k d + pi-ki k d -k u - k r k d pi-k u ) 

ci\ 2 zi(l - zi) + \c 2 \ 2 z d (l - z d ) + \c 3 \ 2 z u {l - z u ) 
-2Re[ci^](l - zi){l - z d ) - 2Re[ci^](l - z t )(l - z u ) 
-2Re[c 2 c* 3 }{\ - z d ){\ - z u )} , (163) 



The differential decay rate follows: 

(pT N c m o / i 



^QEM 2 ). (164) 



dzidz d 2 8 

spins 

where N c = 3 and the kinematic limits are 

< zi < 1, 1 - zi < z d < 1. (165) 



Adopting the technique as used in [27] , the total decay width for the above Feynman diagrams is given 
as: 

jV c m 5 

r = (K*| 2 + \c' 2 *\ 2 + \c' 3 *\ 2 - Re[<£c£ + ci*4* + c' 2 V 3 *]) , (166) 

where c'*, i = 1, 2, 3 are obtained from Cj's by integrating zi, z d in the aforementioned limits. Therefore, 
on simplifying, the result comes out to be 

3.vfm 5 3 / 2 .V-Tf 

r 2 - — « io- 4 vi _ S m3 / 2 ~ 1(T 18 / 2 (167) 

2« • 3.7r3V§m 4 3 5 

2 

In dilute flux approximation, considering f 2 < 10~ 8 , life time of neutralino: 

h 10~ 34 Jsec 10" 34 Jsec 



r io- 18 / 2 io- 26 G e y 

3.3 Two-Body Slepton/Squark Decays 

The decay width for l/q — > l/q + ip^ in our set up is given by: 



> O(10 1 )sec (168) 



m~,_ 

i/9 



r^/^i/^^j-^-V-^, (169) 



implying a lifetime of around 10 25,5 s, satisfying the BBN constraints. 
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3.4 Three-Body Slepton Decays 

As explained in [3T], the following set of effective operators are relevant to three-body slepton decays: 



0i 




e* 






e* 


o 3 


= i'kPg-g c pi 


•e* 


04 


= H'hl-t* Pv 


• G c 




= l' h 7 • Pv 7 • e 


* p r • G c 


6 


= i l' h 7 • Pv 7 ' 


e* p y • G c 


o 7 


= *Jft7-i>G7- 


e* p v ■ G c 


o 8 


= iT h e*-G c 




o 9 


= H'hJ-Pv e* 


■ G c 


Cio 




Pv e* • G c 



(170) 

where is the polarization of gauge boson V . Notice that for an on-shell gravitino in the final state, 
O3 = by using the gravitino equation of motion. 

Using the following volume-suppression factors in various vertices in Fig. 14(a)-(d) above: 

l-G-l vertex : 9AiAi = 1; 

K A 1 A 1 



I- I'-V vertex : ^] V ~ HnV 



I — I — V vertex 



fv- 



4K AiAi 104 ' 



l-l'-G vertex: = 1, (171) 



and the operators of (|170p . the matrix elements for three-body slepton decay [diagrams (a)- (d)] are 

11 

fV"T8 1 
My ~ ^ (201 -0 6 ) 

fV~ 2 1 

^ ~ /n4 A/f 2 7 2 (02 + 03) 

10 4 M p m Z Q + m v - S13 - s 2 3 
~ 2 K 3 (0 4 -09)-4(0 6 + 7 -0io)+m 2 y O 8 ] 

4 S13 S X3 

~ i TS° 8 ' (172) 

where X3 denotes the lightest neutralino in our set up, and 

S12 = (Pg + Pv) 2 , «13 = (p G +py) 2 , S23 = ipv +Pvf ■ (173) 
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Notice that S12 + S13 + S23 = m~ + m 2 -, + m 2 z . 
The differential decay width is 



dT 



(2tt) 3 32m? 



\M\ 2 dml 3 dm. 



23 s 



where M = M^ + M b v + My + M v . 



The sum of the matrix elements can be written as 

M(l h ^ I'VG) = M v + M b v + M c v + M v = M l O l , 



i=i. ..10 



where the Mi can be read off from Eqs. (|172h above as: 



M 
M 

Mi 

The squared matrix element is 
M{l h -> iVG)' 2 



/v- 



1,5 



4,9 



M„S 2 3 



/y- 2 



2,3 



M p( m 3/2 ~ S 13 - S23) 



10/v- 



M p (si3 - 



; ^6,7. 



4,9 



X3^ 



10/v- 



+ 



10 



/v- 



??? 



11 
9 



si 3 - m 



\3 



10 4 



i=1...10 



ij=1...10 i,j=1...10 



(174) 



(175) 



(176) 



(177) 



All the Mi are real except for Ms, which has both real and imaginary components. The only non- 
zero contributions to the last term in Eq. (I177j) come from 9m (MiMg) (i < 8) and Qm(MsM : j) 
[j > 8). Expressions for O^i, 0\ e p Of^ (i < 8) and Cg™ (j > 8), and the corresponding contributions 
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to J J dsizds23^r- are given as under: are: 

4 



O 



2,2 



03,3 
04,4 



3m 2 - 



G m v 



i i 

m G m i 



(pg-piY 



m G ~V G - P\ +P G -PV 
2 



m|my + m-my - 2m v p G ■ pj - (p G ■ p v ) + 2p G ■ p v p~ l ■ p v - (pj ■ p v ) 



2m? 
I 



ds-\ 



13 

2 /™2 



'3/2 



o 

mi 
I 



ds 



23" 



|A4i| 2 



3mi 



1 Ci 1 / ~ ii \^ i 



10 



-16.5 



G -P G -Pi+P G -Pv 



2 2 

ra^ray 



(Pl-Pv)' 



- (p g ■ Pi) 2 + 2 P G - P\P G ■ PV - (P G ■ Pv) 2 + m 2 G (mf + m v - 2p f ■ p v ^j 



2m? 
( 



ds 



2 

'3/2 



13 



in 



T , \M 2 \ 2 2 ,2 
ds 2 3 5 

2 m~ 

V I 



fv- 



V"tMj, ~ l(T 21 M p ; 



10 4 



= 0; 



3m G my l. 



2 2 

m G m v 



(Pg-PvY 



m G m v - m v p G -pj + 2 (p G ■ p v ) + p G ■ pv {3m v - 2pj ■ p v ) 



2m i 



ds 



2 

'3/2 



13 



ds 2 3- 



(lO/V-l) 2 v!m p ~ i(T 18 M p ; 



49 



05,5 



07,7 



O, 



8,8 



0, 



9,9 



010, 



10 



3m% 



2 2 
m G m i 



(Pa-PiY 



,2 _2 
G 



ds 



06,6 = " 



rn 

2m? 

4 3/2 

4 



3m 2 - 



m v - m v p G ■ p\ - 4 (p^ • py) + p^ • p v (-3m z + 4p ; ~ • p v ) 



l , |-M 5 f <%,5 
OS23- 



13 / Q 

2 m- 



2 

V^Mp ~ 10 _16 - 5 M P ; 



m G m ^ ~ {PG 'PvY 



m 



2 ^,2 



m v - m v p G -pj-4(p G - p v ) +Pq-pv (-3m y + 4p,~ • py) 



/•2m? /• 
/ 2 / 



2 

m~ 



"3/2 
4 



23" 



|A^ 6 | 2 06,6 



m 



~ ( 10/V~3 j V"^Mp ~ 10~ 35 Mp; 



3m^m v 



m % m v ~ (Pg'PvY 



{m^my + 2m| (pg • py) 2 + 4p 6 ■ p v - p G ■ P\ m Q my + 4 (pg • py) 
-m G p G ■ p v (m v - 2p~i ■ p v ) | , 



/ 2 dsi3 / 

Jm 3/2 J " 



o 

; 



23" 



|^ 7 | 2 07,7 



(io/v-1) 



2 8 



16.6 



M p ; 



f2m? 



G 



p G ■ Pi +p G -pv 2m G m v + (p G ■ p v ) ' 



/ 2 / 2 

m 3/2 m ' ' 

4 



23" 



|^ 8 | 2 08,8 



10 



-27.5 



.1/ 



2m G m v + (p^ - py)' 



m^m 2 / - m 2 ^ ■ p r - 2 (p^ • py) 2 - p G ■ p v (m v - 2pj ■ p v ) 

9,9 



,2m? ,m? |A^ 9 | 2 99 / 3 \2 5 1S 

/ ds 13 ds 2 3 91 3 9 ' 9 ~ (lO/y-s) V9M p ~ 10- 18 M p ; 

•/ m? ,„ im?- m ~ 



'3/2 
4 



V 



3m G m v 



2m G m v + (p^ • py)' 



{m G m v - 2m 2 G (p G ■ pyf - 4 (p 6 • py) 



+p G ■ pi 



2 2 

-m G m v 



/•2m? rm 

/ ds 13 / 



+ 4 (p<5 • py ) 2 + m 2 G p G ■ p v (3m v - 2p f • py ) | , 
^ (lO/V-i) 



i , 1-MiorOio, 
a«23 5 

2 m~ 

V I 



VaM p ~ 10" 16 - 6 M p ; 



(179) 
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(p g ■ Pi) + Pg ■ PPg -Pv + m^l m f - pi ■ p v 



-m 



]4 + m 2 z p G ■ pj - Pq ■ p v pj ■ p v + fa ■ p v ) ' 



\ ds 13 / ds 2 3 

m 3/2 +m V ^ m V 



m l Re (MiM 2 ) Of 2 



>- ~ 1(T 26 M 



/m 3/2+"4 
Of 4 = 



my 



P' 



G 



2 2 

m Q m 



!-(p G -Pi) ~(P G -Pv) +m 2 G (m v -p rPv ) 



+P G • Pv- (-"if + Pi ■ Pv) + P G ■ Pi (-4 + Pg " Pv + ■ Pv) 

/•2m? /•) 

/ ctois/ 

Jmi,„ Jm 



3/2 
■ire 



rfS23 5 ~ / V 9 



2 10" 16 - 5 M, 



p> 



Oft = Oft = Oft = Ofto = 
Oft = 

— / (r> x. ■ nrt-' 2 

3m 



rn 2 G rny + 2 (mf - m v ) p G ■ p v 



y{(Pg-PiY 

G 

+2 (p G ■ Pi) 3 {m 2 v - pi ■ pv) - 2p G ■ Pl [(p G ■ pv) 2 (mf - Pl ■ pv) 
+m 2 G (mf - pj ■ p v ) (m v - pj ■ p v ) + m 2 Q p G ■ p v p\ ■ Pv 
+m 2 G mf (p G ■ pv) 2 - 2p G ■ p v (mf - pj ■ p v ) 
+m| (-mfm v + fa ■ p v ) 2 ) } , 



ds 



13 



'3/2 



ds 



23" 



Re (M 2 M 5 )O t 2 c 



~ 10- 37 M p ; 



Oft — Oft — Oft — Oft — 

/-/-ire /-/-ire /r\xe /r\xe /r\xe /r\xe r. 

^3,5 — ^3,6 — u 3,7 — u 3,8 ~ u 3,9 ~ u 3,10 ~ u 

0; 
_8 

m 



w 



G 

2m? 
I 



ds 



l G my - (p G ■ p v ) [m v +p G -p V ~ Pj ■ Pv] , 
Re (M 4 M 6 ) Of 6 



13 

2 l m 2 

'3/2 



m? 
I 



ds 



23" 



~ f z V~^M p ~ 10~ 42 M 3 



m G 



i - P G " + Pg " Pv- m|my - (pg • py) 



• /m 3/2 



9 

mi 
I 



ds 



23" 



Re (7W4VW7) e>4 e 7 



my 



(/v-l) 



V"~M P ~ 10 



-3 
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3m G m v 



-m\ (mj - 2p G ■ pjj Pq-Pv+Pg' Pv 



+m 2 ^m v (my - pj ■ p v ) - (p G ■ pyf (m v + pj ■ p v ) 



J mi ,„ J mi. 



Re (M 4 M 8 ) Of 8 



'3/2 



Sm^rriy 



rn^my + m 2 G m v p G ■ p v + m|my (p G • py) 2 - 2 (p G ■ py) ' 



-m v p G ■ pj (m G m v + (p G ■ p v ) J - (p G ■ p v ) (m v - 2pj ■ p v ) 



,2m f ri 
/ 2 / 

Jm i/2 Jm 



n Re(M A M 9 )Ol% 
ds 2 3 5 - 



m 



(lO/V-i) 



2 V3M P ~ KT 20 M P ; 



|m^,my - m| (p G • py) 2 - 3p G • p v + p G • pj -m|my + 3 (p G • py) 
+m 2 G p G ■ p v (3m v - 2p t ■ p v ) } , 

/ dsi 3 / (isss 3 ' ~ (lO/V-2-J V'^M p ~ l(T 24 M p 



'3/2 

/-/-ire /n re ^O re 

^5,7 — ^5,8 — u 5,9 







5,10 



I 





3m| 



m|m y - (p G -p v y 



3m G - 2p G • p-A p G ■ p v + 2 (p G • py) + m G (m y - p f • p v ) 



,2m? , 

/ 2 ^ 13 / 

Jmf,,„ J rr 



mi 
l 



'13 / dS23" 

3/2 



Re (.MeMy) 0£V 



(lO/V"f) 2 V~f M p ~ 10 -33 M P ; 



3m| 



m^,m y - m| (p G • py) 2 + (p G • py) 1 



-p G ■ pj (m^rriy + (p G ■ Pv) ) ~ ™ G P G ■ p v (m v - 2p f ■ p v ) 



2m? 
I 



'3/2 



mi 
I 



dsi 3 / ds 2 3 



Re(M 6 M 8 )Ol% 



~ 10" 27 Mp; 



3m G 



m v 



( 3m l ~ 2p G " Pf) • - 3 (p G • py) v 



(Pg " Pv) 2 (my - 3pj ■ p v ) + 



2 2/2 

m G m v (m v 



PfPv) 



r 2mi 



'3/2 



9 

m~ 



dS23 



Re(M 6 M 9 )Ol% 



m~ 
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3m% 



+m^my (m v - p~ x ■ p v ) + m| (p G ■ pyf (m v - pj ■ p v ) 
Re(M 6 Mi )Ol c l0 



G 
c2m? 



\ ds 13 

J m 2 ,„ J rr, 



% Z/2 



ds 2 3- 



m~ 



(lO/V-i) y~f M p ~ 10" 33 M P ; 



3m 2 ~m 



2 



m G m vPG • PV - ( m | - 2 Pg • Pi) (PG • Pv) 3 - 2 (PG • PvY 



+m^m v (m v - pj ■ p v ) + m Q (p G ■ p v ) (m v - pj ■ p v ) 



G 

,2m 2 rmf Re {M 7 M 8 ) Of s 

/ ds 13 / ds 23 7 3 B; 7 ' 8 ~ 10~ 17 - 5 M P ; 

Jm i /2 Jm v m i 



V 



3m 



111 



-%m v - m| (p G ■ pvY + (p G • py) 1 



G L 



-P G " Pf + (P(3 " f'V') 2 ) - m G P G ■ PV (my - 2pr ■ pv) 
/ dsi 3 / d«23 3 - ~ (10/ V 2 J 



:;m |^r {™S™v + m 4 G m v p G ■ p v + mjm y (p G ■ pv\ 

G z 

2m| (p G ■ pvY ~ 4 (p g " P^) 5 



m G m v + m G m v (p G ■ p v ) - 4 (p G ■ p v ) 



+m| (p G ■ pvY ( 3 ™v ~ 2 P\ ■ Pv) \ , 



/ dsi3 / ds 23 3 — 7, ~ (lO/V-5) V9M P ~ 10- 16 - 6 M ? 

J ml,- J rr 2 m ~ 



'3/2 



3m G m v 



2m G m v + (p G • pv-) 2 [my + p G • pv - P\ ■ pv] , 



2 1 



dsi3 



'3/2 



TO? 



^^W„ 10 -, iM>; 



my 



3m 2 -m 



2 



2m 2 G m v + (p G -py)' 



2p G • p,-) P G • - 2 (p G • pv) 2 + m|( 



GV m v-PfPv) 



/• 2m r , r 

/ 2 / 



^ MWiol^o 18 . 5 



3m 



/// 



G L 

2m 2 
1 



G 



- Pg • Pf + Pg " P^ 2r n G mv + (p G ■ Pv) 



dsn 



'3/2 



dS23 



Re(M 9 M 10 )Ol c 10 
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(m| -Pg-Pi + Pg- Pv) { - (Pf • Pvf + (p g • Pv)(Pi ■ Pv) 



m 



-m 2 v {p G ■ pjf + (p G • pi){p G ■ p v ){pi -Pv-P G -Pv)}, 



2m? 
I 



'3/2 



13 



033 



2 

171- 



ds 23 — 8J h8 ~ l(T 44 M p ; 



"W^ (m G -Pg'Pi+Pg- Pv) (mjmy - (jpj ■ pv) )m G - m v (p G ■ pjf 



-(P G • Pvf Pi ■ PV + {P G • Pl){P G • Pv){m 2 V + Pi • Pv)] , 
,2m? ,r 

/ 2 ^ / 



im 



3,8 

im 
5,8 



'3/2 

oi m 8 = o 



ds 23 3 ~ 10 M M p ; 

2 m~ 



?>m?~ 



r {(pg • pi ~ m2 G )Pi ■ Pv m l ~ p g ■ pi(pg ■ pv? 



+P G ■ Pv 



/-/-iim 
u 6,8 

Hence, 



/•2m? 




l' m \ 


/• 




1 ds 23 


Jm 3/2 




m v 


/-/-iim 
U 7,8 ~ 


/-/-iim 
" ^8,9 


/-/-iim 
— ^8,10 



Pg • PP 1 ! + (Pj • Pv ~ 2mf)m 2 G + (p G ■ pjf } , 
9m (M 5 MS) °58 



l(T 38 M p ; 



Tll^l' + G + V 



~ / ' dsi 3 / 

J ml,,. Jn 



m x , Re(7W 4 A^8)OI e 8 

a«23 3 - 

2 m; 
v l 



~ 10 _15 M p , 



G 



(187) 



(188) 



implying that the associated life-time of this decay is about 10 28 s, which respects the BBN constraints. 

4 Gravitino Decays 

In this section, we discuss two-body and three-body gravitino decays and show that the corresponding 
lifetimes are either of the order or more than the age of the universe 

4.1 Two-Body Gravitino Decays 

We discuss the decays of the gravitino into neutrino and gauge bosons as well as the light Higgs and 
neutrinos. 
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v 

Figure 15: Two-body gravitino decay: W — t v + 7 




Figure 16: Gaugino-(sneutrino) -neutrino vertex 

In M = 1 SUGRA the gaugino-(sneutrino)-neutrino vertex, Fig. 16, is given by: 

vT B — — 

0(a x ) - term m — -== 2 77^ • (189) 

Similarly the gravitino-gauge-boson-(z>) — v- vertex in the M = 1 SUGRA Lagrangian is given by the 

0(ai)-term in gYM9 ai T B X TB Z° /A^^j^ipu, which is: 9Y * l9a i' fB ^ v 10 7" . Using [32], one sees 

W k a 1 a 1 ) 

that the lifetime of the decay ipu — > 7 + v e is given by: 

93 / 11 \ 2 

£ rn . o / t r t- \ 



W/ s 1 A Ml) A m 3/2 



2 / 



10 4 



1 ( (V-l)V Hj/ 2 f M Xl -M X2 \ 2 . 2 „ 2 „ /V~t\ 2 non , 

id iT? — sm ^vycos (9 W -77^4- • (190) 



64^ ^ v ^ / M P 2 V m Ai m A2 y y 10 4 _ 

We have a gaugino mass degeneracy (up to 0(1) factors) implying: M Al ~ M\ 2 . One hence gets a 
very small decay width and an extremely enhanced lifetime. 
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+ A 



Figure 17: Two-body gravitino decay: ip^ ^ Z° + v 




Figure 18: Two-body gravitino decay: W ~> h + v 



Again, using (|189p and [32], the decay width for ifi„ — > Z + v is given by: 



i /(^ 



64vr V <2i) 7 Ml 



1 m 3/2 



3/2, 



10 4 



Mi 



x O(l) + O(l) 



1 + sinj3cosf3 



Mf 



3/2 



M| 



MaA^iZ 2 



^ x 0(1)M P . 



64vr x 10 8 



x OH) 



(191) 



This, for V ~ 10 6 yields r ~ 10 21 s. The mass- like insertion is given by nv~ h , which is the coefficient 
of vh in the M = 1 SUGRA action. To estimate this, we note that near \z±\ ~ \z2\ ~ 0.7V36, |ai| ~ 
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Figure 19: Gravitino-gaugino-Higgs vertex 

_2 _1 13 _U_ 

V s , \a2\ ~ V 9 , 1 0,3 1 ~ V is , 1 0,4 1 ~ V 9 , the potential can be approximated by 



V ~ e K G Tsfs \D Ts W\ 2 

z 72 e -2n°{vol{i:s)+^ 2 )^ V ^ +fi3Z 2 



V + (Vs + (ai 2 Vi8 + Q13V9 + ai 4 V9)ai + (anV 9 a x + a^Vis + ai 3 V9 + ai4Vs)|ai| + ^z 2 ) 2 - (Vis + fi 3 z 2 



(192) 



The coefficient of the OUz- 0.6V3e)(ai - V'a)) I K z±Zx K A^M in (fT ^ 1 turns out to be: 



m% h ~ Vi5m§ /2 x 1(T 16 . (193) 

The gravitino-gaugino-Higgs vertex: will be given by the following term in the N = 1 SUGRA 
action: 

gYMD T '0 At 7 M A = 0, on shell gravitino. (194) 

Hence, the second diagram involving a gaugino NLSP, does not contribute. However, if the gaugino is 
replaced by a neutralino, then the Higgsino-component of this neutralino will contribute via the term: 

g Zl zA Z ^^ V ^Hl + h.c. (195) 
in J\f = 1 SUGRA action. This yields a vertex: 

jjz 1 z£0 ^ j . (196) 

Further, the Higgsino-(z> e )-^ e vertex will be given by 0{a\ — z{) term in e~a" (V ai D Zi W) v e lH^ : which 
after assigning (zi) a value of ~ Vse yields: 

/Vl(i/ e ). (197) 
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(a) 



e jL t 



e-kR 




00 



(c) 



Figure 20: Three-body gravitino decays involving $, p \jk coupling 



Hence, using (fT96l) . (fTWD and [32]: 

1 



3/2 



384?r I M p 2 



+ svnpcosp 



We) 



^ 3— ^777 



~ l(T 3 V~ 6+ i x 10~ 32 i\4 



- 10 _60 M 



v~io 5 



(198) 



which yields a lifetime of around 10 s. 

4.2 Three-Body R-parity Violating Gravitino Decays 

In this subsection we will be considering gravitino decays involving three types of R-parity violating 
vertices that appear in R-parity violating superpotentials 



W# p - -V ,/,/-,/-/ /'.I + \' ijk LiQjD c k + X" jk UfDjDl + mHLi. 
Decays involving Af/& coupling 

• Using the same approach as discussed to calculate R-parity violating interaction term in the 
previous (sub)section(s), the vn — ZjL — £kR vertex corresponding to Fig. 20(a) is given by considering 
the contribution of following term: 

K_ 

e 2 



C = — (V Al D A3 W)x Al X M , 



where 



V ai D Wi W = 8 ai d as W + (d ai d a3 K) W + d ai KD a , 3 W + 8 a3 KD ai W- (d ai Kd a3 K) w - v k aiaA D k w, 

where 03, Z{ correspond to undiagonalized moduli fields. By expanding above in the fluctuations linear 

_ 2 

in a\ — > a\ + V 9 M p , on simplifying 



K_ 

e 2 



V ai D a3 W ~ V~-^ ai 



Utilizing above, the Dn — eji — e k R vertex will be given by 

efv Al D A . A W X Al X M ~ eTV ai D a3 W X Al X M ~ V-%5A lX Al X M 
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and physical vertex will be given as 

_19 _19 
V 72 V 72 7 c 

C u lL e ]L e kR ^_ ~ V~3,for V ~ 10 5 . (199) 

l K - k - k - ^10^ 

^AlAl^AlAl^AzAz 

Similarly, ejL — v%l — ?kR vertex of Fig. 20(b) is same as vn — ejL — e-kR vertex corresponding to 
Fig. 20(a) and given as: 

(je jL u iL e kR ^ QUiLejLe k R ^ y~l (200) 

The ekR — vf R — ejL vertex corresponding to Fig. 20(c) comes from: 

{v M D Al w) x A h Al - 



K 

e 2 



In terms of undiagonalized basis, 

V^D ai W = {da.da.K) W + d a - 1 KD ai W + da.KDa.W - (da.Kda.K) W. 

13 

By expanding above in the fluctuations linear in 03 — > as + V'^Mp by utilizing equations (J2"9"j) and 
(f33|h on simplifying 

K_ 

^-Va.D^W-V'Uas, 
implying that the lepton-squark-quark vertex is given by: 



(201) 



and physical vertex will be given as: 



_i _ 1 

C S kR u? L e jL V ~ -^=L ~ y-I , for V ~ 10 5 . (202) 

V K AiAi K AiAi K A 3 A 3 

Now, the matrix amplitude for all three Feynman diagrams corresponding to Fig. 20 will be given as: 

\M\ 2 = \M a + M b + M c \ 2 (203) 

The analytical results for the full matrix amplitude summed over spins in terms of pure and cross 
terms are given in [33]. Utilizing their results, we will estimate matrix amplitude for all for both pure 
and cross terms to calculate decay width for the process ^V - ^ e j^fc in our set up. Strictly speaking, we 
will be neglecting fermion mass as compared to gravitino and squark masses. Introducing kinematic 
variables 

2p{vi) -p(ej) = (1 - z ek )ml /2 , 2p(ej) ■ p(e k ) = (1 - z Vi )m% /2 , 2p(z/j) ■ p(e k ) = (1 - z ej )my 2 . (204) 
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In view of above, we can express the following definitions in terms of kinematic variables as given 
below: 

m lj = (p(y%) + P(ej)) 2 ~ M"i)-P(ej) = (1 - *e k )ml/2 
m)k = (P(ej) + p(e k )) 2 ~ 2p{e j ).p{e k ) = (1 - z Vi )m\ /2 

m ifc = (p{ v i) + P(. e k)f ~ 2p(i/i).p(e fc ) = (1 - z ej )m 2 /2 , (205) 
utilizing the form of expressions given in the appendix of [33] 

|M Q 



2 l^LyJ J , 2 2,2 w 2 2 2 \ 

= 3 M 2 (m 2 fc -m 2 ) 2 (m3 / 2 " m ' fc + ) ^ " m< * " ™ efc ) 

m ?jfc , (206) 



After simplifying using (I230p . we have 



imp - - ic^^r 4 

|W " ! 3M 2 ((l-z ei W 2 -m?) 2 ^ (1 Z " K 4 )J ' 



~- { - . j tfr.fl -*„.)( — )] (207) 

Here we have neglected gravitino mass as compared to sfermion mass. With same steps of similar 
procedure, we have 



\M h 



1 (c eLVLeR ) 2 
3M 2 (m 2 fc -m? 



2 ) i 2 2,2 \i 2 2 2 \ 

~ F27^2 ^2-^2 ( m 3/2 " m ik + ™ ej ) Kfc " m Vi - m £k , 

^e jL ) 



{ml /2 + m 2 ik -m 2 e .) 2 

Am 2 , m ' lk 

- 2 L [ Z ej (l-z ej ^ (208) 



\2 1 ^ / 2 2 , 2 w 2 2 2 

1 3M 2 (m 2 -m? ) 2 3/2 J fc J 1 3 

VI \ Cut)' 



(m 2 /2 + m| - m 2 J 2 ^ 
4m 2 /2 



1 (C g « p fl ei ) 2 -z 2 



2Re{M a Mh 



1 {C^iL^jL^kR _(je L v L e R \ 



/ 2 2 2 2 2 2 \ 

( m ik m jk ~ m 3/2 m e k ~ m vt m ej ) 



( 2,2 2 2 \ 1 / 2 , 2 2 

(m 3/2 + m efe - m v . - m e ) - — 5— (m 3/2 + m jk - m h 

Zm 3/2 



60 



(ml /2 + m 2 k 



,1/2 2 2 v - 

+ n( m a -m v . -m e .)(m, 



2\'"iJ ""Vi ■■"ejJV b jk 



/• 2 2 2 \ / 2 

2 

m e fc /2 2 2 \2 , 2 2 2 



2 2 



2 \2 



/ 2 2 

-(m ifc - m„ 



2 {C" iLe i LekR .C eLUL£R ) 

3 M^.m? m? 



(1 - Z„.)(l - Z ej )(-1 - ^e fc + 22„. + 2z ej - Z Vi .Z ek ) 



2Re(M b Ml) = 



1 (C* lVlCr .C* RV R eL ) 

3M U m fk - m l jL )( m2 ij - m l kR ) 



'■plV'Hk 
l 3/2 



I 2 2 2 2 2 2 \ 

Kj% - m 3/2 m ^ - m ei m e fc , 



/ 2 2 2 2\ i /2,2 2\ 

(w. i/2 + m,. - m e . - m e J - ^(m 3/2 + m ik - m t . 5 



/ 2 , 2 2\\, ± /2 2 2 \ / 2 2 2 n 

(m 3/2 + m i3 - m e J j + -(m^ - m„. - m e .){m jk - m e . - m e J 



(m ifc - m„. -m e J - — ±( 



2 2 2 n2 

— m e — m„ 



-(m 2 k -m 2 



m 



et/2 2 2 \2 , o 2 2 2 



2(C 



e-LV L e-R 



_C^R u R e L) 



3 M 2 m 2 ~ m 2 - 



(1 - Z e )(1 - Z e J(-l - + 2z ej + 2z ek - Z ej .Z ek ), 



2Re{M a Ml) = 



UJv iL ej L e k R _(je.RV R e L \ 



ZkR> 



(( 2 i 2 2 2 \ 

(m 3/2 + m e . - m„. - m e J - ^ 



(mfjin 



ij jk 



m 3/2 m e 



/■ 2 , 2 2 \ 

(rn z/2 + m jk -m v .) 



( 2 i 2 2 

(m 3/2 + m i3 - m efc 



*" l 3/2 

, \ , 1 / 2 2 2 u 2 2 2n 

) J + ^ m ij ~ m ^ ~ m e 3 K m jk ~ m e . - m ek ) 



(m 



m 



ik 



m 



ek 



m„ 



("4 



m 



Vi i 2 

m 



jk 



m 



2 \2 
e k ' 



m 2 )+2m 2 m 2 m, 



''k 



2 (C° iLe i LekR .C^ RU R eL ) 

3 m^4J4~ l 



(1 - z„.)(l - z e J(-l - + 22„. + 2z efe - z„..z, 
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Utilizing the results from (|206j) - (|212p . (|203j) one gets the following form: 



2 1 I rier.PT en\2 ^ 



1 (C i>iLe i LekR ) 2 Z„. 1 (C^ LULeR ) 2 



\M\ 2 = ^ ' K(l " **)(-?■)] + 2 M 1 - ^e,)(^)] 



3 M>4 — -*a 4 " 3M 2 ( 



1 (C g « p ^) 2 _ < 
+ 3 M 2 m 2 m 2 C 1 " ^iX 1 " Z ej(-1 " ^e fc + 2^ + 2z ej - Z v% .Z eh ) 

Pi V iL gjL 

+ o M 2 m 2 m 2 C 1 " ^)(1 " *ej(-l " ^ + 2^ + 2z efc - Z ej .Z e J + 

a" TTi—2 2 -C 1 " - ^)(-l - z e . + 2^. + 2z efe - z n .z ek ). (213) 

3 M pl m o iL m~ e . L 



The differential decay rate follows: 



spins 



Putting the result of |M| from above, 

d 2 r iV c m 3/2 1 (c^ e ^ e *fl) 2 ^ z 2 1 (C £i ^ e «) 2 ^ 2, 



z Vi (l - z u J(^)] + - , \ \z e Al - z e M^)} 



d Zej dz ek 2V ^ M 2 m 4^ l*V m/V 4 yj 3M 2 (m y 2 L^ e,A 4 

+ 3^KJ 2K(1 " ZeJ( ~ )] 

2 (Qv iL e iL e kR . (je L u L e R \ 

+ 3 M 2 m 2 m 2 C 1 " - - * efc + 2^ + 2z ej - z Vi .z eh ) 

pi ViL ej L 

2 (c^ LVLeR ■ C gRp R eL ) 

+ o M 2 m 2 m 2 (! " ^)(1 " " Z ^ + 2 ^ + 2Z ^ ~ z er z e k ) + 

* M pl m y lL m e jL 

2 (C i>iLe i LekR ■ C iRp R eL ) 

o- M 2 rn 2 m 2 " ^JC 1 " ^)(-! " % + 2 ^ + 2z e k ~ ^.Ze k )\ (215) 

Using the same approach as given in fjl65f) . 

< z e . < 1,1 ~Zj < z £k < 1. (216) 

and using the numerical estimates of masses, m%/2 ~ V~ 2 M p ,m~. L = rrv~ L = m 1 kR ~ V^ m 3/2 hi our 
set up, after integrating, decay width reduces to 



N c m 7 3/2 



(2 9 .3.120)7r 3 .M 2 .V 2 7^ /2 



t(jhLej L e kR \2 _|_ (jje L v L e R ^2 _|_ (^je R v R e L -^ 



_|_^ y^Q^iLejhekR . (J^ L v L e R ^ _|_ ^Qie L u L e R _ (je R u R e L ^ _|_ ^QUiL^jLekR . Q^R v R e ty\^ (217) 
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ViL d jL d kR 




Utilizing the set of results given in equation no ()199|) - ()202p . decay width simplifies to 

N c m 7 3 ^ 2 



Life time will be given as 



" (29.3.120)vr3.M2.V 2 m4 /2 (V 2) 
—11 

1 V 2 .7Tlo/n 

5-^ ~ 10- 45 - 5 Gel/; for V ~ 10 5 . (218) 

10 6 M* ' v ; 



- ~ ~ O(10 21 ) S ec. (219) 

r 10- 45 - 5 GeF v ; y J 



Decays involving A^- fc coupling 

• The Dl — djL — dkR vertex corresponding to figure 3(a) and II — ul — d c L vertex calculated in 
section 4 get identified with same set of moduli space superfields and given by similar interaction 
vertex. Therefore contribution of vl — djL — dkR vertex is same as II — ul — d L given in (1137|) . 

c v L d jL d kR ^ V~%v L dj L d kR , for V ~ 10 5 (220) 

• The djL — — dkR vertex corresponding to figure 3(b) and ul — II — d c L vertex calculated in 
section 4 get identified with same set of moduli space superfields and given by similar interaction 
vertex. Therefore contribution of djL — v%l — dkR vertex is same as ul — II — d c L given in (|152p 

C d jL u iL d hR ^ y-f,for y ^ 10 5 (221) 

• This time, dkR — v£ R — djL vertex corresponding to Fig. 21(c) and d R — II — ul vertex calculated 
in section 4 get identified with same set of moduli space superfields and given by similar interaction 
vertex . Hence, contribution of dkR — vf R — djL vertex is a same as d R — II — ul given in (1145 j) 

QdkRPtRdjL ^ V~ t , for V ~ 10 5 . (222) 

The analytical result of matrix amplitude for all three Feynman diagrams corresponding to Fig. 21 
is same as calculated in Fig. 20 by replacing m^ L — > ray iL (m^), me jL — > m d L ( m u jL ) an< ^ fne hR — > 

m-, fro; ) and C" ihe i hekR -)• C ViLd i hdkR C^ LUlLCkR -> C d ^ LUiLdhR . C^ kRp iR e ^ L ->• C dkRp iR d ^ L . Doine 

(7,1. E> V fill. D / C) 
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so and using same numerical estimates of masses, rn. 3 / 2 ~ V 2 M p ,mf, iL = m~ = m| ~ V 2 m3/2 
our set up, after integrating, decay width reduces to 



m 



N c m 7 3/2 



(29.3.120)vr 3 .M2.V 2 m4 /2 

3 
+ 4 



fQV iL d jL d kR \2 _|_ (Qd jL v lL d kR ^2 _|_ fQd kR v? L dj L \2 
^(jVihdjhdhii . (jdj L v iL d kR \ _|_ ^(jd L u L d kR _ (jd kR Pf L d jL ^ _|_ ^(jUihdjhd kR . (jd kR u^ R d jL ^ 



. (223) 



Utilizing the set of results given in equation no f|199[) - (|202p . decay width simplifies to 



r 



^ m 3/2 



(2 9 .3.120)7r 3 .M 2 .V 2 m 4 /2 



-(V-t) 



16 Q 

X V 3 .m 



10 6 M 2 , 

pi 



3/2 ~ 10~ 44 GeF; for V ~ 10 5 



(224) 



Life time will be given as 



h 10 34 Jsec ^. 9fK 
- O(10 20 )sec 



T 10~ u GeV 



(225) 



Decays mediated via X'-- k coupling 

• The um — djR — d c kL vertex corresponding to Fig. 22(a) arises from: 

K_ 



In terms of undiagonalized basis, 

V a - 4 D a4 W = (d a - 4 d a4 K) W + d a - 4 KD a4 W + d a4 KD a - 4 W - (d a - 4 Kd a4 K) w. 

By expanding above in the fluctuations linear in 04 — > + V _ ~s~ M p by utilizing equations (I29p and 
([33|h on simplifying one obtains: 
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implying that the contribution of quark-squark-quark vertex is given as under: 



(226) 



and the physical vertex will be given as following: 



C ^ndl L „ 6 === „ „ y- § u iR d jR d% L , for V ~ 10 5 . (227) 

^4^4 

• Similarly, the dj R — d kL — Ui R vertex of Fig. 21(b) is same as Ui R — dj R — d kL vertex corresponding 
to Fig. 21(a) and is given as under: 

c d ]R di L u m ^ c u lR d jL d% L ^ v-Md jR d c kL u iR , for V ~ 10 5 . (228) 

Again, introducing kinematic variable 

2p(ui) -p(dj) = (1 - Zd k )ml /2 , 2p{dj) ■ p{d k ) = (1 - z Ui )ml /2 , 2p{m) ■ p(d k ) = (1 - z d] )m\ j2 . (229) 

and 

m ij = (P( u i) +P(dj)) 2 ~ 2p(ui).p(dj) = (1 - z dk )m\ /2 
m )k = (P(dj) +p{dk)f ~ 2p(dj).p(d k ) = (1 - 2 Ui ) m 3/2 

m L = (P(^) +P(4)) 2 ~ 2p(ui).p{d k ) = (1 - z dj )m3 /2 , (230) 

Here also, utilizing the form of expressions given in the appendix of [33] for A^- fe coupling, the full 
analytical result of the squared amplitude summed over the spins for the gravitino decay reaction 

ifin Uidjdk is the sum of the following squared amplitude and interference teramas: 



({m\ l2 + mf fc - m 2 u ) 2 \ \ ( C u lR d jR di R \2 z f 

3/2 A ? — — ~ rncfg 4 j k(i-^)(f)l (23i) 



2Re{M a Ml) 



4iV c ! (C*h*hr<%h)2 
3 M2(m? fc -m| j 2 ^ 3 / 2 " + m ^ )(mife " " ^ 

o M 2 ( 2 d K- 1 " ^ (4 > ( 232 ) 

2iV c ! (C^ iRd j Rd kR ■ C^ RUiRd kR) 



3 M2 (m 2 fe - m?^ ) (m 2 fc - m 2 - r ] 



I 2 2 2 2 2 2 

{m ik m jk - m 3/2 m dk - m u .m d . 



05 



(/ 2 ■ 2 2 2 \ ^ / 2 , 2 2 \ 

(m 3/2 + m 4 - m Ui - m d .) - ^-\ m 3/2 + m jk ~ m Ul ) 



1 



2 



(m 3/2 + m ifc - m dj )j + -{m^ - m u . - m ej )(m jk - m~ e . - m £k 

m ^ Til 2 

/ 2 2 2 \ "Vj / 2 2 2 \2 e j / 2 2 2 \2 

(™ifc - m Vi -m ek )- -^-(rn jk - m e . - m efe ) - — (m ifc - m„, - m e J 

2 

m efc / 2 2 2 n2 i0 22 2 
2~ - m ^ - ^e, ) + m ej m e fe 

2 (C iiiRd 3 Rd kR. ■ C d i RUiRd kR) 

o ^72-^ 2 ( X ~ ~ - z d k + 2z u t + 2z d . - z Ui .z dk ). (233) 



Decay width will be given as: 

N r m 7 



"3/2 

1 - 



(2 9 .3.120)7r 3 .M2.V 2 m4 /2 
Utilizing the set of results given in equation no ()227p - ()228|) , decay width simplifies to 



(234) 



Nrjnl, 

-(V-!t) 



, c" t 3/2 , 43. 



(2 9 .3.120)7r 3 .M2.V 2 m4 /2 



7;i 



X V is .m 3 



106 M 2 



3/2 ~ 10- 42 GeV; for V ~ 10 5 , (235) 



and therefore the life time will be given by: 

h W~ u Jsec 



O(10 18 )sec. (236) 



r W- 42 GeV 

5 Relic abundance of Gravitino dark matter candidate 

As discussed in section 4, explicit calculations yielding the gravitino-decay life time of the same order 
or greater than age of the universe, justifies taking Gravitino to be a viable dark matter candidate in 
our setup. Keeping in mind the fact that relic density of dark matter particle should be within the 
limits provided by recent WMAP observations and other direct and indirect experiments, this section 
is devoted to calculating the relic-abundance of the dark matter particle. Assuming that re-heating 
temperature will be low enough to produce an appropriate amount of dark matter in case of heavy 
gravitino, we focus on the case of gravitinos produced in the decays of co-NLSPs and show that the 
same are produced in sufficient numbers to constitute all of non-baryonic dark matter. The mass scales 
and life time estimates of sleptons and neutralino discussed in section 2 and 3 manifestly indicate the 
same to be valid co-NLSP's which freeze out with appropriate thermal relic density before decaying 
and then eventually decay into the gravitino. Therefore, gravitino then inherits much of the relic 
density of the neutralino/slepton. 

5.1 Neutralino density calculation 

The number density of cold relic density (CDM) from the early universe depends sensitively on the 
annihilation cross section of such particles. Amongst the various approaches to calculate the thermal 
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cross section given in the literature, we rely on the partial- wave expansion approach used in |13j to cal- 
culate the annihilation cross section for each possible process. Since the //-split SUSY set up discussed 
in section 2 includes only neutral light Higgs boson; heavy Higgs boson and similarly superpartners of 
neutral particles (i.e only Neutralino), we consider possible annihilation channels which proceed via 
these particles only. The important channels that we will be discussing are: X3X3 ~^ hh, X3X3 ~~ ► ZZ, 

xlxl^ff- 

The general procedure which is being followed in partial wave expansion approach, is given as 
follows. As given in [34] . in QFT, the general cross-section is given as 

<™ M0 i = j dLIPS\M i2 

where 



(TV U0 i = J dLIPS\M\ z (237) 

d 3 p 1 



dLIPS={2^5\p 1+P2 -Y,P 3 )X{ j J^L ■ {2: ^ 
product of 

velocity (ci>M0l) is given as p~3l |3l] 



^2tto ' 

The thermally-averaged product of the neutralino pair-annihilation cross section and their relative 



f d 3 p\d 3 p 2 M0I e~ El l T e~ E2 l T 

^M l)(T) = J d3pid3p2e - El/Te -E 2/ T > ( 239 ) 

where p\ = (E\,pi) and p 2 = (E 2 ,p 2 ) are the 4-momenta of the two colliding particles, and T is 
the temperature of the bath. Given the complexity of the general cross section, it is difficult to solve 
analytically. Alternatively, one finds a way to get the solution by expressing {<jvm0i} in terms of 
x = jj. For this, one defines: 

w(s) = 1 j dLIPS\M\ 2 = E 1 E 2 av M0 i- (240) 

Incorporating the value of (TVM0I in equation (1239[) . 

Jd 3 p 1 d 3 p 2 w(s)e- E ^ T e- E2 / T 
{aVM » l){T) = jd^p 2 E lE2 e- E ^e- E ^ - (241) 

By defining change of variables discussed in [34] i.e expressing momentum and energy in terms of x, 
carrying over the integration in terms of x, ((TVM01) takes the form given below: 



(o"«M0l) = — 2 



m x 



3 

w - - (2w -w')x + 0(x 2 ) 



a + bx + 0(x 2 ). (242) 



2 

The coefficients a and b summed over all possible final states fif 2 are defined in |13| . and given as: 

= cd (4M| ~ + m h?) v hh "fifa 
hh 



b= J2 ce ( 4M I ~ ^ m /i + m f'^ 2 ) V hh \ h hh + a hh 
hh I 



mi + mi {mi — m" f 
h h 1 v h h 



-3 + -v~* + 



4 Kb \ 2M? 8M$ 



The analytic expressions of a and b are given for the s, t and u channels of various possible 
annihilation processes are given in [13]. For the paper to be self-contained, we directly quote the 
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Figure 23: Feynman diagrams for X3X3 ~~ ^ hh y i a s— channel Higgs exchange and t-channel exchange. 



forms of analytical expressions and utilizing the same, we calculate the numerical estimates of a and 
b for all kinematically possible annihilation processes in our set up. 



The idea is to first calculate the required vertices corresponding to different annihilation processes 
in the context of M = 1 gauged supergravity action and then use their estimates to calculate a and b 
coefficients . Following the same formalism as used in section 3 and 4, utilizing the gauged supergravity 
action of Wess and Bagger, we obtain the numerical estimates of required vertices. 

The physical eigenstates of neutralino mass matrix in the context of gauged supergravity action 
are given as: 



A2 



/A u + 



mass ~ V 36 M v > m a . 

y 2 



CP 



(i H9 + H9 „ 35 w 

U H ; mass ~ V seM^ras: 

v2 2 

A + / (Hi + H$) ; mass ~ V~*M P > mz . 



(243) 



Utilizing (M) . 



(7x3X1^ 



G X»H»h + fC l 
j£,\°\°h _|_ p(j\ H h _|_ JqH°H ) 



(244) 



where -ff 



" ~- 2 is physical Higgsino and the physical light Higgs is defined as h = 1 2 Q 



Working in the sublocus, where position moduli z 1 and Z2 are considered to be equivalent, for notational simplifica- 
tion, we will write Higgsino superfield \ v^' Zl+2 ' 2 - ) ~ X Zi ■ 
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Higgsino-higgsino- Higgs vertex 



K 



C = ^- (V Zl D Zl W) xf^xf + igz^xt 1 



7 • d X f 1 + rg^ 7 • dZixf + i 0^7 • Zi ~ cc.) x? 



Xl/Xji corresponds to left-/right-handed components of the Higgsino. In terms of undiagonalized 
basis, 

V zl D Zi W = (d Zi d Zi W) + {d Zi d Zi K) W + d Zi KD Zi W + d Zi KD Zi W - (d Zi Kd Zi K) W + T« z p K W. 

Since SU(2)l symmetry gets spontaneously broken for Higgsino-Higgsino-Higgs vertex, therefore the 
idea is to expand above term quadratic in Zi such that one of the two ZiS acquires a VEV. Utilizing 
Zi — > Zi + VseMp and thereafter solving with the help of equations (J29]) and ([33]), one has 

K_ 

e -Lv Zi D z% W ~ V-t (zi)6zi 
Followingd35]), e^V 2i D Zl W ~ 0{l)e^V Zr D Zi W , implying: 

e^VzVztWxfxj? ~ (v-"(2i>) ^afxf ~ (v-^(2<)) ^ixfxf • (245) 



Using x^ ~ V rri 3/2 f° r the Higgsino mass and m 3 / 2 = V 2 M p , one obtains: 

37 

— — 



v-3xf ^xf ; 



_ 4 

-V-Ixf^xf- (246) 
Incorporating results of (|245p and (|246|) . the physical Higgsino-Higgsino- Higgs vertex will be given as 

= (V^) + (V- 2 + V-i +V~i)l ~ Vi (247) 



c H°H°h ^ 1 



Gaugino-Higgsino-Higgs vertex 

C = gYM9BZ l X B XL^L + c.c. (248) 
The physical Gaugino-Higgsino-Higgs vertex works out to yield: 

c hS°x L = V^f\ _ „ f( w s V -l) „ /v-f . (249) 



69 



Gaugino-Gaugino-Higgs vertex 



C = iXL^-AKz^Zi - c.c.)A L , 



(250) 



where A^ corresponds to gaugino. Here also, the aforementioned vertex does not preserve SU(2)l 

2 

symmetry - one has to obtain the term linear in (zi). In terms of undiagonalized basis, d Zi K ~ V _ s (zi), 

2 

and using dz t K ~ 0(l)d Zi K, we have: d Zi K ~ V^^(Zi), incorporating the same 



10 & V~36/iA L 4f A L ~ V~36hX L 



' M, 



~ 10 5 V"i/iA L ^A L ~ V~i/iA L A L 



l-(Px L +P\ L ) 

Mr, 



Xl 



(251) 



Higgs-Higgs-Higgs vertex 

Expanding the effective supergravity potential V = e K G TsTs \Dt s W\ 2 in fluctuations: Zj = 5Zj + 
V»M p , the contribution of the term cubic in Higgses is of the order V~3er < zi >. Utilizing the 
same the physical Higgs-Higgs-Higgs vertex will be given as: 



C 



hhh 



1 



V-^(Zi)bZ l bZnz k ~ 0(V- L )5Z l 5Z J 5Z 



Now, using the set of results given in equation no (|247p . (|249|) . (|251|) . (|330p . the contribution of 
vertices appearing in equation (|244p are as follows: 



V"36. 



(252) 



Since now we have got the estimates of coupling required to calculate partial wave coefficients for 
X3X3 ~~ ► hh annihilation process, we are in position to calculate the contribution of partial wave a and 
b coefficients in our set up just by using (and quoting verbatim below) the form of analytical results 
provided in [13J . 

s-channel Higgs-boson (h, H) exchange: 



The analytic expressions of a^h and bhh corresponding to Fig. 23 are given as: 



(h,H) 
hh 

7AKH) 
'hh 



0, 

3 
64?r 



r=h,H 



Qhhr (jx%X%r 



4 Ml - ml + iY r m r 



Expanding the summation 



(253) 
(254) 



J hh 



64 7T 



Qhhh QX%X%h 



+ 



C hhH C X 3 X°3 H 



4:m^o-m z h + iT h m h Am z ^- mjj + iT H m H 



Mi 



(255) 
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Utilizing the value of mass nih = 125GeV,iriH ~ V"m3, mo ~ V»m2 ~ V~5M and C hhh 

2 X3 3 

10^" (V -1 ), C x ° x 3 )/l ~ (7 X 3X3 H ^ y-gs from above, after simplifying, we have 

2 



64 7T 



V _1 .V"36 V _i .V"36 

+ 



M 2 - 



25 \ 2 



P 64 7T i m^j 



M 2 
p 



64 it 



x O(10)- 2y GeU~ 2 for V ~ 10 5 



(256) 



Here we assume that i\# < m^jj in our set up. 
• Neutralino (xi) exchange: 



r(x°) 
l hh 

J hh 



0, 



1 3 

1_ (c^^) 2 ^.' 



16 7T 



A L 



x |4 m 2 o (m 2 o - m/j 2 ) 2 + 4m x « (m 2 o - mh 2 ) (m x o + m x o) A 
+3 (m x o + m x o) (m x o + m x q) A hi A hj 



(257) 



(258) 



where = m^ 2 — m? — m 2 . Utilizing the values of masses given above, A^i = m^ 2 — m 2 — m 2 

X3 Xi X3 Xi 

V 2 m 2 3 , A h2 = m h 2 - m\ - m 2 ~ V 2 m% , A h3 = m h 2 - m 2 - m 2 ~V3m 2 3) and 

2 A3 X2 5 X:i X3 5 



16 7T 



(C*i x 3>) 2 ((7x¥xS,h*)J 



A L A L 



4m x o (m x o - mj + 4m x o (m x o - m h ) (m x o + m x o) A w 



+3 (m x o + m x? ) (m x o + m x o) A w A J + {CP&& h )\CP&&> h *f 



A 2 A 2 

^hl ^h2 



2\2 



4m x§ ( m x§ ~ m h) 



+4 mo (m 2 - m|) (m o + m o) A hl + 3 (m o + m o) (m o + m o) A hl A 



h2 



+ 



{C x\xlh ) 2 {C x%xlM ) 2_ 1 _ x ^ _ mi) 2 + 4m ^ (TO ^ _ m , } (mxg + m ^ Am 



*hl "/i3 



+3 (m x o + m x o) (m x o + m y o) A hl A h3 ] + (CP&& h f{CP&&> h *f 



xi 



A L A h2 



4m x o (m x o - m h 



2\2 



+4m x o (m 2 o - m\) (m x o + m x o) A h2 + 3 (m x o + m x o) (m x o + m x o) A fe2 A 



(C^X§ ft ) 2 (C X 3 X 3' ft *) 2 1 



X3 X2 

x 



A L A L 



X3 X2 ' ^ X3 

2 l^Jl _2\2 



"X^ 



/)2 



+ 



4m x o (m x o - m h ) + 4m x o (m^o - m h ) (m x o + m x o) A A2 



0y°/i\2/^v°v° /i*\2 



+3 (m x o + m x o) (m y0 + m x o) A h2 A h3 + (C*^) 2 (,7X3X3 



X3 X3' X3 X2 

,2 2 



A 2 A 2 



2\2 



4m (m — m h ) 

x 3 v X3 h/ 



+4m x o (m x o - m h ) (m x o + m x o) A h3 + 3 (m y o + m y o) (m x o + m x o) A h3 A h3 



X3 ' X3 ; V X3 



X"j,> 
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1 

16 vr 

1 

16 7T 



(/ 4 vw + (/ 4 v)- 



x O(10) _36 GeV _:i for V ~ 10 5 and / ~ 10 -4 . 



Al A 2 



+ (/ 2 V2V" 



35 . 
18 



' +/ 4 v)^ + (/ 2 vhH§; 



m v om v o 

Al A3 



2 
X2 



m v 8m v o 

A2 A3 



(259) 



Higgs (/i, H) neutralino (Xi) interference term: 



Ah,H- X °) 



"hh 



°hh 



0, 



1 3 

Sir 



16 vr 



E 



. 4 m 2 n — m 2 + i T r m r 
r=h,H \ xi 



[2 m x o (m 2 g - m 2 h ) + 3 (m x o + m x oi) A hi ] 



A hi 



(260) 



(261) 



Re 



E 

r=h,H 



(jhhr (jXXr 



4 m 2 — m 2 + i T r m 

X3 



y I I Of* 

35 



4m 2 - m\ + i T h m h 

A3 ✓ 



4m 2 o - mf^ + iT h m H 



(V).V-^M p (V).V"36M„ V36M„ 

o 1 o ~ 



m 



Expanding this, we get 



11 



2 
x° 



(262) 



L(h,H- x °) 
J hh 



1 V 36 M p 

16 7T m 2 n 

X3 



(C* ft ) 2 (2m z » Ko - m 2 ) + 3 (m x o + m^o) A ftl ) 



V3 



A 2 



+ 



(C^^) 2 (2m x o (m 2 g - m 2 ) + 3 (m x o + m x o) A h2 ) (C7^ h ) 2 (2 m x o (m 2 g - m 2 ) + 3 (m x o + m x o) A M ) 



1 VaeMp 

16 7T m 2 n 

X3 

1 V^eMp 

16 7T m 2 

X3 



"X 2 



+ 



X3 v "x u 3 



X3 



Xz> 



A 2 

(M) 2 (3m y0 ) (/V?) 2 (3mJ ) (V"i) 2 (5m y0 ) 

Xl _| X2 _| X3 

4m x? 

(V-i) 2 (5m y0 ) 

X3 

4 

X3 



4m 4 n 

x§ 



4 

x° 



x O(10)~ 26 GeV- 2 for V ~ 10 5 . 

16 7T 



(263) 



Utilizing results of equations (l256jl . ([259]) . (12631) . 

6m = ?/ ) +C ) +^ H - X ° ) ~O(10)-^eF- 2 

0. 



J hh 

a hh = a hh 



hh 
(h,H) 



. ~(x°) , 
+ Hh + «M 



(264) 
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(JxlxlZ = j c H°H°Z^ 

C xlx%Z = j c \°\°Z + f C H°H°Z^ 

C X%xiZ = C A°A°Z + p c S°HOz (2g5) 

Z Boson-Higgs-Z Boson vertex 

C 3 G TbTb X Tb X t ~ b A^A v (266) 

where corresponds to gauge boson and \i = 0, 1, 2, 3 is space time index. The required vertex will 
be accommodated by dzrd Zi G TB f g . As given in appendix B 

dz- k d Zl G TBfB ~ 0(l)%^G Tflfs V" 2 (^)5^ (267) 

and X = X b &b = — Ylmol k^^QbQt b ~ V~3. Incorporating values from above, the physical Z 
Boson-Higgs-Z Boson vertex is proportional to 

VWz~U o(io)- 8 ; 1 j 

Higgsino-Higgsino- Z Boson vertex 

£ = fe^X*' £ Im (X B iT + iD B ) x Z ' , (269) 

X 21 corresponds to Higgsino, X = X b 8b = — Ylinol ^^QbOtb corresponds to the killing isometry 
vector and the corresponding D term generated is given by: 

D B _ ^a'n\ii 7 QBV B ^ 270 ^ 

Utilizing gz x z- ~ ~ I'' -3 as given in (|124p as well as v B ~ V» and Qb ~ Vs (27ra') 2 /, yields 

value of physical Higgsino-Higgsino- Z Boson vertex as: 



(v-i)/ 
^1 ~o(io~ 5 ) 



C B ° 6 ° Z ~ V 9 J (o(10 5 )/V-i) ~ fV'Ts. (271) 
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Gaugino-gaugino-Z Boson vertex 

C = \ L t Im (X B K + iD B ) X L , 



A L 7 • A \ 6nln 7 2Tra'Q B K + 



V~~ 



B 



(272) 



A corresponds to gaugino. Again, Utilizing values of v ~ Vs and Qb ^ V3(2tv(x) /, yields value 
of physical Gaugino-gaugino-Z Boson vertex as: 



V-I8 /- /V 



(273) 



Using set of results given in equation no (|27ip . (|273p . (|268p . the contribution of vertices appearing in 
equation (|274p are as follows: 



(274) 



Again, using (and quoting verbatim) the analytical expressions needed to calculate partial wave coef- 
ficients for X3X3 ~~ * annihilation process, we will estimate the values of same in our set up. 
• Higgs boson (h, H) exchange: 



a 



(h,H) 

zz 
'zz 



0, 

3 
64?r 



r=h,H 



(jZZr (jxlxl* 



m 



2 + i T r m r 



3mz 4 — 4m^ 2 ra v o 2 + 4m v o 
^3 -^3 

m z 4 



(275) 
(276) 



Expanding the summation 

r(h,H) _ 3 



64 vr 



+ 



Am 2 -mi+iT h m h 4m z - m z H + iT H m H 

x 3 a x 3 n 



3 mz 4 — 4 mz 2 f7i x o 2 + 4 rn x o A 



(277) 



Utilizing the value of mass = 125GeV,mH ~ V 72 M p , m x o ~ V 3M p ,mz ~ 90GeU and C 



~ / 2 y- 4 ,C*3X^ „ v _ 36 from above, after simplifying, we have: 



b zz ~ 



64 vr 



/ 2 v-iv- 



4.m o 

X3 



+ 



/ 2 v-iv- 



77? 



H 



2 4 

x 3 2 

— T-"*pi ~ 
my 



64 7T \ 4.TO q 



10 \ 2 4 

f/G n 

X3 2 



rrir 



~ O(10)- 23 GeU" 2 for V ~ 10 5 . 
Here we assume that # < h in our set up. 



(278) 
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neutralino (xi) exchange: 



r(x°) 

l zz 

1 

47T 



J_ ^ \C x * x ° z \ 2 \C x >° z f 



A/.A 



| CX ;xg^|2| C X?xg^|2 m X§ 2 + | C X?X^|2| CX ° X °Z|2^XrL + | g X^ |2 1 gxgx^ 1 2 ^ + 

AziA^i A^iA^2 A^iA^ 



|C^xg^|2 [( gxgxg^|2 + | C X 3 °xi2|2| (7X X 0Z|2_^x|!_ 

A Z2 A^ 3 A^ 3 A^ 3 



(279) 



where A^j = m^ 2 

^Z2 



m 2 — m 2 . For the given neutralino mass eigenstates, Azi = mz 

X3 Xi 

V 2 m 2 , A^2 = mz 2 — m 2 n — m 2 n ~ —V 2 m 2 , Azx = mz 2 — m 2 



mz" - my - m z ~ -V*m\ , A Z3 
3 2 
values and couplings given in (|274p . ()279p reduces to 



2 







X3 



m 



(1 



X3 



2 — m 2 

x 

— 9 

2V 3 m 3 . Using above 



m" — m 2 „ ~ 

x§ x? 



Ax°) 
l zz 



1 



1 2 

47T X 3 



/ 4 v- 



V 4 m| 

2 



ft /» v 
+ 



V 4 m 4 , 



+ 



112 

3 



10 

V 3 m 



+ 



10 

V 3 m 



¥ / 4 v 
— + 



22 
9 



8 



22 
9 



— m o- — s 

Att *3 yfm 4 



O(10)- 54 GeF 



(280) 



As from [13] , the analytical expression of b zz is defined as following: 



r(x°) 
} zz 



1 3 1 

leTr^m^Al-Al. 

x|C^^| 2 |^°^| 2 bg) A|, + D§> A Zi A Zj + 4 3 > A 2 ^ + D<« A^A^ 



(281) 



where 



(i) 

ij 
(2) 



(3) 

1 1 



D 



(i) 



16 m y o 2 mz 4 (mz 2 — m x o 2 ) 



4m o 2 (rng 2 — m v o 2 ) 2 (3m/ + 4m v o 2 m v o m v o 

A3 x A3 7 L A3 Ai Xj 

+mz 2 [4 m v o 2 + 4 m v o m v o — 6 m v o (m v o + m v o )1 } , 

X3 Xi Xj X3 v Xi Xj /JJ ' 

-4m v o (m^ 2 — m v o 2 ){4m v o 3 m v o m v o + (7m v o + 3m v o + 6m v o) 

A3 v A3 7 L A3 Ai Aj v A3 Ai Aj 7 

|2m/m v (i f4m^o 2 — m^o m^o — m^o (m v o + 5m^o)]\, 

X3 X3 Xi Xj X3 v Xi Xj /JJ 

8m/ + 2mz 2 m v o 2 [4m v o 2 — 6m v o m v o — 5m v o (m v o + m v o)\ 

A3 X3 Ai Aj A3 v Xi Xj 

+4m v o 4 [2m v o 2 + 3m v o m v o + 2m v o (m v o + m v o)l 

X3 X3 Xi Xj A3 x Ai Xj 7 

+m/ [3m v o 2 + 9m v o m v o + 5m v o (m v o + m v o)l. 

X3 Xi Xj X3 v Xi Xj /J 



(282) 
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Solving above, 



n (l) n (l) n (l) n (l) n (l) _ 4,,4 4 

2 

n(2) ^(2) ,,22 10 _,(2) „(2) ,,21 in ^(2) ,,22 i 

2 22 



n(4) _ n (4) 
^11 ~ ^12 



'13 ^23 
,(3) = n (3) _ -.it « n 1 -" 



. 16 



2 22 

V^m 6 3 , D<? = ~ Vfm 6 3 , 4} ~ V 4 m 6 3 . 

2 22 



(283) 



Now, incorporating results of equations (|283p . Az 1,2,3 and ()274p in (|28ip . the same reduces to: 



~ O(10)" 14 Gey- 2 . 
Higgs (/i, ff) neutralino (x^ ) interference term: 



(284) 



Xh,H- X °) 



J zz 



1 3 

~l IT * -^ 



16 vr 



i=l 

2 _ 2\ 



E 



4m o 2 — m 2 + i T r m r 
h,H X3 



4 A 2 s 



x <j 2 m x o (m^o 2 — m^ 2 ) [—3 m^ 4 — 4m x o 3 m x o + 2 m^ 2 m x g (m x o + m x o)]. 



-A; 



-4m x o 4 (2m x o + 3m x o) + 2m^ 2 m x g 2 (5m x o + 6m x o) — m^ 4 (5m x o + 9m x o) (285) 



Utilizing results for Az 1,2,3 and (|274p . the same reduces to 



r(h,H- x °) 



1 m x° 3 



3 3 



16 7r mi 

z i=l 



E 



„ 4m v o 2 — m 2 + i T r m r 

n,iJ *3 



1 V2?71 | / C^Zfi (C*§X§ ft ) 2 C ZZH ( C X 3 X°3 H ) 2 



+ 



16 7r y4m 2 — + i 1~\ m^ 4 m 2 — m^- + i T# m# 



M„ 



(286) 



Utilizing the value of mass = 125GeV,rriH ~ V 72 M p , m x o ~ V sM p ,mz ~ 90GeU and C 



ZZft 



ZZff 



-1 ^y,/0^0 



/ V 3 , C x s x 3 V se from above, after simplifying, we have: 



b 



th,H-x°) 



ZZ 



1 V 2 m| //2 V -|. v -§ /V-i.V" 

~ i I = 1 o 



16 7r m\ 



4.m 2 

x§ 



1 



p 167rmi I 4. 



xS 



- O(10)- iU GeU- 2 for V ~ 10 5 . 



(287) 



Here we assume that T^h < m h,H m our set up. Utilizing results of equations (|278p . (|280p . (|284p . 

(EST}, 



bzz = bfP + + bf z H ~^ ^Om-^GeV^- 



(288) 



76 




Higgs-Fermion-Fermion Interaction 



(289) 



Fig. 25 involves the s-channel Higgs boson (h, H) and Z boson exchange and the t- and u-channel 
sfermion (/„) exchange. To get the estimate of thermal cross-section for this process, one needs to 
evaluate the following partial wave coefficients: 
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~{h,H) ~{Z) -(_/) + ~(h,H-f) + ~(Z-f) 

a ff + a ff + a ff + 7/ + a ff ' 



i ff = b[Y ] +W + + b { t H - f) + b ( i z - f) 



(290) 



7/ 1 7/ 1 ff 1 7/ 1 7/ • 
We have calculated below the interaction vertices corresponding to s and t channel of X3X3 — ► //■ 



.4/ 



£ 9 *fUi^X.L 



ftrf' + r^U^xf + j (d Zl K0z t - c.c.) x£ 



1 



.4/ 



(291) 



X^ 7 , 1 = 1,2 corresponding to first generation of Left-handed leptons and Left-handed quarks. 
Higgs-lepton-lepton interaction 

Utilizing z\ = bz\ + VseMp, for lepton mass ~ 0(l)MeV and m 3 / 2 = V~ 2 M p , g^a-i 

' 1 2;a, ~^ M„ '^z*-" M„ ' 



"M^ ' OIle S etS: 



2 



V 36 J , 



V 36 _ J 



(292) 
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Incorporating results of (|292|) in (|29ip . the physical Higgs-lepton-lepton vertex will be given as 



V~36 



10 1 



V36 



V 180 

~m7 



QA\ h A\ L h _ 1 
Higgs-quark-quark interaction 

Utilizing z\ = bz\ + VseM p , for quark mass m x A 2 ~ 0(5)MeV and m 3 / 2 = V _2 M p , ga 2 a 2 



(293) 



AJ2 Af p ~^ M p ' OIle § etS - 



, V 9 (ji) 



i 



M„ 



^2 



v-i 



(294) 



Incorporating results of (|294p in (j291|) . the physical Higgs-quark-quark vertex will be given as 



qA 2l A 2l \ 



v-i 

M„ 



10" 



V _ 36 



We would estimate: 

= Max (jA\ L A\hh ^ (jAihAihh 
fermion-fermion- Z Boson vertex 



V 180 



(295) 
(296) 



£ 3 9A 1 A 1 X Al l ■ A Im [X B K + iD B ) ^ + s^X^T ■ A Im (X B i^ + iD B ) X A \ (297) 

^1,2 correspond to first generation of leptons and quarks, X = X b 8b = — 12ma' k\^Qb9t b corre- 
sponds to the killing isometry vector and D term generated is given by: 

B _ iira 1 'k\h 7 Q b v b 



V 



(298) 



Utilizing g AlA \ ~ ^W&nai ~ V5, 5^ 2 ^ 2 ~ 0(l)g a2 a 2 ~ V s as given in[124|as well as v B ~ Vs and 
Qb ~ V3(27ra') 2 /, the physical fermion-fermion- Z Boson vertex is estimated as: 



C^ z ~ Max 



V-4 / 



/V"36. (299) 
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Quoting the analytical expressions as given in |13| : 
• Higgs-boson (h, H) exchange: 



a ff 
If 

Expanding the summation 
TiKH) _ 3 



7/ 



64 7T 



0, 
3 

47T 



Cff r C x ° x ° r 



4 m v o 2 — m 2 + i T r m r 



4m 2 o -m^ + irtmj, 4 m 2 q - m\ + i r# m H 



(300) 
(301) 



(302) 



Again, utilizing the value of masses as given in above cases and f h ~ V s , C^xi' 1 ^ (7x3X3^ y ~ 
from above, after simplifying, we have: 



V 5 .V 36 V 5 .V 36 



4.V 3m 2 z 



+ 



// 64vr 



V 36m 2 ; 



Ksj a -mJ) 



V" 



64^ V y-fm 2 



V"3M 2 



pi. 



~ 0{W)-^GeV- z for V ~ 10 
Z— boson exchange: 



(303) 



l ff 



1 



c //z c xx^ 



4m x o 2 — m^ 2 + iTz mz 



2 to 2 (mz 2 - 4m x o 2 ) 2 
m z 4 



(304) 



Utilizing the value of mass as given above, and C^ z ~ (V sa)^C x ^ x ^ z ~ V 3 e from above, after 
simplifying, we have 



;(M0 

77 



64 vr 



_23 _35 _23 _35 

V 30 .V 36 V 30 .V 36 

V 36m 2 z 



4.V sm 2 z 



2 m 2 (m^ 2 — 4m x o 2 ) 2 



V36 \ mj(4m x ») 



64 vr V V"!m 2 



pi, 



O(10)- 43 Ge1/" z for V ~ 10°; 



(305) 
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W) = — 



ff 



2vr 



Cff z C^xlz 



4 m x o 2 — m z 2 + iTz ~m z 



m z 2 ((4 m x o 2 - m z 2 ) 2 + (T z m z ) 2 ) 



2vr 



1 



2 \Cff z \ 2 [m z 2 (m x o 2 - m 2 ) (m z 2 - Am^f 

+r|[m x o 2 m z A + m 2 (24m x o 4 -6m^ 2 m x o 2 -m^ 4 )]} 
+m z 2 \C^ Z \ 2 {(2m x ,, 2 + m 2 ) [(4m x o 2 - m z 2 ) 2 + m z 2 r|]} 
C ff z C^3 Z 



4 m x o 2 — m z 2 + iT z m z 



m z 2 (16 "i< x o 4 



2.16 V^tI m z 2 m x o 6 +V"S m z 2 (16m x o 6 ) 



2vr 



4 m v o 2 — m^ 2 + i T z m z 

*3 



1 



~ V 15 .V 3ffl 

2 7T 



pi 



35 \ 2 
V 30. V 36 \ 

4.V~fM| / 



~ O(10)" 39 GeT/- 2 for V ~ 10 5 . 
sfermion (f a ) exchange: 



r(/) 

77 



1 ^{m f Cl + m x oDl){m f C b + + m x aD\ 



32 7T 



a,6 



Ar At 

/a Ji) 



(306) 



(307) 



a is the index for sfermion mass eigenstates so that a = 1, ..6 for squark and charged sleptons and 
a = 1, ..3 corresponds for sneutrino, and: 



A a 1 2 _i_ I a 6t |2 



(308) 
(309) 



D% = A} L (A a fR y±(A} L yA} R ; 

AJ L corresponds to the Neutralino-fermion-sfermion interactions mediated by L- handed squarks / sleptons 
and AJ R corresponds to the Neutralino-fermion-fermion interactions mediated by R-handed squarks /sleptons. 

Using results from section 3: C x ^ 1l1l = fV~^, C x ° UlUl ~ fV~s. With exactly similar procedure, 
we find: C*§ w « = fV~T5 ,C*% UlU ~ r ~ /V~i. Utilizing above 



Max 



(JX3 U LU~L 



(JXi u LUR 



J2 D ± = Max r<7^^((7 x 3^ r «)*,(7^^^((7 x 3^^)^ ± c .c. ~ f 2 V 



(310) 



and 



A — 2 2 2 2 

At =m f — m v o — m~ ~ — m v o . 

/a J A3 / a A3 



Considering only first two generations quarks/sleptons and assuming the universality in scalar masses 
for first two generations quarks/sleptons, equation (|307p reduces to the following simplified form: 



Af) 
77 



l^m// 2 V 1 + m x o f 2 V w){rriff 2 V 1 + m x o f 2 V io) 



7T 



4 

m o 



~ 10" 52 GeU- 2 ,for m x o ~ V - sM p and V ~ 10 5 . 



(311) 
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Now, utilizing the numerical estimates of coupling summed over first two generation of squarks/sleptons 
given above, we will evaluate the value of 6V/. Using (and quoting verbatim) the form of expression 
from 1 131: 



ff 



E 



64 7T ^ A 3 ~ A 3 ~ 

a ,b fa fb 



cic b 



2 It 2 2\ a 2 

m f m x o (m x o - m f ) A? 



fa 



-4m v o 2 (m 4 f + m^o 2 m 2 f — 2 m^o 4 ) Ar Ar 

*3 V J A3 J X3 > fa fb 



+D1 D\ 



+4mV (mi + 2my)A% A 7 + 4 (mV - mi) A z ~ A 7 

A3 •> X[i ' fa fb V A3 J' fa f b _ 

8m v o 4 (m^o 2 — m 2 f ) A% + 4m. v o 2 (m v o 4 + m v o 2 m 2 f — 2mj)Ar Ar 

A3 v X3 ' fa A3 ^ X3 A3 J 1' fa fb 



+4m v o 2 (5m v o 2 - 2m 2 f ) A| Ar - 3 (m 2 , - 3 m v o 2 ) A 2 ~ A| 

A3 V X 3 J' fa fb V J X 3 ' f a fb 



+ctc b _ 



8m v o 2 (m v o 2 - m 2 ,) 2 Ar Ar + 8m v o 2 (m v o 2 - m^)A| Ar 

A3 v A3 J ' fa fb A3 v A3 J ' f a fb 



+2(m 2 f + 2m x a 2 )A 2 ~A 2 ~ 



A3 > fa fb] 

+C% D b + 2m f m x o [8m x o 2 (m x o 2 - m 2 ) A| + 12m x „ 2 (m x o 2 - m 2 ) A Jq A^ 
+3A| A| -2(m 2 f -4m Y o 2 )A f A 2 ~ - 6 (m 2 , - 2 m Y o 2 ) A| Ar 

fa fb 1 X 3 ' fa f b V J X 3 ' f a f b 



A3 /a / b 



for m x o ~ V 3M p , Aj ~ — m x o 2 , after solving, one gets: 



1$ = l x lX^~0(W)- 48 GeV- 2 . 



ff 



7r 



(312) 



Higgs (h, H) sfermion (/ a ) interference term: 



a 



{h,H-f) 
ff 

ff 



0, 



. 7T 



E* 



E 



°5 °5 



rr 4 m v o 2 — ml + i F r m r 

r=h,H X3 r 



(m x o 2 -mj) 



A 2 ~ 

fa 



[c a + 2m f m xl + Dl (2m x o 2 + 3A /a 



(313) 



(314) 



expanding summation 



;7r 



Cff h c x ° x °3 h 



+ 



c hhH QxlxlH 



Am\-mi+iY h m h Am 2 -mi\ + iT H m H 
A3 n x 3 n 



(m x o 2 - mj) 
A 2 

fa 



C\2m s m x o + Dl (2 m x o 2 + 3 AjJ 

23 35 23 35 

V 36 .V 36 m p l V 36 .V 36 m p i 



TT 



4.V 3 ffl 2 



■> 1— — 

V 36 771 



13 ~o 1 ~ 13 

V^to/ 2 ~ -f 2 y-w 

TT 



2 



O(10)- 42 Gei7~ z for V ~ 10 s . 



(315) 
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Z sfermion (f a ) interference term: 



a 



(Z-f) 
If 



4 7T ^ 



C ffz c^xlz 



4 m v o 2 — m^ 2 + i Y z mz 

*3 



m/K 2 -4m x » 2 )K^ + m x »D;) 



utilizing the value of mass mz = 90GeU, m^o ~ 
and equation no 13101 , after simplifying, we have 



AZ-f) 
If 



1 

47T 



23 ~ H 
V"36 /V~18 

4m v o 2 

*3 



"V (m-^ 2 — 4m x o 2 ) (m/ + m x o Z?" ) 



mz' 



A 7. 



o(io)- 42 G e y- 2 



(317) 



(4m x o 2 -mz 2 + irz mz) 2 | mz 2 A~ 



+ 



C " Z C - { 2 ™z 2 p z A fS 2m x% Ki 2 + A fJ + m ) (" 2 -xg 2 + A /J] } 

c ffZ^ C a 2m 2 m ^2 ^ m ^2 _ m 2^ ^2 _ 4m ^j p ^ m x o 2 [m 2 f m z 2 + 2m x o 2 (m z 2 - 6m})} P z 
+ 2m z {-m z (m x o 2 - mj) (mz 2 - 4m x o 2 ) + i T z [mz 2 "\o 2 - mj(mz 2 + 3m x o 2 )]} A~ 

+mj m x o D" 4m x o 2 (m x o 2 — m 2 ) (mz 2 — 4m x o 2 ) Pz + 2 [6mz 2 m x 9 2 — 16 m x « 4 — to 2 (3 mz 2 — 4m x o 2 )] Pz Aj 
-3 [(m z 2 - 4m x o 2 ) 2 - im z T z (m z 2 - 8m x o 2 )] A?J | ; 

(318) 

where Pz = 4m x o 2 — mz 2 +i Tz mz- Again using the numerical values of masses and relevant couplings 
and assuming T zmz « m x o 2 > the above expression reduces to 



Xiz-f) 
If 



O(W)- 50 GeV-' 2 . 



Utilizing the results from equation no (13031) . rf306j) . rf3T2j) . (l3T5l) . rf3T9j) and ([305jl . (I3T7) : 



r(ft,H) r(z) r(/) 
7/ ~ °/7 "// +b /7 



//;;*'•• + />'".' t />•/; + +bff H ~ f) + ~ O(10)- 42 GeF- 2 ; 



(319) 



(320) 



» 7/ = ag^ + 3$ + 3$ + +4^-^ + 4f ~ O(10)- 42 GeU- 



Relative velocity w/j/a is defined as 



j _ ( m /l + m / 2 )' 

4m 2 



1/2 



i _ (m A -m /2 ) ; 
4m 2 



1/2 



(321) 



(322) 
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For /i,/ 2 = hh,Vhh = 1, /i,/2 = ZZ,v Z z = 1, /i,/2 = ff,v f f = 1. Having estimated the partial 
wave coefficients for each possible annihilation processes, summing up their contribution as according 
to (p33|) : 

a = a hh + a Z z + a ff - = O(10)- 29 GeV- 2 

b = b hh + b zz + b ff = O(W)- 10 GeV- 2 (323) 



f x f f x f x 2 f 

J(xf)= / dx(avyi0\)(x) = / dx(a + bxf)=axf + b— (324) 
•/ii ./ii 2 

where x = T/m x . The value of Xf can be calculated by solving iteratively the equation 



where g* represents the effective number of degrees of freedom at freeze-out {yfgZ — 9). Solving above 
equation, Xf = Tf/m x a comes out to be around 1/33 

The analytical expression of relic abundance is given as |35| 

n x h 2 = - 1 (326) 

where /i = 1.2x itfGeV For J(x f ) ~ a(x/) + 6^ ~ Kr^^GeV -2 and x/ = ^, ^ = 9, 

n h 2 2 ' l 33 ^ 2 ^ = 168. (327) 

x 3 1.44 x 10 10 • 9 • 10- 10 v ' 

For m± ~ 10 8 GeV and m v o ~ Vsm3 ~ 10 11 G ! eV A , relic abundance of gravitino is given as : 

2 A3 2 

ms 2 

ft =,/i 2 = O oh 2 x — = o o x V"3 ~ 0.16, (328) 

* 3 m^o * 3 



X3 

clearly a very desirable value! 



5.2 Slepton relic density calculations 

For the case of slepton NLSPs, the dominant annihilation channel possible in our set up are: £ a l^ — > ZZ, 

Ul zhJJt -> w», 4^ 77, 4i* -> 7 /i, 4^ h. 

The analytical expressions for u>(s) are given in [TJ]. Once again, the approach is to first calculate 
required vertices in the context of M = 1 gauged supergravity and then utilize the same to calculate 
partial wave coefficients. 

(a)£J* -> ZZ 

Slepton-slepton- Higgs vertex 

Expanding effective supergravity potential V = ex.p K G TsTs \Dx s W\ 2 in the fluctuations around 

1 2 

— >• Zi + V^Mp , .Ai — > Ai + V~vM p , contribution of term quadratic in Ai as well as Z{ is of the 

-89 

order V as < -Zj >, which after giving VEV to one of the Zi, will be given as: 
(*Mh. \ 

^(K zai ) 2 (k AiAi f 



V 



X!) 
:5(j 



< Z % > Z l A A 



i 4*1 



0(V 



15 



(329) 
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Figure 26: Feynman diagrams for l a l^ — > ZZ via s-channel Higgs exchange and t-channel £ c exchange. 



Similarly 



! {K Zi2i f{K AlAl f 
Slepton-slepton-Z Boson- Z Boson vertex 



1 A*l 



0(V 



18 



(330) 



In the context of supergravity action, contribution of required vertex will be given by: 
d Ax d Ax G TB f B X TB X rfB A^A V . As given in appendix B 

dAM G T B f B ~ 0{l)d- ai d ai G TBfB ~ V-f^i, 



(331) 



X = X B d B = -Vlmoi 'K(fi 7 Q B dT B ~ V~s. 

Incorporating values from above, the physical Slepton-slepton-Z Boson- Z Boson vertex is propor- 
tional to 



c U t zz_. V */ 2 V 



V^amJ 2 °( 10 ) 4 

Slepton- slepton- Z Boson vertex 



(332) 



The gauge kinetic term for slepton-slepton- Z Boson vertex, relevant to the second Feynman graph 
in Fig. 26 will be given by ^■G t - bTb Vu^bV^b- This implies that the following term generates the 
required squark-squark-gauge boson vertex: 



k 2 V 2 



7 

V9 



K | M7 (2™') 2 C li ~V 



IT 



K 



Ax A, 



~ 0(l(T 4 )/vi ~ /, for V ~ 10 5 (333) 
Quoting directly the analytical expressions given in [14] . we obtain a numerical estimate of w: 



w £ a e;-^zz 



w 



(h,H,P) 
ZZ 



(h,H,P-£) 



(334) 



Higgs (h,H) exchange (+ Point interaction): 



(h,H,P) 
ZZ 



eft 



laZZ 



s 2 — 4m 2 z s + 12m z 
8m| 



~f 2 V--im 2 pl f 2 V-tm 2 pl 



+ 



s — m 2 + iT r m r s — m H + iTiimH 



s 2 — Am 2 z s + 12m z 



8m| 



(335) 
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slepton (£ c ) exchange: 



w 



(!) 
zz 



— — (jqe c z c e*e a z c qe d z* c ey a 
mi ~^ 



z* 



z c,d=l 



Ti - 2(mf + mi + 2m z )T 3 + [mi + mi + 4m~ rrd + 2m z (m~ + mi + 3m z )]T 2 

- 2K+ 4)( m k-^)-ik + k-Hk + ^ ' 

+(m? -m|) 2 (mf - m|) 2 T - ^ + [s(ml + mi -2m 2 z )-2{mi - m 2 z )(mj - m 2 z )]y 2 



\s 2 (mi — m 2 z )(mi — mi z ) + s{—m~ m~ (mi + m~ ) + 3m z (m~ + mi ) — 3m z (m~ + m~ ) + 2m z } 

,2 „,2\2/„,2 _2 \2l 



Considering only first-generation squarks, 

2 



(336) 



(337) 



c,d=l 



(|336l) yields: 



ZZ 



VTLr 



T 4 - Am} T 3 + 6mi T 2 - 4mf Ti + mf T - J 4 + [2smf - 2mi ]3> 2 



[s 2 mi -2sm& +mi~]y Q 

; auxiliary functions define< 
QQGeV , on solving and simplifying the same, we get 



(338) 



71, 7~3, etc are auxiliary functions defined in the appendix of [T5j . Given that mi ~ V 3 M P , m z 



42 / 21 \ 21 8 

■vt - 4v 5 + (y- +s\v~ + V5 



4vV*-4v^vt + (-4vt + sV f + 1) log - 2V ;-v s - 4V }y- 4V ; 



4 / 17 13 \ / 21 

4Vs (-4VT +sV^ + 1) \/s-4Vt 



r 2 




4VV*-4VT (V* + sVf) Vf + (-8Vf + 2sVf) log gz^V^EbZfzE| 



4 / 21 / 21 8 

4V5 \/s - 4V~ I 16V 5 - 4sV~ - 4Vs 



4 / 4 / 21 



2 6 (4V¥ - aV * ) log - 2V !- V r 4V !V' 4V 7 V 10 + Js - Wljs - Avf (-2V 11 - 6sVf + s 2 vf) V§ 



4 / 4 / 21 

s-2VH+Vs-4VHVs-4V^" 



s - 4V5 y s - 4V^ (l6V 5 - 4sV^ - 4V§ 
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4 log 



4 / 4~ / 21 



r 4 = 



4 / 3" 

s-2V5 + ys-4V5 



21 
4VTT 



4vt 



4VT 



1 

+ 3 



3vf 



-4VT + sV 13 / 5 + 1 



+ 18V 42 / 5 - 7sV 17 / 5 V 4 / 5 + s 2 



2 log 



4 / T I 21 



y = 



4 / r / sr 

s-2V5 + Vs-4V5 y s -4V3- ^ 



^/s-4vt (s - 2vf 1 \/s - 4VT 



2(s-Vl (vt +2)) log 



3^2 = 



4 / I" / 5T 
s-2V* +yj a-iV*/ 5 ^ s-AV^ 



4 / 4\ / 21 

s-4V5 U-2V5J ys-4V~ 



3^4 = g (-s 2 + 2vt (sVtt + 2) s - 2 (6vf + 8V 5 + 3V§ )) . 
• Higgs (h, H) (+ Point) - slepton (£ c ) interference: 



w 



(h,H,P- 

zz 



2mi - , 

z c=l 



r ZZr r i* h t a r 

^ n 1 - C e b eaZZ I C e b icZ C e * e ' Z 



,r=h,H 



s — m~ + iT r m r 



\s 2 + s(m~ + m~ - 2m~ - 4m|) + 2m|(m~ + m~ - 2m~ + 2m|) 

£a -t-Jj <-c ' *-6 *-c ^ 



-2[s(mj -mj -m 2 z )(mj - m 2 ~ - m|) 



,2 r^,4 , _4 



4 4 
,2 ^,2 1 ^,2 /^,2 



+2m z {m~ + m~ - m~ m~ + m~ (m~ - m~ - m~ - 2m z ) 
-m|(m| + m| - m|)}]^] . (339) 

Using the universality in squark masses as obtained in Appendix B, the simplification leads to: 



where 



wf z H ' P - J) ee * Re 
zz 2m 4 



E 



7 2 v 3 1 p 



s — mj. + iT r m r 



z \r=h,H 

(s 2 - 4sm 2 z + 4mf ) - 2(sm| - 8m z m~ )T 



(340) 



log 



T = 



4 / 4 / 21 

S-2VT, - y s-4VH y s— 4V1T 
4 / 4~ / 21 



Y^s - 4V5 \/s - 4V^ 
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Figure 27: Feynman diagrams for l a l* h — > Zh via s-channel Higgs exchange and t-channel £ c exchange. 



Summing up the contribution of w for all s, t and u-channels as according to equation ()334[) : 



~ 122 122 

WZZ\ s={i ml ) ~ V 45 ,W ZZ \ s = (4m 2 ) ~ V 45 . 



The "reduced" coefficients azz and bzz will be given by El 

1 



azz 



32vrm| 



-^22 ■ (3 122 

V^§- ~ 10- 9 Gel/- 2 ;6zz ~ — -V^ ^ lO^GeF" 2 



647rm| 



(341) 



(342) 



Zh 



(b)4^ 

The process l a t^ —> Zh involves the s-channel Z-boson exchange, and the t— and u-channel slepton 
t a , a = 1, 2) exchange: 



~{z) . ~{t) . ~{z-t) 
w Ut^zh = w zh + w zh + w zh ■ 



(343) 



Z exchange: 



w 



(Z) = 1 

z/ ' 12m| 



QZZhQili a Z 



s — m 2 7 + iTzmz 



s 2 ^3(m~ - m~J 2 + m z \ - 2s^3(m~ - m~ ) 2 (m 2 t + 2m|) + m%(rn~ + m~ + — 5m|) j 
+m%\ (m\ — m 2 z ) 2 + 4(m~ + m~ )(m? — 5ml) > + (m~ — m~ ) 2 (3mf l + Qm 2 z m\ + 19m z ) 

^ 2 f /" 2 2 \2/o 4 2 2, 4 \ , 2/2 2n2/2, 2n1,^4/2 2\2/2 2 \2 

--m z j(m~ - m~ ) (3m ft - 2m, z m h + m z ) + m z (m h - m z ) (m~ + m~ ) j + -^m z {m h - m z ) (m~ - rrijj 



(344) 



For C zz/l ~ f 2 V-3,C e *b £aZ ~ / and assuming universality in slepton masses; after simplifica- 
tion, above expression reduces to: 

2 



(2) 



1 



12m? 



s — m z + iTzmz 



24 < 2 4 , o 4 2 2 ^ 4 2 4 

s m z — 'ism- m z + om z mj m h m z m~ m h 



.(345) 



3 Here, we are considering: / ~ 10 
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slepton (£ c ) exchange: 



w 



(t) 

Zh 



111% 



y c ete a z c qe c h c iy a z* c i t e dh * ^ _ 2(m 2 + m 2 )7 * + {m 2 _ m | )2 ^ 



c,d=l 



+4- C e ^ h C e ^ z C e y« h *C e ^ z * T 2 U - 2(mf + m|)71 u + (mf - m 2 z ) 2 l% 



z c,d=l 



+ —R e c e*J a z c qe c h c ey a h* c e;e d z* 

z c,d=l 

2 y 2 + i ( . _ ml + m |) K _ m | )3 ,, + { s(m | + m | _ 2m |) _ K + m |) _ 2m | m j 
+ (3m| - m 2 h )(mj a + m|) + -(m^ - m|)(m| - m|) 2 - — (m^ - m|) 2 (m| - mj b ) 2 \\y 



(346) 



Assuming universality in slepton masses of first two generations 

2 2 

c itIah c e* b I c z c I d I a h* c e* b I d z* ^ y^ c l *- lah C^ lcZ C^ lah *C l ^ dZ * ~ ^ c W a z c IiI c h c t%£ a h* c It& d z* _ p v ~% M 2_ 



1 c,d=l 

Therefore, (|346p reduces to: 

?2->)-2i 2 



c,d=l 



'pi 



rrir 



m 2 ~ (7? + 77) + mi (7? + % u ) - 2y 2 + (2smsli 2 - 2mf )y c 



(347) 



Z ~ slepton (£ c ) interference: 



it- 



[Z-t) 
Zh 



1 2 / c zzh C titaZ y 

^" R(3 / 9 ^ 

m Z \ s ~ z zrnz I 



C £*e a z c t b 



+2m 2 z (m~ —m~ + m 2 z ) - -m 2 z (m\ - m 2 z )(m~ —m~)+ s(m~ — m~ +m 2 z )(m~ - m~ - m 2 z ) 



ia 4 



+m z {-2m~ + m~ (3m~ + m~ + 3m|) - 3m 7 

+3m~ m~ — 2m~ + m|(4m~ + m~ ) — m z } — m^im- — m~ + m 2 z ){m~ — m~ + m^) T l \ 
+C^ a/l C^ cZ { (s - m 2 h )(m~ -m~ )+ 2m 2 z (m~ - m~ + m|) + -m|(m 2 - m 2 z )(m~ - m~ ) 

t *-a 4 4 g "a 4 

+ — s(m~ — m~ — m 2 z ){m~ — m~ — m 2 z ) + m z {—2m~ + m| (3m| + m~ + 3m 2 z ) — 3m~ + 3m| m| 



-2m~ +m 2 z (Am~ + m~ ) - m|} + m 2 h {rrij — m~ - m z )(m~ - m~ +m|) J 7 " > 
For C ZZh ~ / 2 V _ 1,C F *^ Z ~ f,C*& h ~ V"ft, after simplification, (p48]l reduces to 



(348) 



(Z-i) 
Zh 



z c=l 



s — m 



| + iT z m z J 



Ql c laZ Qt* b t c h 



4m| + (-2sm| + m?; m|)(^ + J^) 
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Figure 28: Feynman diagrams for £ a £t — > 77 via point interaction and t-channel £* exchange. 



where 



3 log 



3s _ 8y f_j3 f± 2(-4 ± ytI)F5 



/ 3s - 2 (4+ VTE ) V 5 V s - 4V ~F 



3s _ 8yK j3 f±£ (_4 ± yT5)v5 



^3s-2(4+Vl5)V^ V s-iV ^ 



3s+2(-4+v / 15)Vs' 



3s - 2 (4 + ^15) v s y s - 4y- 



Summing up the contribution of w for all s, t and u-channels as according to equation (|348p 

^Z/i| s =(4m2 ) ~ Vis, ^ Zh | s=(4m 2 ) ~ Vis. 



(349) 



The "reduced" coefficients azh and bzh will be given by 

1 



azh 



32vrm^ 



■ViK ~ lO-^GeV-^jftz,, ~ 



647rm~ 



(350) 



(c)« -> 77 

The four-point contact interaction needs to be included along with the t- and u-channel slepton 
,, a = 1, 2) exchange: 



«j;-m 77 ^ W 77 ^ W ll 



Contact interaction: 



slepton (£ a ) exchange: 



77 



(C £ ^)% b 4(T 2 + 2mf Ti + mj To) - (s - Am} ) 2 y 



/ 4 



4(T 2 + 2mf Ti + mf %) - (s - Am} ) 2 y 



where 



1 _ Vv^+y.-4y2i/5 y 1 

F 21/5 (V21/5 + s ) _ y42/5 ! 7l ~ _ 4y21/5 + s 



(351) 



(352) 



*9 



•7(P4 



Figure 29: Feynman diagrams for 



hrj via t-channel £* exchange. 



2 log 



^-y/s-AV 21 / 5 



2] ( ^-y/s-W^ \ y2l/5 

/2 = / = 1 r i; J4) — 

vW* - 4F 21 /5 a s 3 /Vs-4F 21 /5 



v ^+^/ s _4y 21 /5 



Contact - slepton (£ c ) interference: 



w 



(P-I) = 2C l ^{C U ^) 2 eH ab [-4 + (s - 8m 2 ~ )?] ~ (/ 4 V- 5 )e 4 <U [-4 + (a - 8m? ) J"" 



and 



log 



V3V S - 4 V 21/5 - 



" 2 2 



% / s _ 4 y2i/5 v ^ 



.F = — 

vWs - 4y 21 /s 

Summing up the contribution of w for all channels: (|35ip 



w 



the "reduced" coefficients a 77 and 6 77 will be given by 



fl77 ~ 32vrm 2 



V ~ W- 18 GeV- 2 ;b yy ~ 



77 



647rm~ 



V ~ W- 18 GeV~ 2 . 



(c)ej* h -> 7h 

• slepton (£ c ) exchange: 
With 



w 



7/1 ■ 



~(0 



2 r 



+(s + mj[-2mf -2m 2 r )y 
„ (/ ~ v -34/i5 )2 r (7? + 2 7^) + (7r + m | 75.) + (s + m 2 _ 4mf 
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where 



7^ = V 



7? = v 



42 / 01 , r 4 \ 21 / 21 8 4\ n1 ' 

sVt + s ( V 21 / 5 - Vs + sj VT + (-2sV~ +Vr> - sV^j V 21 / 5 
7-u = v f 

42 / 21 4 \ 21 / 21 8 4 \ 21 ' 

sV~ + s [V- - V~> + sj VT + f -2sVT + V5 - sVs J V~ 

(s-v* -\[s-vi \J 1-— \/s-4vir \ / 4\ 2 / rr / — / 
. I— 'I +5V /l-3 ? \/ 5 -4VT V /,-V 



1 "L ; 

5 



77 = v 



-tHWH '" w +n /i-^-w<- 



4 

V5 



(*-vt) 5/2 yrw^v~ 



2 log 

3^o = V 



(*-vi) 3/2 yr^r^ 



For s = 4m 2 - 



W~/h\ s =(4m 2 7 ) ~ V 15 ^Zzls=(4m| ) ~ V « 

The "reduced" coefficients a 7 ^ and will be given by 

a lh ^V-tI ~ Kr 55 GeV -2 ; 6 7fe ~ ^V _ S ~ KT^GeV -2 . 

32irmf 647rm~ 

(d)£J* -> hh 
Now, 

W l a ll^hh ~ W hh + w hh ^ W hh 

where 

• Higgs (h,H) exchange (+ Point interaction): 
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Hps) 



(h,H) 
> < 



h(p 4 ) 



■ Hp±) 



Figure 30: Feynman diagrams for i a i^ hh via s-channel Higgs exchange and t-channel I* exchange. 



w 



(h,H,P) 
hh 



fihHr fil%t a r 

^2 —4 - c e b £ahh 



r=h,H 

.59 2 



s — m$ + iT r m r 



(361) 



For m h = 125GeV,m H ~ V^ms, m x o ~ V^mg ~ V"3M P . and ~ (V" 1 ), C £ai » h ~ 



V is , (jtatbhh ^ y lg f rom a bove, after simplifying, we have 



V 15 m£ 



+ 



49 

V~15 



s - + iT h rrih s - m 2 H + iT H rriH 



41 

V~18 



(362) 



slepton (£ c ) exchange: 



c,d=l 



c,d=l 



-2Re ^ c e * eah C e b ecH C e *d e * H *C e b edh *y . 



(363) 



c,d=l 



Strictly speaking we are considering first generation of squarks which get identified with same modulus 
Al. Therefore, 



£ cniahQ^ioHfjr^fjtitiH* _ (v-it) 4 ~ V~ 9 M 4 . 



c,d=l 



Hence 



42 = V- 9 (T * + V-^o)M p 4 . 



(364) 



The contribution of 7^,7^\3^o for — > hh is of the same order as t a t* h — > ZZ. 
• Higgs (h, H) (+ Point) - slepton (i c ) interference: 



c=l \r=h,H 



s — m~ + ir r m r 



c e*j a h c e* b e c H p + c e* c e a H c e* b e c h -pu 



, (365) 
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4Cpi)" 



1(P3) 



i(P4) 



Figure 31: Feynman diagram for l a l\ — > hh via t-channel i* c exchange. 



which expanding and incorporating in the values for the couplings, 



w 



(h,H,P- 
hh 



= 2V"« 



+ 



V 15 



log 



where P = = 



s — mj, + iT r m r s — m 2 + iT r m r 

-4V 2 -2V 2 / 5 +s--y/s-4V 2 /-V5-12V 2 
-4V 2 -2V 2 / 5 +s+v / s~4V 2 /Vs-12V 2 



41 

- V~T8 



y/ s - 4V 2 /Vs - 12V 2 
Summing up the contribution of w for all channels: 

Whh\s=(4ml ) ~ V"^>^| s=(4m 2 j ~ V~^>, 

and the "reduced" coefficients ahh and bhh will be given by 
dhh 



32-irml 



' V - n ~ 10- 25 GeV~ 2 ; b hh ~ - 3 „ ~ l(r 25 GeF- 2 



6ATTm~ 



(366) 



(367) 



(368) 



(e) £ a 4 

This process proceeds only t- and u-channel neutralino (x°, i = 1,2,3) exchange {wj 



w 



E 



neutralino (x^ 1 ) exchange: 



-(mj -mj)(mj - m 2 e )T 



LLLL" L xV"-x\ 

me.m x o 



r-.abij r-r~ i j-^abij r-r~ 

D LLRR m x° m x J 7 + LRRL'l 
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LRLR 
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DTAdTi (mj ml)%} + D^ LL {71 (m\ m$)T } 



1 3 



2Da LLLL m x1 m X°S S ~ 2m i)yo + ^ D ThRR m x\ m X° m iyv ~ 2D L°RLR m e( m j a + m \ - 2m^)^0 



abij 2 ( 2 



+2L> 



abij 
LRRL 



-y 2 - (mj mj - mj)y 



+Dl% L {yi + {s- ml + ml - Amj)y } 



+D R b H L {yi -(s-mf+ ml 4mj)y } 



la lb 

,2 , 2 



m<m y o 



— mim x o 

-.abij 



Dl% R {yi + (s + mj ml - 4mj)y } 



DIIllWi -(' + mj a m\ ~ Amj)y } 



(369) 
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where 
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(370) 

Given that m x o ~ V^M p ,m x o = m x o ~ V _i M p , Cg* a " ~ fV~5,Cg" a " ~ /V~it, by expanding 
summation, the above expression reduces to 



sWs - 4Vir 8sV z - 160V 5 - 2 -160V sV~ + 9s z + s8V~ log 



21 



^ " 



For s = 4m~ 



„ /„ 29 / 21 

2s 3 / 2 V M/s -4V^ 



~ i —12—/ | _ !£ 

w«| s =(4m2 ) ~ V s ,u; tf | s=(4m 2 j ~ V s . (371) 



The "reduced" coefficients and 6^ will be given by 



-^-g-V - ^ ~ 10~ 41 GeV- 2 ;6 a ^V - ^ ~ 10~ 41 GeV^. (372) 

327rm~ 647rm~ 

Summing up the contribution of partial wave coefficients for each possible annihilation processes: 

a i a q^hf2 = ® zz + + Qhrf + a-y-y + ^hh + ail = 0(W)~ 9 GeV~ 2 

hji^hh = bzz + b h z + b h ^ + b^ + b hh + b u = 0(W)- 9 GeV- 2 , (373) 

f X f f X f Xl 

J{ x f)= \ dx(avM0\)(x) = / dx(a + bx f) = ax f + b-J- , 
Jo Jo 2 



and 



where x = T/mj . 

Using the analytical expression of relic abundance as given in (|326j) 



% = 2 (374) 
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where fi = 1.2 X 10 5 GeV. For J (a;/) ~ a(x/) + b^- ~ lO^s/GfeV -2 and x/ = | 



33 ' 



33 

"t ~ l-44xl0"».»-10- = °' 25 ' <375) 



for m-± ~ 10 8 GeU and mj ~ V 2 TO3, relic abundance of gravitino will be given as : 

2 *-a 2 

TO3 1 

SI* = n 7 x —2- = 0.25 x V"2 ~ 0.001 for V ~ 10 5 . (376) 

From (|328|) and (I376p . it appears that relic abundance of gravitino turns out to be suppressed in case 
of slepton (NLSP) (co-)annihilations as compared to relic abundance of gravitino in case of neutralino 
(NLSP) aannihilations, for almost similar value of thermal cross-section. The reason, as explained in 
[8], is based on the fact that in the context of partial wave expansion approach, neutralino annihilation, 
because of the Majorana nature of the neutralino, occurs in P-wave while slepton annihilation occurs 
in S-wave. Therefore relic abundance of sleptons gets suppressed as compared to that of neutralino 
by a factor of — . 

6 Results and Discussion 

In most of the supersymmetric models, gravitino appears as most natural candidate of dark matter, 
the stability of which is governed by R-parity conservation that was initially proposed to explain the 
non-observation of proton decays in collider physics. However in split SUSY models, the constraints 
appearing from R-parity conservation are relaxed because of the fact that high scalar masses help to 
avoid the fast proton decay. Therefore it is interesting to see if gravitino appears as a stable dark 
matter candidate by taking into consideration trilinear R-parity violating couplings. 

Large volume Swiss-Cheese type IIB compactifications [39] have been of wide interest, both from 
cosmological as well as phenomenological points of view. In this paper, we have investigated, in detail, 
the possibility of gravitino as a viable dark matter candidate in the framework of "L(arge) V(olume) 
li- split SUSY" scenario. The presented scenario includes four Wilson line moduli on the world volume 
of space-time filling D7-branes wrapped around the "big divisor" and two position moduli of a mobile 
space-time filling Z)3-brane restricted to (nearly) a special Lagrangian sub-manifold. We have shown 
that the M = 1 gauged supergravity fermionic mass term corresponding to fermionic super-partners 
of two (A\ and .A3) of the aforementioned four Wilson line moduli matches very well with the order 
of Dirac mass of the electron, and the M = 1 gauged supergravity fermionic mass term corresponding 
to fermionic super-partners of the remaining two (A2 and A4) of the aforementioned four Wilson line 
moduli, matches very well with the order of Dirac masses of the first generation SM-like quarks, such 
that the fermionic superpartners of A\ and A3 get identified, respectively with ei and en, and the 
fermionic superpartners of A2 and A4 get identified, respectively with the first generation quarks: u/di 
and u/d,R. In doing so we also show that the RG-ffow equations for the "effective" Yukawa couplings 
relevant to the fermionic mass calculations, change by 0(1) in flowing from the string to the EW scale; 
it has been shown in [10J that with some fine tuning, it is possible to RG flow the value of (zj) (identified 
with Higgs VEV) ~ V^M p at string scale to (zi) ~ 2A6MeV at the EW scale. In fact, using the 
RG-flow arguments of [10], one can also show that the Weinberg-type dimension-five Majorana-mass 

generating operator: 0((zi) 2 ) coefficient in , = =| {x^t^i) or in fact - , = ^ =i — (xt 1 z i 



produces the correct first-generation neutrino mass scale of slightly less than leU for (zi) ~ 0(l)Vs6 . 
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We also show that it is possible to generate the massive gauge boson mass scales corresponding to 
the W/Z-bosons, under a similar RG flow. Building up on this set up, we have also evaluated masses 
of all supersymmetric as well as soft SUSY breaking parameters obtainable from bulk F-terms in the 
context of gravity mediation. Using RG-flow solutions of scalar masses, similar to the ones for a single 
Wilson line modulus set up of one can again show that values of same do not change from string 
scale down to EW scale. On evaluation of the soft supersymmetry breaking parameters, we find a 
universality in the squarks/sleptons masses and assuming non-universality in the D3-brane position 
moduli ( to be identified with the neutral components of the Higgs doublets) masses, on diagonalizing 
the Higgs mass matrix we obtain, at the EW scale, one light Higgs of the order 125GeV and one 
heavy Higgs while the Higgsino mass parameter comes out to be heavy. The lightest neutralino is 
predominantly gaugino(Wino/Bino)-type formed by linear a combination of gaugino and Higgsino. 
The mass scales of various SM and superpartners are given in table [TJ After calculating the masses of 
various SM and their superpartners, it appears that gravitino is the Lightest Supersymmetric particle 
(LSP) for Calabi-Yau volume V ~ 10 5 which motivates the query: can we have gravitino DM in gravity 
mediation scenarios ? 



Quark mass 
Lepton mass 


M q ~ 0(5)MeV 
Mi ~ 0{l)MeV 


Gravitino mass 
Gaugino mass 
Neutralino mass 


D-^-l 
7743 ~ V 2 

2 2 

Mq ~ V ^rnz 
2 2 

M v ~ U 3 777 3 

A3 2 


.D3-brane position moduli 

(Higgs) mass 

Wilson line moduli mass 


1 2 
1 

771 7 ~ V 2 m3 
Ai 5 

7 = 1,2,3,4 


A-terms 


A vgr ~ U 1/36 772,3 
2 

{p,q,r} e {Ai,Zi} 


Physical //-terms 
(Higgsino mass) 


~ V 36 771 3 

2 


Physical fiB-terms 


(f lB )z 1 z 2 ~ Vfiml 

2 



Table 1: Mass scales of first generation of SM as well supersymmetric and soft SUSY parameters 

In big bang cosmology, gravitino population depends on two kinds of mechanisms: thermal as 
well as non-thermal production. We assume that reheating temperature is low enough to produce the 
effective relic abundance of gravitino in agreement with experimental observations, therefore almost 
all of the gravitinos are produced by electromagnetic as well as hadronic decays of unstable NLSP. 
The scale of masses of various superpartners suggest that because of an 0(1) difference between the 
masses of the sleptons / squarks and the lightest neutralino and given that in our calculations we have 
not bothered about such 0(1) factors, both can exist as valid NLSP candidate if life-times of the same 
are such that they do not spoil the bounds given by Big Bang Nucleosynthesis (BBN). Based on this 
fact, by considering the contributions of all required couplings in N = 1 gauged supergravity action 
and assuming that similar to the aforementioned effective Yukawa couplings, there is only an 0(1) 
multiplicative change in the three-point interaction vertices in RG-flowing from the string to the EW 
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scales, in section 3.1 we have discussed in detail the leading two body decays of gaugino/neutralino 
B — > ip^Z,^ and three body A° 2 — )■ ip^ t uu,ip fl W + W~ decays and determine that life-time of these 
decays comes out to be too short to affect the prediction of BBN. In 3.1.2, we have (re)calculated 
the tree-level as well as one-loop gluino decays into neutralino as well as Gravitino/Goldstino in four- 
Wilson-line-moduli set up, large life time(s) of which satisfy one of the important phenomenological 
features of "/i- split SUSY" but can not be considered as an appropriate NLSP's because late decays 
of the same can surely elude the constraints coming from BBN. Furthermore, since we are considering 
R-parity violating couplings into account, in addition to NLSP decaying into LSP, there might be 
chances that neutralino(NLSP) directly decays into SM particles via R-parity violating couplings 
and hence affect the relic abundance of gravitino produced by neutralino. To ensure that this does 
not happen, in 3.2, we have studied in detail the three-body decay width of three-body decay of 
neutralino into SM particles (x-j — > ulSrIl), which comes out to be very much less than the two- 
and three-body decay widths of ordinary neutralino decays into LSP. Therefore, we show that R- 
parity violating couplings do not affect relic abundance of gravitino. On a similar ground, we have 
calculated three body decay width of sleptons I — > Gl,l — > I'GV, which similar to neutralino decays 
are short enough to effect bounds on BBN. In order to satisfy the requirement of DM to have life 
time of the order of/greater than the age of the universe, in section 4, we have estimated the mean 
lifetime by calculating two-body and three-body decay widths of Gravitino decaying into SM particles, 
which in case of two-body decays, comes out to be 10 17 seconds and in case of three-body decays, 
comes out to be 10 21 seconds, i.e., greater than the age of the universe. The numerical estimates 
of various N(LSP)'s are provided in tabled! In short, the explicit calculation of life times of both, 
LSP and NLSP, confirms that gravitino can be considered as an appropriate dark matter candidate. 
However this is not the end of the story. The viable dark matter candidate should also have the right 
order of relic abundance in the range provided by WMAP data and other direct as well as indirect 
experiments. With the assumption that the next-lightest supersymmetric particle (NLSP) freezes out 
with its thermal relic density before decaying to the gravitino and then eventually decays to gravitino, 
relic density of gravitino will be completely given in terms of relic density of NLSP. Therefore, in 
section 5, utilizing the partial wave expansion approach, we explicitly calculate thermal annihilation 
cross-section of sleptons and neutralino which ultimately give gravitino relic abundance from sleptons 
to be 0.001 while the same from neutralino comes out to be 0.16, almost in agreement with the value 
of n c h 2 suggested by the WMAP 7 -year CMB anisotropy observation [2]. 

An interesting outcome of our current project is that for trilinear R-parity violating couplings A^- fc 
(to be interpreted as an effective Yukawa couplings in M = 1 gauged supergravity evaluated in 4) 
~ V~3 and X'-j k (to be interpreted as an effective Yukawa couplings in M = 1 gauged supergravity 

evaluated in 4) ~ V~5o, according to the analytic expression given in |37| . the rough estimate of 
proton decay width comes out to be 

2 

r p _ e+7r o ~ m 5 proton - — ~ 10~ 73 GeV, (377) 



mi mi 



thus giving life time of about 10 42 years, which explicitly governs the stability of proton in the presence 
of small R-parity violating couplings discussed in /x-split SUSY scenario. 

To summarize, we conclude that the gravitino qualifies as a potential dark matter candidate in 
Large volume"^ split SUSY" scenario. 



97 



Decay channels 


Average life time 


Gaugino decays: B — >■ i^^Z/j 


10- 30 s 


B ip^uu 


10~ 13 s 


W ibuUU 


10~ 15 s 


w° ->• v u tu + tu~ 


10" 66 s 


Gluino decays: (7 — ^ XrS-iQj 


10 _5 s 


o — >■ y2o 


10 2 s 


o — >■ ib,,q T q ' 


10 4 s 


9 -> ^7*5 


10 n s 


RPV Neutralino decay: X3 ~~ ^ "ude - 




Slepton decays: Z — >■ Z'GV 


10- 28 s 


l/q-tl/qipp 


10 -25.5 S 


Gravitino decays: ^ M — > hu e 


10 i7 s 


ipn ->• 1/7, 


10 21 s 




10 21 s 


% -> LiQjDl 


10 20 s 


% UfD^Dl 


10 18 s 



Table 2: Life time estimates of various N(LSP) decay channels 
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A Geometric Kahler Potential 

The crux of the Donaldson's algorithm is that the sequence 



—didj J In ^ h a/3 s a Sf. 



kir 



on P({zi}), in the k — > oo-limit - which in practice implies k ~ 10 - converges to a unique Calabi-Yau 
metric for the given Kahler class and complex structure; h a p is a balanced metric on the line bundle 
Cp({ Zi }) (k) (with sections s a ) for any k > 1, i.e., 



Wa P = V 2^ 7^ ( \ 7 \ - n a/3> 
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i*(J 3 ) A 

where the weight at point pi, Wi ~ oao w ^ n * ne embedding map i : P({zi}) ^ WCP 4 and the 
number of sections is denoted by N^. The defining hypersurface of the Swiss-Cheese Calabi-Yau in 
the X2 = 1-coordinate patch in WCP 4 [1, 1, 1, 6, 9] is given by: 

1 + zf + zf + zl + Z 2 A - iJZlZ 2 Z 3 Z4 - &j>zl4 = °- ( A1 ) 

In the large volume limit, the above can be satisfied if, e.g., l+z^+Zg 8 ~ — zf, z\ ~ i^ziz^z^za+^z^z^. 
For zip ~ V36, one sees the same are satisfied for zs : a ~ Vs provided ^V 1 * ~ 30. Therefore: 

hlZi ~ 2-2 2 1 hlz 4 ~ 2 _ 2 2 

/l 2 4 z 4 V 3 /i, 2 4 2 4 V 3 

Vis yt 

^1z 2 ~ — — 2"; ^lz 2 ~ — ^ — 2~ 

4 hA^ivl 4 ^lyf 

1,1 1 

Ve + 36 Vis 

hlziZA ~ I"' ~ 2^2 I" 

/i z 4 z 4 V 3 /( Z 4 2 4V 3 

J_4- 1 J_ 



h Zi Z4 ~ 2-2 2! hziZjZ k ~ 2-2 2 

/j Z 4 Z 4V3 /l Z 4 2 4Vs 

V36 + 3 yis + 6 

^2i2~ 2 ~ 2^2 2~) ^2i2,'24 ~ 2~^2 2~ 

14 h z l z lv* h z l z lVz 

1 1 I 1 

V3 Vl8 + 6 

/iz„z„ ~ ^ 5", /i,.?2 ~ 



2 4 2 4 2=2., 2' '"242/ ,J=2,,2 
/l^.^2 ~ — — TT, /l 



2 4 2| , 2=2 2 ! '"242i24 2=2 2 

/l 2 4 2 4V3 Al 2 4 2 4 )/3 

V5 yT8 + l 

h, ~ ~ ~ rv-' n, - o r\_- 

1 lj Z.^Z-;ZuZi o o 1 Jfc^.^.^Z 



ZiZjZ k Z[ 2=2. 2 5 "ZiZjZl 2=2. 2 

/l 2 4 2 4V3 h 4 iV't 

1 i 1 

VI2 + 6 1 

hz;ZAZi~ZA ^ n o o"i 



ZiZjZ k Z 4 2=2., 2! ,[ -2|2| ,^2=2 

/l 2 4 2 4V3 ft 2 4 2 4 

J_ , 1 _1_ , 1 

\?36~ r 2 V 1 * 3 

^2|2iZ4 ~ , -2=2. ,2 ) hziZ^Zi ~ , 9=2. ,2 ' (^ 2 ) 

/l 2 4 2 4 V 3 ft, 2 4 2 4 J/ 3 
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which on being inverted gives: 



( 


h 11 


h lzt 


h U 4 


h lz ? 




f l lz i z 4 


\ 






h ZiZj 


foZiZi 




h ZiZ % 


J^Z{ZjZ4 






h z 4 l 


foZAZi 


^2424 




_2 


■foZAZiZA 








J^ZiZj Zfc 


J^ZiZj Z4 


UZ^Zj Z\~Z\ 


_2 

Jl Z i Z 3 Z 4 


y^ZiZjz^z^ 






h z l x 


2 - 
fl Z 4 Zi 


foz\z4 


2- - 


fl Z 4 Z 4 


2 — - 
fl z 4 z iZ4 




\ 




foZiZAZk 


J^ZiZ4Z4 


LZi Z^ZjZf, 


-2 

JjZiZ^Z^ 


f l Z i Z 4 Z j Z4 


J 


( 


h zlz%y2/3 h z 


2 z 2 y23/36 




V h z 


fzfyll/18 






^z 2 z 2 y 23/36 foz 




h44y 


17/36 ^z 







h zZzlyll/18 

h z l z l 



)l z 4 z 4 

2 -2 
}l z 4 z 4 

h z\zl v \\im 

fl z 4 z 4 



yl7/36 
y7/12 



yi/9 



h Z 4 Z 4 

h z 2 z 2 y 4/9 

h z t z i yv 

\l z \ z \ 



- 2 yA/9 



yll/36 



h z i 

2 -2 
fl z 4 z 4 

h zply5/18 
\l z \ z \ 



yb/9 
yb/12 



]l z 4 z 4 
1 Z 2 Z 2 

\l z 4 z 4 
}l z 4 z 4 

frzply 5/36 



y-11/36 

V 5 / 18 



h z%4yl7/36 ^ 
h z 2 zlyi/9 

h z %- 

iff. 



yll/36 
y5/12 
y5/36 
! y5/18 y 



(A3) 



Using (|A3p . one hence obtains the g Kahler potential ansatz of ([25 



B K T f . the Massive Gauge Boson mass and Soft SUSY Parameters 
Miscellania 

The relevant term in the N = 1 SUGRA action of [12], assuming that for multiple D7-brane stacks, 
the non-abelian killing isometry of the moduli space gets identified with the SM gauge group, the mass 
squared of the VF/Z-bosons are given by: 9t b t b {X Tb ) ■ Now, 



1 3[T B + T B - Cjjcuaj - /X3PVT8 

G T B f B - 1 * ^ '- (Bl) 

>T B + f B - CjjaTaj - fi 3 l 2 VS E 



3 3 



where E = ^Tg + Tg — C 7J -a/aj - W lW - ^Tg + T B - #,rV" ) + E/3 *$(■•■)■ So, to get a 
mass of 90GeU or so at the EW scale, one needs to show that (IB1|) can RG-flow down to the required 
very small value of the square of the same; lets call this small value as GjT^ ~ 10~ 20 . This would 



imply the following relation between Tg = Tg + T B — Cjjaiaj — ^3! Vis and E: 

37g 



whose valid solution is given by: 




3TJ + 97^Gg^ + C(Gg^ 2 ). (B3) 
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The above is equivalent to: 
which can be satisfied, e.g., if 

Now, 



T B {EW) ~ 0.6VI 



(ai)(EW) 



0(1 



■<a/>(M s ). 



Similarly 



Therefore, 



da,>d ar G Ta f n ~ 



^z^TbTb 



Cjja 1 9C u -a J , n TlC 7j a J 



2 OV2 

27^ s ^ 



+ 9- 



s 3 



A* 3 J 2 ^ 9^ 2 ^ , n T* (J, 3 l 2 z~i 



3 ny2 

2TJ S /2j 



+ 9- 



£ 3 



(B4) 
(B5) 

(B6) 
(B7) 



C 



IK 



+ 



3C IJ a J C LR a L , VC lS a J C LR a L 9C lR , 27 C l3 a J C Lk a L y/T5 , 9C /J? T, 



+ 



2rjs 4 rj 4 s27 " B 



2S 2 



+ 



s 3 



+ 



3 

, 2 

/if 'B 



81 T^Cjj^C^ + 27 VnCjja J C LR a L 



E 4 



e 3 



s 3 

(B8) 



Now, for 7 = K = a\, CLai(a L )(Ms) ~ Vs; which because of (|B5j) . implies that at the 1<W scale, 



da R d ai G TB f B 



3 2 

E(£W)~3T^(£VK),T S (-EVK)~0.6V3 



v~». 



(B9) 



Similarly 



dz- k d Zi G Tg f B 



fi 3 l 2 3 (y 3 l 2 ) 2 zh l 9 (ji 3 l 2 ) 2 zh l /i 3 / 2 27 Qu 3 / 2 ) 2 ^ L Tj^f 
■ A ^| + 4 E 2 ^ £ 2 2 s 2 r B 



2T B 2 S 



T 2 

' B 



E 3 



+9 



TJ/^ 2 81 T^ 2 ) 2 ^ , 27 V7h(^ 3 l 2 ) 2 zh l 



+ 



E 3 2 E 4 
Now, for i = k = z\, at the £W scale will be 



E 3 



Tb(£W)~0.6V§ 27gE 



V 



-2 



(BIO) 



(Bll) 



T B (EW)~0.6V3 



We will briefly describe evaluation of various soft supersymmetry breaking parameters in our 
current setup involving four Wilson line moduli. To begin with, in order to evaluate the gaugino masses, 
one needs to evaluate the bulk F-terms which in turn entails evaluating the bulk metric. Writing the 

3 3 

Kahler sector of the Kahler potential as: K ~ — 2ln (<tb + &b — 7-^geom) 1 — (^s 1 + &s — lK rrrigeom )^ 



i S(n m) cos ' ~~ G)9s + m k • (G + G)) , and working near 
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components of the bulk metric: G r , 



the following inverse: G Tl 



sin (ink ■ (G — G)g s + mk ■ (G + G)) = - corresponding to a local minimum - generates the following 

/ , . 37 59 \ 

/ V~36 y-36 \ 

y-§ y~i o o 

o o 0(1) 0(1) 

\ 0(1) 0(1) J 

/ vl vi \ 

13 4 

Vis y 3 o o 

0(1) 0(1) 

\ 0(1) 0(1) J 

F m = e%G mn D n W, one obtains: 



, which therefore produces 



Given that bulk F-terms are defined as: 



implying that the gaugino mass is: 



(B12) 



F m d m T B < 2 



(B13) 



where the gravitino mass, 7713/2 ~ V 2 M p . 

To calculate the mixed double derivatives of the Kaler potential (|25p with respect to (|31|) . we use 
((52)1 and: 



E 2 



^l 2 C jI aj VTBVTsuzPziCf^ VjsVnC ri aj 



nc IJ aj 



£ 2 



s 2 



d ai da,K ~ 



C lR a R C Lj a L T B C lR a R C Lj a L VTbC iJ 

S 2 s : 



(B14) 



we see that: 



<9z fc dz. 



<9aj cfe; 



2 

V-3, 



(B15) 
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implying: 



d <y B [dzidz.K 



da B d aB [d Zi dz.K 



X 2 



4 

V"3. 



9a B Ba s \ d Zi d Z K 



v~W 2 , 

_ 43 9 

dG°da B , s (d Zi 8z.K) 



d G ad & idz,d z K 



nn(ink-(G-G)g B +mk-(G+G))=0 



sin(ink-(G-G)g s +mk-{G+G))=0 



V"3. 



From ()B12p and (|B16p . one obtains: 



F as F* B d as d„ B ln (dzMK) ~ V 



F G ° F* 5 ' 9 G a ^6 Zn ( d^. % A 



sin(ink-(G-G)g s +mk-(G+G))=0 



implying that 



Similarly, utilizing: 



iriZi ~ V 72 m 



3/2- 



oAi dA\ oAi dA\ 3 oA\ dA\ oA\ dA\ 



daj dz 



dzj da 



dz; dzj 



along with: 



yields: 
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\F^\ 2 d as d- as ln [d Al 8 Al K)~V 
\F° B \ 2 d aB d- aB ln {d Al d M K) ~ V"f , 

_ _ q r 

F° a F & 'd G ad Qb ln (d Zi d Zi K 



V 



sin(ink-(G-G)g B +mk-(G+G))=0 

m Al ~ VVm 3/2 . 
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Due to the logarithmic derivatives in the definition of the open moduli masses, in fact, one can show 
that there is a universality in the open moduli masses and that m Al ~ VVm 3 /2- 

We now show that the universality in the trilinear ^-couplings that was seen in the case of the 
D3-brane position moduli and a single Wilson line modulus in [9], is preserved even for the current 
four- Wilson- line moduli setup. 

Using, 

Atjjc = F m (d m K + d m lnY XJIC + d m ln [K xl K j3 K K £)) , (B22) 

and noting: 

35 

F m d m K ~ V"T8, 



F m d m lnY ZiZjZk 



V 36 . 



sin(ink-(G-G)g a +mk-(G+G))=0 

F m d m ln (k? .5 K? K, 



m 1 "- \ lv ZiZi 1 ^ ZjZj lv Z k Z k 



4 



sin(ink-{G-G)g s +mk-{G+G))=0 

F m d m lnK MAi = F m d m lnK MM = F m d m lnK MM = V" 2 , (B23) 
one sees that one obtains a univeral trilinear ^4-coupling: 

37 

M.JK ~ V36m 3 / 2 ~ (iz x z 2 - (B24) 



One can show that 



A* 



^ ZlZ2 = e ~ vim 3/2 . (B25) 

V K z 1 z 1 J ^z 2 Z2 



Further, 



A* 



e -iarg(W)+ 2 



(f lB )z 1 z 2 = i tv" _ K =i?m {dmKfx Zl z 2 + d m [i Zl z 2 - fiz 1 z 2 dmln {K ZxZx K z ^) 

\J "~Z\Z\"~Z 2 Z 2 

19 , 37 9 9 

~ V is+36 m^ /2 ~M| l22 , (B26) 



an observation which was very useful in obtaining a light Higgs of mass \2bGeV for a two Wilson line 
moduli set up of [TT] - a conclusion which is still true using the same as shown below. 

The Higgs masses after soft supersymmetry breaking is given by {m 2 z . + fi 2 z .) l l 2 (where m^s cor- 
respond to mobile D3- Brane position moduli masses (to be identified with soft Higgs scalar mass pa- 
rameter)) and the Higgsino mass is given by pLz v Had the supersymmetry been unbroken, Higgs(sino) 
masses would have had been degenerate with coefficient [iz i - Nevertheless we have defined SUSY 
breaking but we are still justified to use RG flow equation' solutions because fiz t » rnz v However, 
due to lack of universality in moduli masses but universality in trilinear couplings, we need to 
use solution of RG flow equation for moduli masses as given in [38] . 

m 2 Zl (t) = m 2 (l + 6i) + m\ /2 g{t) + ^S p, (B27) 
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where 



So = Tr(Ym 2 ) = m% 2 - m\ x + ^{m] L - 2m\ R + m| - m? + m\ R ) 



(B28) 



in which all the masses are at the string scale and n g is the number of generations, p is defined by 
p = ^[1 — ( a°(if ) )] where di = g 2 /(47r) 2 and 51 is the *7(l)y gauge coupling constant. Further, 

3 

™l 2 (*) = m-o A ^2 + "ii/2 e W + A om 1/2 f(t) + m 2 /i(t) - k(t)A 2 D - -S p 
where A Z2 is given by 



(B29) 



= g 2 + 3 + 4 ^ + 2 ' = ~ ^MtJ 



(B30) 



Here = (Ms)/(47r) 2 where 1^(M S ) is the physical top Yukawa coupling at the string scale which 
following [?] will be set to 0.08, and 



E(t) = (1 + fat) 3 "3 (1 + fat) b 2 (1 + fat) 9 "i 



(B31) 



In equation ()B3ip fa = oti{M s )bi/Ai: {a\ = (5/3)ay), £>, are the one loop beta function coefficients 
defined by (h,b 2 ,b 3 ) = (33/5,1,-3), and F(t) = /„* E(t)dt. 

In the dilute flux approximation, gl(M s ) = gf(M s ) = g 2 (M s ). To ensure E(t) G R, the 5*7(3)- 



valued 1 + fat > 0, which for t = 57 implies that 53 (M s ) < 



3x57 



0(1). Hence, sjj(M s ) = 0.4 is 



what we will be using (See [IT]). From (|B27|) . the previous results in this appendix, one sees that: 



1 19vr 



m Zi (M EW ) ~ mi x (M s ) + (0.39) V^m^ + — x — 



x S , 



(B32) 



and 



m% 2 (M EW ) ~ (0.32) Vtm 2 /2 + (-0.03)n s /} Zl 2 2 Vim3 /2 + (0.96)m^ - (0.01)(n s ) 2 /i^ 2 



197T 

2200 



x S , 
(B33) 

where we used Az i z i z % ~ n s fiz x z 2 (|B24[I . The solution for RG flow equation for /} 2 to one loop order 
is given by 



wherein 
Ci = 

C 2 = 



1 

tan 2 j3 — 1 
tan 2 f3 



mgCi + + ml C 3 + miA C 4 - -Mf + 



1 9 19vr (tan 2 (3 + 1 



3go_-l^ 2 



tan 1 fa) + 



1 



tan 2 /3 - 1 



2200 V^ 2 /3 - 1 



(B34) 



tan 2 /? - 1 



fc(*);C 3 



<Ji - 5 2 ton 2 /3 5_ — (<5 2 + £ 3 + ,5 4 )ton 2 /3 ) ; 

tan 2 /? 



tan 2 /3 - 1 



(o(t) - tan 2 fa{t)) ; C, = - tan2p _^ it) 



(B35) 



and where the functions e(t), f(t), g(t),k(t) are as defined in the appendix of [UJ. In the large tan/3 
(but less than 50)-limit and assuming 81 = 82 = 0, one sees that: 



-2 

a z x z 2 



19vr 



-m 2 - (0.01) (n 5 )^^ + (0.32)V3m 2 /2 - 1/2M^ + (0m)V*n s (, ZlZ2 m 3/2 + —S 

(B36) 



105 



From ()B33|) and (|B36P one therefore sees that the mass-squared of one of the two Higgs doublets, 
m^ 2 , at the EW scale is given by: 

2 1 197T 

m 2 H2 = m| 2 + fi%. z . = 2m 2 ~ (0.06)V3n s /^ 2 ™ 3 /2 + -M 2 EW - —So- (B37) 
From [9j, we notice: 

V^(i Zl z 2 m^/ 2 ~ m 2 ,, (B38) 
using which in (|B37p . one sees that for an 0(1) n s , 

m 2 H2 (M EW ) ~ lMm 2 + Im 2 ew -^S (B39) 
We have assumed at m 2 z , 2 {M s ) ~ m 2 , (implying 82 = but 5%^, 4 7^ 0). Further, 
m 2 Hl (M EW ) = (m 2 Zl + jl 2 ZlZ2 ) {Mew) ~ "4(^(1 + 5i) + -A/j^ + m 2 - (0.03) V^n s fi Zl z 2 m 3/2 

+(0.01)(n s ) 2 A| i22 ~ (1.97 + 8 x )m 2 + \m 2 ew + (0.01) (n 5 ) 2 ^ (B40) 

By assuming (/ifl) 2 ^ ~ frz z 2 ' see flB26[l - to be valid at the string and EW scales, the Higgs mass 
matrix at the £W-scale can thus be expressed as: 



pB m z H2 j \ Hi 2 m\ 2 



m Hi \ ( m 2 Hi ifi 2 



(B41) 



£ being an appropriately chosen 0(1) constant - see (|B45|h The eigenvalues are given by: 

\ (m 2 Hl + m 2 H2 ± ^J(m 2 Hi -m 2 Hi y + 4^) • (B42) 

As (for 0(1) n s ) 

"»Hi + m H 2 ~ (3-91 + ^m 2 , - O.O6S0 + ... 

" ™ir 2 ~ (0-03 + 5i)mg + 0.06,So + -, 
Al^a ~ 0.97m 2 , - 0.03S + ... (B43) 
one sees that the eigenvalues are: 



(3.91 + 5i)ml - O.O6S0 + .... ±y ((0.03 + d^m 2 + 0.065 + ...) 2 + £ 2 (l.94m 2 - 0.06S ) 2 . 

(B44) 

Hence, assuming non universality w.r.t. to both D3-brane position moduli masses (mz 12 ) an d 
squark/slepton masses, if So and £ are fine tuned as follows: 

(0.03 + Sx)m 2 ~ -O.O6S0 and £ ~ 2 + IH&t (B45) 

8 rriQ 

one sees that one obtains one light Higgs doublet (corresponding to the negative sign of the square 
root) and one heavy Higgs doublet (corresponding to the positive sign of the square root). Note, 
however, the squared Higgsino mass parameter fi Zl z 2 then turns out to be heavy with a value, at the 
EW scale of around 0.01Vm3/2 i.e to the order of squark/slepton mass squared scale which is possible 
in case of /i split SUSY scenario discussed above. This shows the possibility of realizing \i split SUSY 
scenario in the context of LVS phenomenology named as large volume "//-split SUSY" scenario. 
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